AD#  445 1:2 


l|^VV.  :jv  - 

;*•».*  *  i  '  *  .. 

iv  -VO  »  '  *  .'  >i.  • 

utvfr*  *■  ’C .  (  J  ,  '  ,  V  ''  i  ’ 

S^*  *  ■  i'  '*>.  *’  '  v”. 

iy~  Z  ~  ~i  ’*«  Sv  •  ,  , 

\  ''  >  * -  .  \  »  .’*  -  ■■/*  __ 

$*•*,  /  y,  r  *  V'  l  V  v 

pj  £  £<■  '  '  o  ' 

'■  ' 

£  1  s 

f'  '  ,-  *  ,iAo **Y*  '.*.***'  ‘  ^  *, 


.  •"••••  . 


♦***£•% 


|!?V*  "v  ^  n 


Ij&'cv,  ', ': 

V 


M,  • 
$ 
v 

jSo : 

•//  * 


\x 
i': ' 
&•& 

lr’ 


m  ANALYSIS 


m 


Hardcopy 

$  $vQO 


PANEL 


N  tt  H  OUSE 

»  *  ^ntific  and 

^'iVjpfeplf^  '  ■*  .MATION 


Microfiche 

—  ■  1  1  "”M‘ 

133 

i  pp 

f/ 

Everett  L  Cook 


-3j 


Department  of  Aeronautical  Engineering 

School  of  Engineering 
Wichita  State  University 
October  1966 


NOV  3  01966 

LiUil^LbU  U  LhlHi 


n 


Distribution  of  this  document  is  unlimited. 

Tins  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  ot  the  Army  position,  unless  so  designated  ’*v  *,nr 
authorized  ts. 


AN  ANALYSIS  OF  AN  INTEGRALLY  STIFFENED  PANEL 


by 

Everett  L.  Cook 


Aeronautical  Report  66-6 


Final  Report  of  Research  Performed  Under 
U.S.  Army  Research  Office  -  Durham 
Grant  No.  DA-ARO (D) -31-124-C482 


Department  of  Aeronautical  Engineering 
School  of  Engineering 
Wichita  State  University 
October  1966 


Distribution  of  this  document  Is  unlimited. 

_  in  this  renort  are  not  to  be  construed  as  an  official 

Department*  of  the  Army  position,  unless  so  desisted  by  other 
authorized  documents. 


INTRODUCTION 


Matrix  methods  are  used  extensively  in  the  analysis  and 
design  of  aircraft  structures.  Although  methods  have  been 
developed  for  both  isothermal  and  thermal  analyses,  experimental 
verification  of  these  methods  for  thin-wail  structures  is  not 
available.  The  purpose  of  the  investigation  described  in  this 
report  was: 

1.  To  design  and  test  on  integrally  stiffened  rectangular 
shear  panel  (Figure  1)  , 

2.  To  develop  analytical  procedures  for  shear  panels 
subjected  to  thermo  mechanical  loading,  and 

3.  To  compare  the  analytical  and  experimental  results  to 
evaluate  the  accuracy  of  the  analysis  methods. 

The  panel  has  been  designed  and  tested  under  mechanical 
loading.  There  is  good  agreement  between  the  analytical  and 
experimental  results.  It  is  also  shown  that  any  of  the  matrix 
methods  of  structural  analysis  investigated  have  the  same  order 
of  accuracy  when  they  are  applied  to  the  same  mathematical  model. 

The  body  of  the  report  presents  a  summary  description  of 
the  methods  of  analysis  which  were  investigated,  the  design  of 
the  test  panel,  the  experimental  program,  and  the  analytical 
and  experimental  results.  The  details  of  the  analysis  methods 
and  the  digital  computer  programs  developed  to  implement  the 
use  of  these  methods  on  the  test  panel  are  given  in  the  appendices. 

The  thermal  investigation  has  not  been  completed.  The 
progress  on  the  various  phases  of  this  investigation  is  presented 
in  the  appropriate  sections  of  the  body  of  the  report.  Appendix 
D  presents  the  significant  analytical  results  to  date.  This 
investigation  is  being  continued  as  a  thesis  project  by  Mr. 

R.G.  Merritt. 
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the  analytical  results  which  are  compared  with  the  experimental 
data.  This  program  does  not  at  present  include  provisions  for 
a  thermal  analysis. 

Both  the  force  analysis  and  the  displacement  were  performed 
with  the  Force/Displacement  Program  described  in  Appendix  B. 

It  is  also  shown  in  this  appendix  that  not  only  do  the  force 
and  displacement  methods  yield  essentially  identical  results, 
but  these  results  agree  with  those  given  by  the  direct  stiffness 
method.  This  excellent  agreement  is,  of  course,  directly  related 
to  the  fact  that  identical  mathematical  models  of  the  test  panel 
were  used  for  all  three  analyses.  This  means  that  the  accuracy 
of  the  results  of  a  matrix  structural  analysis  depend  only  on 
the  accuracy  of  the  mathematical  model  and  not  on  the  method 
of  analysis.  In  fact,  Przemieniecki  (Reference  6)  suggests 
that  if  the  substructure  method  is  used,  different  parts  of  a 
structure  may  best  be  analyzed  using  different  methods.  A 
force  method  for  calculating  the  displacements  and  internal 
forces  due  to  initial  strains  is  given  in  Appendix  C.  This 
will  be  expanded  to  include  the  displacement  method  for  the 
thermal  analysis. 


TEST  PANEL 

The  test  panel  is  of  monolithic  construction  and  is 
symmetrical  with  respect  to  the  plane  of  the  web.  Monolithic 
construction  was  selected  because  it  eliminates  both  joint 
friction  and  joint  thermal  resistance;  neither  of  which  can 
be  accurately  accounted  for  in  the  analyses.  Although  most 
actual  thin-wall  structures  are  not  symmetrical  about  the 
web  plane,  the  analysis  methods  assume  such  symmetry.  Therefore, 
this  construction  is  ideally  suited  for  obtaining  experimental 
data  for  comparison  wiuii  uie  anaxytical  results. 

The  geometry  of  the  test  panel  was  dictated  primarily  by 
the  results  on  the  heat  transfer  study.  Mr.  D.E.  Hull  conducted 
this  study  and  presented  a  paper  on  the  results  as  a  partial 
requirement  for  his  undergraduate  degree.  Mr.  Hull's  report 
is  given  in  Appendix  D  of  this  report.  The  purpose  of  this 
study  was  to  obtain  a  geometry  that  would  have  nearly  constant 
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temperatures  over  the  area  of  each  stringer  and  over  the 
length  of  each  panel  and,  at  the  same  time,  a  relatively 
large  difference  in  temperature  between  the  outboard  stringer 
and  the  centerline  of  the  panel.  The  geometry  selected 
satisfies  these  requirements. 

The  only  geometric  considerations  which  influenced  the 
design  for  mechanical  loading  were  the  overall  length  of  the 
panel  and  the  lug  design.  The  panel  length  was  limited  to 
approximately  twenty-four  inches  by  the  open  height  of  the 
testing  machine.  The  lugs  were  designed  for  a  limit  load  of 
13,000  pounds.  However,  the  applied  loads  were  limited  to 
approximately  7,500  pounds,  giving  a  factor  of  safety  of  2.4. 

The  material  selected  for  the  test  panel  was  6061-T6 
aluminum  alloy  plate.  Aluminum  alloy  tool  and  jig  plate  was 
investigated  and  found  to  have  not  only  low  physical  properties, 
but  also  a  nonlinear  stress-strain  curve. 

The  panel  was  machined  in  the  Wichita  State  University 
Engineering  Machine  Shop.  The  largest  available  mill  was  not 
large  enough  to  allow  milling  one  entire  side  without  moving 
the  pane.1  .  It  was  also  necessary  to  fill  the  panels  with 
water  during  the  final  cuts  on  the  web  to  prevent  thermal  buck¬ 
ling,  However,  the  desired  thickness  of  the  web  was  maintained 
within  a  tolerance  of  plus  0.0015-  and  minus  0.0005-inches, 
and  the  overall  thickness  was  within  plus  0.001-  and  minus 
0.000-inches  of  the  values  specified  in  Figure  1. 

TEST  FIXTURES 

The  mechanical  loading  was  performed  in  a  160,000  pound 
Tinus-Olsen  Universal  Testing  Machine.  The  test  panel  is 
shown  mounted  in  this  machine  in  Figure  2.  During  the  preliminary 
tests,  it  was  discovered  that  a  constant  load  could  not  be 
maintained  with  the  loading  system  of  the  testing  machine.  It 
was  then  necessary  to  design  a  seperate  loading  system  for 
the  final  tests.  Two  Blackhawk  hydraulic  cylinders  with  a  total 
capacity  of  28,000  pounds  were  mounted  on  the  upper  crosshead  of 
the  testing  machine  (Figure  3) .  A  linkage  was  designed  to 
transmit  the  load  from  the  hydraulic  cylinders,  through  the 
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opening  in  the  crosshead,  to  the  whiffle-tree  system.  The 
load  was  applied  and  maintained  with  a  hand-operated  hydraulic 
pump , 

The  whiffle-tree  systems  are  both  shown  in  Figure  2  and 
the  upper  system  is  also  shown  in  Figure  3.  The  upper  system 
was  designed  to  apply  equal  loads  at  all  four  of  the  lugs  along 
the  top  of  the  panel.  Its  compact  design  was  necessary  because 
of  the  limited  vertical  clearance  of  the  testing  machine.  Loads 
are  applied  only  at  the  inboard  stringers  along  the  bottom  of 
the  panel. 

The  thermal  tests  will  not  be  performed  in  the  testing 
machine,  but  in  a  special  fixture  being  fabricated  for  these 
tests.  Two  Conrow  SNH  scrip  heaters  will  be  used  to  radiate 
heat  to  the  edges  of  the  outboard  stringers.  All  other  stringer 
and  rib  faces  will  be  insulated.  The  web  panels  will  be  cooled 
by  forced  convection  to  produce  the  desired  temperature  distribu¬ 
tion  (Appendix  D) . 


STRAIN  GAGES 

The  strain  gages  used  on  the  panel  are  Budd  mediun. 
temperature  (MT)  foil  strain  gages.  The  axial  gages  on  the 
stringers  are  Type  C12-124-A  and  the  rosettes  are  Type  C12- 
124-R3C.  These  gages  have  a  one-eight  inch  gage  length  and  are 
recommended  for  temperatures  up  to  400°F  for  static  loading 
conditions.  The  strain  gage  mounting  pattern  is  shown  in  Figure 
1. 

The  cement  finally  selected  for  bounding  the  strain  gages 
was  Bean  BAP-1.  Budd  GA-50  .  as  used  for  the  gages  installed 
prior  to  the  preliminary  tests.  A  previous  experience  had  con¬ 
vinced  the  author  that  both  the  strain  gages  and  the  cement 
should  be  purchased  from  the  same  manufacturer.  It  was  discovered 
subsequently  to  the  preliminary  tests  that  the  cured  cement 
contained  microscopic  bubbles.  The  manufacturer's  sales 
representative  was  consulted  and  modifications  to  the  installation 
procedure  were  made.  When  these  modifications  did  not  eliminate 
the  bubbles,  a  new  batch  of  cement  was  obtained.  The  results 
were  the  same.  A  technical  representative  of  the  Budd  Company 
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was  finally  consulted  and  it  was  learned  that  the  GA-50  cement 
had  been  discontinued .  The  decision  was  then  made  to  use  the 
BAP-1  cement. 


INSTRUMENTATION 

A  Budd-Datran  20-channel  Automatic  Strain  Indicating  System 
was  originally  proposed  for  strain  gage  data  collection.  A 
Clary  digital  printer  owned  by  the  University  was  to  be  used  as 
the  output  device  for  the  system.  However,  this  printer  could 
not  be  used  without  major  modifications.  This  would  have 
increased  the  total  cost  of  the  system  to  an  amount  well  in  excess 
of  the  budgeted  amount.  As  a  result,  the  decision  was  made  to 
go  to  a  manual  system. 

The  manual  system  selected  consisted  of  2  Budd  Model  P-350 
Portable  Strain  Indicators  and  13  Budd  Model  SB-1  10-channel 
Switch  and  Balance  Units.  This  system  was  purchased  and  was 
used  in  the  preliminary  tests.  Subsequently,  Mr.  D.E.  Hull  and 
Mr.  Marvin  Davidson,  our  wind  tunnel  director,  modified  one  of 
the  P-350  Strain  Indicators  so  that  it  could  be  coupled,  through 
a  Dayton  Instruments  Model  DI-6-54  digital  converter,  to  an 
IBM  526  Printing  Summary  Punch.  Two  views  of  the  system  are 
shown  in  Figures  4  and  5.  Manual  switching  and  balancing  is 
still  necessary;  however,  the  strain  indicator  readings  are 
punched  into  cards-ten  per  card. 

Each  of  the  switch  and  balance  units  is  connected  to  nine 
gage  eiements--the  nine  axial  gages  on  one  side  of  a  stringer  or 
three  vertically  aligned  rosettes  on  a  web  panel.  Due  to  the 
large  number  of  gages,  the  load  must  be  applied  and  released 
four  times  for  a  complete  set  of  gage  readings.  To  facilitate 
rapid  and  accurate  switching  between  runs,  each  switch  and 
balance  is  permanently  wired  to  an  Amphenol  No.  26-190-32 
connector  and  each  set  of  nine  gage  elements  is  wired  to  a 
mating  No.  26-159-32  connector. 

It  was  originally  proposed  that  the  temperatures  would  be 
measured  with  standard  thermocouples  bonded  to  the  test  panel. 

The  possibility  of  measuring  the  temperatures  with  a  surface 

gage  are  now  being  investigated  in  order  to  avoid  adding  additional 
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wires  to  the  panel. 


EXPERIMENTAL  RESULTS 

The  mechanical  loads  were  applied  as  shown  in  Figure  2. 

Four  equal  loads  wore  applied  at  the  nodes  at  the  top  of  the 
panel.  These  were  reacted  by  two  equal  loads  at  the  inboard 
stringers  at  the  bottom  of  the  panel.  The  loads  were  applied 
in  increments  of  3,000-4,000  pounds  up  to  a  maximum  load  of 
approximately  15,000  pounds.  The  load  was  then  decreased  using 
approximately  the  same  increments.  Due  to  the  difficulty  in 
applying  the  load  uniformly  with  the  hand  pump,  no  attempt  was 
made  to  use  precise  increments;  however,  it  was  possible  to 
maintain  a  constant  load  within  plus  or  minus  100  pounds  while 
a  set  of  nine  strain  gages  was  being  balanced  and  recorded. 

The  experimental  results  are  shown  in  Figures  6-8,  along 
with  the  analytical  results  of  the  direct  stiffness  analysis 
(Appendix  A) .  The  experimental  points  in  these  figures  represent 
the  strains  in  the  vertical  direction-- e  in  Tables  A-3  through 

yy 

A-5 .  Strains  are  plotted,  rather  than  stresses  for  two  reasons: 

1.  The  stresses  at  the  location  of  a  strain  rosette 
cannot  be  calculated  if  any  of  the  elements  in  the 
rosette  are  inoperative. 

2.  Since  strain  readings  normally  have  only  two  or  three 
significant  figures,  while  at  least  eight  significant 
figures  are  carried  in  the  analysis,  it  is  more 
accurate  to  manipulate  the  analytical  results. 

A  large  number  of  the  strain  gages  were  damaged  beyond  repair 
during  the  installation  of  the  lead  wires.  The  stringer  gages 
were  particularly  easy  to  damage  although  brackets  were  used 
to  support  the  panel  while  the  lead  wires  were  being  installed. 

As  a  result  the  number  of  experimental  points  in  the  figures 
is  fewer  than  is  desirable.  The  solid  circles  represent  the 
average  strain  for  back-to-back  gages,  while  the  open  cirles 
represent  the  strain  at  points  where  one  of  the  back-to-back 
gages  was  inoperative. 

Each  of  the  plots  in  Figures  6-8  represents  the  strains 
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along  a  horizontal  line  of  strain  gages.  The  locations  are 
specified  in  the  sketches  on  the  figures.  The  ordinates  in 
the  figures  are  the  microstrain  per  pound  of  load  applied  to 
each  inboard  stringer.  The  analytical  results  (Tables  A-2 
through  A-5)  are  tabulated  in  essentially  this  form.  The 
experimental  results  were  reduced  to  this  form  by  using  the 
slope  of  a  least-squares  straight  line  through  the  data. 

There  is  generally  good  agreement  between  the  analytical 
and  experimental  results.  The  only  explainable  differences  are 
those  for  the  line  of  gages  nearest  the  bottom  of  the  panel. 

It  appears  that  the  loads  are  transfered  from  the  stringers 
to  the  web  more  rapidly  than  the  analytical  solution  can  predict 
with  the  subelement  size  used  here.  It  would  be  better,  therefore, 
to  have  smaller  subelements  in  this  region.  This  is  not  possible 
with  the  present  computer  configuration. 

CONCLUSIONS 

The  following  conclusions  can  be  drawn  from  the  results  of 
this  investigation: 

1.  The  mathematical  model  of  the  test  panel  used  for  the 
direct  stiffness  analysis  was  sufficiently  accurate 
to  provide  good  agreement  between  the  analytical  and 
experimental  results. 

2.  The  force  and  displacement  methods  could  have  been 
used  with  equally  accurate  results. 

3.  The  derivation  of  the  recurrence  relationship  for 
initial  strains  extends  the  recurrence  method  of 
Pestel  and  Leckie  to  include  thermal  loading. 
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Figure  3.  Loading  System  and  Upper  Whi f f le-Tree 


Figure  4.  Strain  Gage  Data  Collection  System 
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Figure  5.  Modified  Strain  Indicator  and  Card  Punch 
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Figure  6.  Analytical  and  Experimental  Strains-Upper  Panels 
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MICROSTRAIN/POUND 


Figure  7.  Analytical  and  Experimental  Strains-Center  Panels 
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Figure  8.  Analytical  and  Experimental  Strains-Lower  Panels 


APPENDIX  A 


THE  DIRECT  STIFFNESS  ANALYSIS  PROGRAMS 

The  Direct  Stiffness  Analysis  Programs  ,  arbitrarily  named 
DORA,  were  developed  to  analyze  the  test  panel  using  the  method 
presented  by  Turner,  Clough,  Martin,  and  Topp  in  Reference  7. 
There  are  four  versions  of  the  program:  however,  only  the  version 
(DORA  IV)  used  to  obtain  the  results  to  be  compared  with  the 
f orce/displacement  results  (Appendix  B)  and  the  experimental 
data  will  be  described  in  detail.  The  only  inputs  required  are 
the  geometry  and  physical  properties  of  the  structure  and  the 
external  loads  and  constraints  on  the  structure.  The  outputs  of 
the  program  are  the  displacements  of  the  nodes.  An  auxiliary 
program  was  written  to  calculate  the  strains  in  the  stiffeners 
and  the  webs  from  the  nodal  displacements. 

The  direct  stiffness  method  of  Turner,  Clough,  Martin,  and 
Topp  is  the  simplest  of  all  the  methods  of  matrix  structural 
analysis.  It  is  based  on  the  following: 

1.  The  force-displacement  relationship  for  each 
element  in  a  structure  can  be  expressed  as 

{ F }  =  [K] {u}  ( A-l) 

where 

{ F }  =  the  column  matrix  of  the  nodal  forces 
on  the  element.  All  nodal  forces  are 
included.  In  the  displacement  analysis 
of  Appendix  B,  only  a  subset  of  independent 
forces  are  considered. 

{u}  =  the  absolute  nodal  displacements  of  the 
element.  In  the  displacement  analysis 
only  the  displacements  relative  to  a 
set  of  reference  axes  in  the  element  are 
considered . 

[K]  =  the  element  stiffness  matrix.  This  matrix, 
sometimes  called  the  unreduced  stiffness 
matrix,  is  singular.  The  reduced  stiffness 
matrix  used  in  a  displacement  analysis  is 
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derived  from  the  unreduced  stiffness 
matrix  by  the  elimination  of  specified 
rows  and  columns . 

2.  The  force-displacement  relationship  for  each  element 
in  a  structure,  Equation  (A-l) ,  can  be  expanded  to 
the  order  of  the  assembled  structure  stiffness  matrix 
by  adding  rows  and  columns  of  zeros  for  the  nodal 
displacements  that  are  irrelevant  for  the  element 
in  question.  These  expanded  element  stiffness  matrices 
can  then  be  added  to  obtain  the  stiffness  matrix  of 
the  assembled  structure.  The  resulting  force- 
displacement  relationship  is  of  exactly  the  same 
form  as  Equation  (A-l) ,  but  the  force  vector  re¬ 
presents  the  external  forces  at  the  nodes.  Martin 
(Reference  4,  Section  2.6)  has  an  excellent  discussion 
of  the  assembling  of  the  stiffness  matrix. 

The  method  of  expanding  each  element  stiffness  matrix  to  the 
order  the  structure  stiffness  matrix  is  not  practical  for  machine 
computation.  Even  if  only  one  element  was  considered  at  a  time, 
the  storage  required  would  be  cwice  that  required  for  the 
structure  stiffness  matrix.  It  is  relatively  easy,  however,  to 
determine  the  position  in  the  structure  stiffness  matrix  for 
each  element  in  the  element  stiffness  matrix.  Thus,  when  an 
element  stiffness  matrix  has  been  obtained,  its  elements  can  be 
added  into  the  proper  position  in  the  structure  stiffness  matrix. 
This  latter  method  is  used  in  all  versions  of  DORA. 

The  stiffness  matrix  of  the  assembled  structure  is  singular. 
It  must  be  made  nonsingular  so  that  an  inversion  may  be  per¬ 
formed  to  obtain  the  displacements  in  terms  of  the  external 
forces ,  i . e . , 

{ u }  =  [C] {F}  (A-2) 

where  [C]  =  the  flexibility  matrix  of  the  structure. 

There  are  two  practical  procedures  for  obtaining  the  flexibility 
matrix: 

1.  The  rows  and  columns  of  the  stiffness  matrix 

correspond: ng  to  the  external  constraints  on  the 
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structure  can  be  eliminated.  This  nonsingular 
reduced  matrix  can  then  be  inverted.  Finally, 
the  flexibility  matrix  can  be  obtained  by  ex¬ 
panding  back  to  the  original  order  and  inserting 
zeros  in  the  rows  and  columns  corresponding  to 
the  constraints. 

2.  All  but  the  diagonal  elements  in  the  rows  and 

columns  corresponding  to  the  constraints  can  be 
set  to  zero  and  the  diagonal  elements  can  be 
set  to  unity.  This  matrix  is  nonsingular  and 
can  be  inverted.  Then  the  ones  on  the  diagonal 
corresponding  to  the  constraints  can  be  set  to 
zero . 

The  second  procedure  is  used  in  DORA  primarily  because  it  is 
easier  to  program  for  arbitrarily  specified  constraints. 

Once  the  flexibility  matrix  has  been  determined,  the  nodal 
displacements  may  be  calculated  from  Equation  (A-2) .  These 
displacements  may  then  be  substituted  into  the  element  force- 
displacement  relationship,  Equation  (A-l) ,  to  obtain  the  nodal 
forces . 

Several  versions  of  the  program  were  necessary  because  of 
the  limited  size  of  the  available  digital  computer.  An  IBM 
1620  with  40,000  positions  of  core  storage  and  card  input- 
output  was  used  (all  w^ork  prior  to  January,  1966,  was  done  on 
the  same  machine  with  only  20,000  positions  of  core  storage). 

The  programs  are  all  written  in  PDQ  Fortran.  This  programming 
system  was  selected  because  it  requires  less  core  storage  for 
the  subroutines  than  do  any  of  the  IBM  Fortran  systems.  It  is 
also  faster  than  the  IBM  Fortrans,  but  this  is  a  secondary 
consideration.  In  the  near  future,  a  1311  disk  drive  will  be 
added  to  the  computing  system.  The  entire  set  of  programs  can 
then  be  combined  and  enlarged. 

As  mentioned  above  only  DORA  IV  will  be  described  in  detail, 
however,  the  other  versions  can  be  summarized  as  follows: 

DORA  I 

This  program  was  written  to  demonstrate  the  feasibility  of 
generating  the  assembled  structure  stiffness  matrix  as  outlined 
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above.  Only  structures  with  bar  and/or  beam  elements  can  be 
analyzed  with  this  program.  These  elements  were  chosen  because 
there  are  many  truss,  continuous  beam,  and  frame  examples  in 
the  literature  which  provided  results  for  checking  and  debugging 
the  program.  This  program  definitely  proved  the  feasibility 
of  the  method.  However,  in  order  to  make  it  a  one-part  program, 
the  structure  is  limited  to  fourteen  nodes. 

DORA  II 

The  program  was  broken  down  into  parts  and  triangular  and 
rectangular  plate  elements  were  added.  The  maximum  number  of 
nodes  is  limited  to  sixteen. 


DORA  III 

Trapezoidal  plate  elements  were  added  to  the  program  and 
the  beam  elements  were  eliminated.  The  elimination  of  the  beam 
elements  reduced  the  number  of  degrees  of  freedom  at  each  node 
from  three  to  two.  However,  thv_  part  of  the  program  which 
generated  the  flexibilities  of  the  trapezoidal  plate  elements 
was  very  large,  so  that  the  maximum  number  of  nodes  is  limited 
to  twenty. 


DORA  IV 

The  trapezoidal  plate  element  was  eliminated  and  the  maximum 
number  of  nodes  was  increased  to  twenty-five.  This  is  large 
enough  for  the  analysis  of  the  test  panel.  Since  this  is  the 
version  used  for  the  analysis  of  the  test  panel,  it  will  be 
described  in  more  detail.  The  program  is  divided  into  six  parts — 
Initialization,  Bars,  Triangular  Plates,  Rectangular  Plates,  and 
Matrix  Inversion  and  Output.  A  program  listing  is  given  in 
Table  A-l .  The  purpose  and  use  of  each  part  is  as  follows: 


Initial ization.  The  number  of  nodes,  Young's  Modulus, 
Poisson's  Ratio,  and  the  x-  and  y-coordinates  of  the  nodes  are 
read  and  stored.  They  are  also  repunched,  with  appropriate 
headings,  for  identification  and  debugging  purposes.  The  size 
of  the  structure  stiffness  matrix  is  calculated  and  the  elements 
of  this  matrix  are  set  to  zero. 

Bars .  The  data  cards  for  all  of  the  bar  elements  in  the 
structure  are  the  input  to  this  part  of  the  program.  If  the 


A-4 


structure  does  not  contain  bar  elements,  this  part  of  the 
program  is  not  loaded.  Each  bar  is  described  on  one  card.  The 
numbers  of  the  two  nodes  which  the  bar  connects,  the  area  of 
the  bar,  and  a  control  code  must  be  specified.  The  control  code 
for  each  card  except  the  last  is  a  zero.  A  one  in  the  last  card 
causes  the  next  program  part  to  be  read  after  the  last  element 
has  been  processed. 

The  bar  stiffness  matrix  was  taken  from  Reference  7.  The 
order  has  been  changed  so  that  the  x-  and  y-displacements  of 
the  first  node  preceed  those  of  the  seconl  node.  There  are 
no  restrictions  on  the  numerical  order  of  the  nodes,  i.e.,  the 
number  of  the  first  node  may  be  higher  or  lower  than  the  number 
of  the  second  node. 

Triangular  Plates .  The  input  to  the  triangular  plate  pro¬ 
gram  consists  of  cards  containing  the  three  node  numbers,  the 
thickness,  and  the  control  code  (See  "Bars") .  This  stiffness 
matrix  was  also  taken  from  Reference  7.  Rearrangement  was  not 
necessary.  The  r ades  may  be  specified  in  any  order. 

Rectangular  Plates .  The  numbers  of  the  four  nodes,  the 
thickness,  and  the  control  code  (See  "Bars")  must  be  specified 
for  each  rectangular  plate.  The  order  of  the  nodes  is  important! 

The  required  order  is  lower  left,  lower  right,  upper  left,  and 
upper  right. 

The  rectangular  plate  stiffness  matrix  was  taken  from 
Reference  1,  P.  50.  The  shear  and  direct  stress  stiffness 
matrices  were  combined  and  rearranged  in  the  node  order  specified 
above.  Argyris'  stiffness  matrix  was  chosen  over  the  one  given 
in  Reference  7  because  it  is  compatible  with  the  trapezoidal 
plate  stiffness  matrix  used  in  DORA  III.  Belirgen  (Reference  2) 
has  shown  that  either  of  these  stiffness  matrices  gives  satisfactory 
results  for  the  type  of  structure  being  considered. 

Matrix  Inverse  and  Output .  The  primary  purpose  of  this  part 
of  the  program  is  to  invert  the  stiffness  matrix  of  the  assembled 
structure.  However,  this  matrix  must  first  be  made  nonsingular. 

The  constraints  are  specified  with  a  data  card  for  each  constrained 
node.  The  data  cards  contain  the  node  number;  the  x-  and  y- 
constraint  codes,  and  the  control  code  (See  "Bars").  The  con¬ 
straint  code  is  1.0  for  a  constrained  displacement  and  0.0  for 
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an  unconstrained  displacement.  After  each  card  is  read,  the 
row(s)  and  column{s)  in  the  stiffness  method  corresponding  to 
the  constrained  displacement ( s)  are  modified  by  the  method 
outlined  previously.  When  all  constraints  have  been  imposed, 
the  stiffness  matrix  is  inverted  to  obtain  the  flexibility 
matrix . 

The  loads  on  the  structure  are  specified  similarly  to  the 
constraints.  A  card  is  read  for  each  node  with  at  least  one 
external  load.  The  card  contains  the  node  number,  the  x-  and 
y-loads,  and  the  control  code  (See  "Bars") .  Both  loads  must 
be  specified  although  one  of  them  may  be  zero.  The  loads  are 
placed  in  the  load  vector  [F] ;  and  after  the  last  card  is  read, 
the  displacements  are  calculated  using  Equation  (A-2) .  The 
displacements  are  then  punched  into  cards.  The  x-  and  y- 
displacements  (u  and  v)  of  each  node,  including  the  constrained 
nodes,  are  punched. 

The  stiffness  matrix  is  replaced  by  the  flexibility  matrix 
during  the  inversion;  therefore,  the  stiffness  matrix  must  be 
reassembled  in  order  to  study  the  effects  of  different  external 
constraints  on  a  structure.  This  will  not  be  necessary  when 
the  disk  drive  becomes  available  since  the  stiffness  matrix 
may  be  stored  on  the  disk  and  recalled  when  needed.  The  effects 
of  any  number  of  loading  conditions  can  be  studied  by  generating 
a  new  load  vector,  i.e.,  by  reading  in  a  new  set  of  load  cards. 

Examples 

Four  panel  configurations  were  analyzed  using  DORA  IV 
(Figure  A-l) .  Since  the  test  panel  was  loaded  symmetrically, 
only  one  half  of  the  panel  is  analyzed.  The  roller  constraints 
along  the  centerline  of  the  panel  maintain  the  symmetry  of 
lateral  displacement,  but  allow  longitudinal  displacements.  The 
constraint  on  the  lower  edge  represents  the  load  applied  with 
the  lower  whiffle  tree. 

Panel  Configuration  I  is  also  analyzed  by  the  force  and 
displacement  methods  and  the  results  are  compared  (Appendix  B) . 
The  results  are  given  in  Table  A-2.  The  other  three  configura¬ 
tions  could  nor  be  analyzed  with  the  Force/Displacement  Program 
due  to  the  large  number  of  internal  forces  and  deformations. 


Panel  Configurations  II-IV  were  analyzed  to  obtain 
stiffener  and  web  strains  for  comparison  with  the  experimental 
data.  Due  to  the  size  limitations  of  the  computer,  it  was  not 
possible  to  use  one  configuration  with  small  enough  subelements 
to  provide  strains  which  would  compare  favorably  with  the 
experimental  data.  Therefore,  Configurations  II,  III,  and  IV 
were  used  to  obtain  strains  in  the  upper,  center,  and  lower 
panels,  respectively.  The  vertical  dimensions  fo'-  the  smaller 
subelements  (Nodes  1-16)  were  selected  so  that  the  line  of 
strain  gages  was  along  the  horizontal  centerline  of  the  elements. 
The  results  are  shown  in  Tables  A-3  through  A-5. 

The  present  versions  of  DORA  do  not  have  provisions  for 
calculating  stresses  or  strains.  A  special  program  was  written 
to  calculate  the  axial  strains  in  the  stiffener  elements  and 
the  direct  and  shear  strains  at  the  center  of  the  web  elements. 

A  listing  of  this  program  is  given  in  Table  A-6. 

The  method  used  here  for  obtaining  the  strains  for  comparison 
with  the  experimental  data  yields  satisfactory  results  and 
could  undoubtedly  be  used  successfully  with  most  other  types 
of  structures.  However,  a  substructure  method,  such  as  the  one 
presented  in  Reference  6,  would  certainly  yield  more  accurate 
results.  The  capability  has  not,  as  yet,  been  incorporated  into 
DORA. 
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C  C 


TABLE  A - i 


DORA  PROGRAM  LISTING 

C  DORA  IV  -  BARS  ANO  PLATES  (EXCEPT  TRAPEZOIDAL) 

C 

C  INITIALIZATION  -  25  NODES  -  2  DEGREES  OF  FREEDOM 

C  *************  **********************  **************** 

c 

C  E i L »  COOK  -  MAY  1966 

C 

COMMON  NO,  NOX lEt V,STIFF<50*50)*X<25)«Y(25> 

PUNCH  AND  PRINT  DORA  IV  HEADING 

CONTROL  102 
PUNCH  944 
PUNCH  999 
PRINT  944 

PRINT  PROGRAM  SWITCH  2  OFF  MESSAGE 

if  sense  switch  2  is  *on*,  the  messages  at  the  end  of  the  program 
are  BYPASSED. 

CONTROL  102 
PRINT  950 

PUNCH  NODES*  MODULUS ,  AND  RATIO  HEADINGS 

PUNCH  999 
PUNCH  955 
PUNCH  956 
PUNCH  999 

READ  NUMBER  OF  NODES,  YOUNG*S  MODULUS,  AND  POISSON*S  RATIO 

10  READ  901 , NO , E , V 
PUNCH  90 1 , NO, E , V 

PUNCH  COORDINATE  HEADING 

PUNCH  999 
PUNCH  999 
PUNCH  960 
PUNCH  999 

READ  NODE  NUMBER,  X  COORDINATE*  AND  Y  COORDINATE. 

DO  20  J= I *NO 
READ  902, I  .X<  I >  « Y  <  I  > 

20  PUNCH  902  *  I»X(  I  >,Y(  I  ) 

C 

NOX  =  2*NO 
DO  25  I  =  I  *  NOX 
DO  25  J= I  *  NOX 
25  STIFF! I , J>=0.0 
C 

AL INK=L I NK ( 1 .0 ) 
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c 


TABLE  A-l  (CONTINUED) 


901  FORMAT* I3.2E16.8) 

902  FORMAT* 13.2E16.8) 

944  FORMAT (42H  DORA  1 V  -  25  NODES  -  2  DEGREES  OF  FREEDOM) 
950  FORMAT (26H  TURN  PROGRAM  SWITCH  2  OFF) 

955  FORMAT ( 5HN0  OF.3X.8H  YOUNG* S . 7X .  1  OH  POISSCN*S) 

956  FORMA  T  <  5HN0DE  S , 3x  «  8H  MODULUS . 9X . 6H  RATIO) 

960  FORMA T ( 5H  N0DE.6X.2H  X.14X.2H  Y) 

999  FORMAT ( 1H  ) 


o  o  r>  n  ci  ci  o  ci  Ci  ci  n  Ci  Ci 


TABLE  A~ I  (CONTINUED) 


C  DORA  IV  -  BARS  AND  PLATES  (EXCEPT  TRAPEZOIDAL) 

BARS  -  25  NODES  -  2  DEGREES  OF  FREEDOM 

***********>.***************************** 

E.L.  COOK  -  MAY  1966 

COMMON  NO'.  NOX  (E'ViST  IFF  <  50 , 50  > . X ( 25  >  « Y ( 25  > 

DIMENSION  ESTIF(4.4> 

PUNCH  HEADINGS 

PUNCH  999 
PUNCH  999 
PUNCH  901 
PUNCH  999 
PUNCH  902 
PUNCH  999 

READ  THE  NUMBERS  OF  THE  TWO  END  NODES.  THE  AREA,  AND  THE  CONTROL 
CODE.  THE  CONTROL  CODE  IS  O  FO  ALL  MEMBERS  EXCEPT  The  LAST  ONE. 
FOR  WHICH  IT  IS  1 . 

C 

30  READ  903, I , J, AR.kk 
PUNCH  903. I . J. AR 
C 

XJI=X( J)-X< I ) 

YJI=Y(J)-Y< I ) 

Sp I = sqrtf (XJI*XJI+YJI*YJI  ) 

C  A  =  LAMBDA  AND  U  =  MU 

A=XJ I/SJI 
ASQ=A*A 
U=YJ I/SJI 
USQ=U*U 
AU=A*U 
AEL  =  AR*E/S J I 
EST I F (  5  » I  ) =AEL*ASQ 
EST I F ( 1 ,2) =AEL*AU 
ESTIFf 1 »  3 ) =— EST I F (  1,1) 

EST IF ( 1  ,4 ) =-EST IF (  I , 2  > 

EST IF (2 ,2 ) =AEL*USQ 
ESTIF(2,3) =ESTIF( I ,4) 

ESTIF(2«4)=-ESTIF<2.2) 

ESTIF(3«3)=ESTIF(  I  .  1  ) 

ESTIF (3,4) iESTIF ( 1 ,2) 

ESTIF(4,4) =EST!F(2.2  > 

DO  40  K= 1 ,4 
DO  40  L=K,4 

a  n  ffcr  ic  ii_  .if  )  -ccttc;  ; 

DO  50  K= I , 2 
K  I  .-K  +  2*  (  I  -  I  ) 

K J  =  K  +  2* ( J—  I  ) 
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TABLE  A- I  (CONTINUED) 


C  C 

DO  50  L=  I  ,  2 
LI =L+2* ( I -I  ) 

LJ=L+2*( J-I ) 

STIFF(KI  ,LI  )=ST!FF(KI .LI  )+ESTIF(K,L) 
STIFF(KI,LJ)=STIFF(KI,LJ>+ESTIF(K,L+2> 
ST!FF(KJ»LI )  =  ST  I FF (KJ«LI )+ESTIF(K+2«L) 

50  STIFF(KJ,LJ)=STIFF(KJ,LJ>+E3TIF(K+2«L+2) 
IF(KK)890i30,60 
C 

60  ALI NK=L I NK ( I .0 > 

C 

Q90  STOP 
C 

901  FORMAT (4HBARS) 

902  FORMAT (6H  NODES, 6X*5H  AREA) 

903  FORMAT ( 2  I  3 ,E 16.8, I  3 ) 

999  FORMAT ( !H  ) 

C 

END 
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non  n  n  o 


C  C  TABLE  A— I  (CONTINUED) 

C  DORA  IV  -  BARS  AND  PLATES  (EXCEPT  TRAPEZOIDAL) 

C 

C  TRIANGULAR  PLATES  -  25  NODES  -  2  DEGREES  OF  FREEDOM 

C  **************************  ****  ##*#■*•***  *«■*■**#***#####■**# 

C 

C  E.L.  COOK  -  may  1966 

C 

COMMON  NO«NOX«E« V«STIFF<50«50) .X(25) * Y(25) 

DIMENSION  ES  T  I  F ( 6  «  6  ) 

C 

NDOF  =2 

PUNCH  HEADINGS 


PUNCH  999 
PUNCH  999 
PUNCH  901 
PUNCH  999 
PUNCH  902 
PUNCH  999 


read  three  node  numbers,  the  thickness,  and  the  control  code 

80  Read  906. I . J.K.T.KK 
PUNCH  906. I . J.K.T 


XJI=X( J)-X( I ) 

XKJ  =  X(K  )-X( J) 

X  I  K  =  X  (  I  )~X(K> 

YJI=Y( J)-Y( : ) 

YKJ-Y(K ) — Y ( J ) 

YIK=Y( I ) — Y ( K ) 

ET  =  0.5*E*T/( (  1 ,0-v**2)*ABSF(XJI*Y(K)+xiK*Y( J)+XKJ*Y(  I  )  )  ) 
VI  =0.5* ( I .0-V  > 

V2=0 • 5* < I .0+V> 


EST IF (  I 
EST I F (  I 
EST IF( I 
EST IF( l 
EST  IF(  I 
EST IF( 1 
EST I F ( 2 
EST IF(2 
EST IF( 2 
EST IF ( 2 
EST I F (2 
EST I F ( 3 
EST IF( 3 
E3T N  C  3 
EST I F ( 3 
EST 1F( A 
EST  IF( A 
EST I F ( A 
EST I F ( 5 
EST I F ( 5 
EST I F ( 6 


1 )  =V 1 *XK J*XK J+YKJ* YK J 

2)  =-V2*XKJ*YKJ 

3)  =V1*XKJ*XIK+YKJ*Y  IK 

A)  =  — V 1  *  XK  J* Y IK-V*XIK*YKU 

5)  =  V  I *X J I *XK J  +  Y J I *  YKU 

6 )  =  — V  I *XKJ*Y J I - V*  X J I *YK J 

2 )  =XKJ*XKJ+V I # YK J*  YK  J 

3)  =~V 1  *  X  I K*YKJ-V*XK J*Y IK 
A) =XKJ*XIK+Vl*YKJ*YIK 

5)  =-Vl*XJI*YKJ-V*XKJ*YJ! 

6) =xji*xkj+vi*yji*ykj 

3)=VI*XIK*XIK+YIK*YIK 

A)=~V2*XIK*vik 

f>  I  it./  i  r  r  1  x/  (  I  «  V  I  i/ 

«J/—  i  4  • ' 

6)=-VI*XIK*YJ!-V*-XJ1*YIK 

а) =xik*xik+vi#yik*yik 

5)  =-Vl*XJI*YIK-V*XIK*YJI 

б) =XJI*XIK+V1*YJI*YIK 
5>=Vl*XJI*XJI+YJI*YJI 

6 )  =-V2*X J I *Y J I 
6>=XJI*XJI+VI*YJI*YJI 
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C  C  TABLE  A- 1  (CONTINUED) 

DO  85  L= I .6 
DO  85  M=L «  6 

EST IF(L  «M> =ET*EST IF(L.M) 

85  £STIF(M«L)=ESTIF(L«M) 

DO  90  L  =  1  «  2 
L I =L+N00F* (  I  -  I  ) 

I _ l=L+NDOF* ( J-l ) 

LK=L+ND0F*(K-1  ) 

DO  90  M=  1  ♦  2 
MI =M+NDOF* (  I- I  ) 

MJ=M+NDOF* ( J- I ) 

MK=M+NDOF*(K-I ) 

STIFF (LI  «  M I  )  =ST  I FF ( L  5 «  Ml  )+ESTIF(L.M) 

ST IFF(LI «MJ)=ST IFF !L I .MJ)+EST IF (L.M+2 ) 

ST IFF(L I  «MK 1 -ST  I FF (L 1  «MK  >  +  £ST iF ( L  *M  +  4 ) 

ST  I FF ( L J  <  M I  )  =  ST  I FF (  L  J  «  M  I  )+£STIF(L+2*M) 
STIFF(LJ.MJ)=STIFF(LJ,MJ)+ESTIF(L+2.M+2) 
STIFF(LJ,MK)=STIFF(LJ,MK)+ESTIF(L+2»M+4> 
STIFF ( LK  «  M I  1 =  ST  I FF ( LK  «  M I  )+ESTIF(L+4,M) 
STIFF(LK.MJ) =STIFF(LK.MJ)+ESTIF(L+4«M+2) 

90  ST IFF(LK. MK 1 =ST IFF (LK.MK  >  +ESTIF ( L  +  4 . M  +  4 ) 

C 

IF(KK>890.80, 100 
C 

100  A  L I NK  =  L I NK (  I  • 0 ) 

C 

890  STOP 
C 

901  FORMAT!  I 7HTR I  ANGULAR  PLATES) 

902  FORMAT (7H  NODES.6X.I0H  THICKNESS) 

906  FORMAT(3I3.EI6.8. 13) 

999  FORMAT ( IH  ) 

C 

END 
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TABLE  A— 1  (CONTINUED) 


C  DORA  IV  -  BARS  AND  PLATES  (EXCEPT  TRAPEZOIDAL) 

C 

C  RECTANGULAR  PLATES  -  25  NODES  -  2  DEGREES  OF  FREEDOM 

C  ****************************  ***********  *****  #•*#**#**■*■»■*■* 

C 

C  E.L.  COOK  -  MAY  1966 

C 

COMMON  NO.  NOX<E<V<ST  1  FF  (  50 . 50  )  .  X  (  25  )  '.  Y  (  25  ) 

DIMENSION  EST I F ( 8 . 8  > 

C 

C  PUNCH  HEADINGS 

C 

PUNCH  999 
PUNCH  999 
PUNCH  901 
PUNCH  999 
PUNCH  902 
PUNCH  999 
C 

C  READ  THE  NUMBERS  OF  THE  FOUR  NODES.  THE  THICKNESS.  AND  THE  CONTROL 

C  CODE.  THE  NODE  NUMBERS  MUST  BE  IN  THE  FOLLOWING  ORDER  -  LOWER 

C  LEFT.  LOWER  RIGHT,  UPPER  LEFT.  AND  UPPER  RIGHT.  THE  CONTROL  CODE 

C  IS  0  FOR  ALL  CARDS  EXCEPT  THE  LAST  ONE.  FOR  WHICH  IT  IS  I. 

C 

60  READ  905. I . J.K.L, T.KK 
PUNCH  90S, I . J.K ,L . T 
C 

AM=(X(J)-X(1)>/(Y(K)-Y( I)) 

ET=E*T/(2A.0*( 1 • 0-V**2 )  ) 

A1 =ET*8.0/AM 
B1=ET*8.0*AM 
C1=ET*6.0*V 
A2=ET*4.0*AM* ( 1 . 0— V ) 

B2  =  ET  *4 .0* (  1 .0-V)/AM 
C2=ET*3.0* ( 1 .0 -V  ) 

EST IF (  I  ,  1  ) =  A  1+A2 

EST 1 F ( 1 ,2) =C1+C2 

EST IF ( 1,3) =-A 1+0.5* A2 

EST1F( 1 , 4 ) =C i — C2 

EST IF( 1 ,5) =0.5*A1-A2 

EST IF ( 1.6) =-EST IF<  1  , A ) 

EST IF (  1  »  7 ) =— 0 • 5*EST I F (  1.1  ) 

EST I F ( 1  » 8  > =-EST IF< 1,2) 

EST 1F( 2.2) =8  1+82 
ESTIF(2,3)=ESTIF(1 .6) 

EST 1F( 2, 4 ) =0 . 5*B 1 -82 
EST IF( 2.5) =ESTIF( 1 .4 ) 

EST IF ( 2.6 >  =-B 1  +  0.5*52 

ps+IP'B-^^ESTIP’'  1  .Hi 
ESTIF(2.8)=-0.5*ESTIF(2,2) 

EST I F ( 3 . 3 ) =EST I F ( 1.1) 

EST IF( 3.4) =EST1F ( 1 .8) 

EST IF ( 3,5) =EST IF ( 1 ,7 ) 

EST I F ( 3 , 6 ) =ESTIF( 1.2) 

EST IF( 3,7) =EST IF ( 1 ,5 ) 

ESTIF(3.8)=ESTIF( 1,4) 
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TABLE  A—  1 


(CONTINUED) 


C 

c 


EST IF(4.4)rESTIF(2.2) 

ESTIF(4,5)=ESTIF( I ,2) 
ESTIF(4»6)=ESTIF(2,8) 

ESTIF(4,7)=ESTIF( I ,6) 
EST1F(4»8)=£STIF(2«6) 

ESTIF^  5  )  =E£T  I  F  (  )i  I  ! 

ESTIF(5,6>=ESTIF( I .8) 

EST IF(5,7) =EST IF( i ,3) 

ESTIF(5.8)-ESTIF( I .6) 
EST1F(6,6>=EST!F(2,2) 

ESTIF(6.7)=ESTIF( I .4) 

EST IF ( 6. 8) =E ST  IF (2.4) 

ESTIF(7.7)=ESTIF( 1.1) 

EST IF (7,8) =  E ST  I F ( I ,2) 
ESTIF(8.8)=ESTIF(2.2) 

DO  65  M= I ,8 
DO  65  N  =  M , 8 

65  EST IF (N.M) =EST IF( M.N) 

DO  70  Mr i,2 
MI rM+2*  <  I-  I  ) 

MJrM+2*( J-I ) 

MKrM+2* (K— I ) 

ML=M+2*(L-1 ) 

DO  70  N=I ,2 
Nl  rN+2*  <  I -  I  ) 

NJ=N+2*! J-l  ) 

NKrN+2*<K~t > 

NLrN  +  2*  (l.-  I  ) 

ST  IFF (MI  ,NI  )r ST  IFF (MI ,NI  )+ESTIF(M.N) 
STIFF (MI ,NJ)rSTlFF(MI , N J ) +  EST I F ( M , N+2 ) 

ST  IFF ( MI  ,NK ) r ST  IFF ( MI  .NK ) +EST IF ( M,N+4 ) 
STIFF ( MI ,NL ) =ST IFF(MI . NL > +EST I F ( M . N+6 > 

ST IFF( MJ.NI  ) r ST  IFF (MJ ,Ni  ) +EST IF ( M+2.N) 

ST  IFF(  Mvl.NJl  rST  IFF  (MJ.NJ)  +ESTIF  (  M+2.N+2  ) 
STIFF ( MJ, NK ) rST IFF (MJ, NK) +EST IF ( M+2.N+4 ) 
STIFF( MJ.NL )  =  ST IFF (MJ.NL  >  +EST IF ( M+2.N+6 ) 
STI FF( MK.NI  ) r ST  IFF ( MK ,NI  ) +EST IF (M  +  4 , N) 
STIFF (MK ,NJ) r ST IFF( MK.NJ) +EST IF ( M+4 ,N+2) 
ST  IFF ( MK.NK ) =ST !FF( MK.NK) +EST IF  <  M+4.N+4 ) 
ST  IFF  (MK  ,NL  )  =ST  >  FF  <  MK<:NL)  +  EST  IF  (  M+0  ,  N  +  6  ) 
STIFF(ML,NI)rSTlFF(ML,NI)+ESTIF(M*6,N) 

ST IFF(ML.NJ) rST IF  F ( ML, NJ) +EST IF ( M+6 . N+2 ) 
ST IFF( ML.NK ) rST IFF(ML.NK) +ESTIF ( M+6 ,N+4 ) 

70  ST  IFF ( ML.NL ) =ST IFF (ML, NL) +ESTIF (M  +  6 , N  +  6 ) 

IF (KK) 890 , 60 , 100 
100  AL I NK=L I NK ( I . 0 ) 


C 

890  STOP 


901  FORMAT ( I 8HREC T ANGULAR  PLATES) 

902  FORMAT (8H  NODES, 8X.I0H  THICKNESS) 
905  FORMAT ( 4  I  3 ,C i 6.8 ,  I  3) 

999  FORMAT ( IH  ) 

ENO 
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C  c 


TABLE  A- I  (CONTINUED' 


C  DORA  IV  -  BARS  AND  PLATES  (EXCEPT  TRAPEZOIDAL) 

C 

C  MATRIX  INVERSE  AND  OUTPUT  -  25  NODES  -  2  DEGREES  OF  FREEDOM 

C  *•##**#•»****#*#******##■*****#■**#*#*■**■***■*•*■*■»*■»■»»■*■»***«*»»*«***# 

c 

C  E  »L  o  COOK  -  may  1966 

C 

COMMON  NO<NOX<E<V<ST IFF(50,50).X(25)»Y<25) 

DIMENSION  P ( 50 )  <  0  <  50 ) 

C 

N02=2*N0 

PUNCH  CONSTRAINT  HEADINGS 

PUNCH  999 
PUNCH  999 
PUNCH  901 
PUNCH  999 
PUNCH  902 
PUNCH  999 

READ  THE  NUMBER  OF  A  NODE  WITH  AT  LEAST  ONE  CONSTRAIN'.  THE  X  AND 
Y  CONSTRAINT  CODES<  AND  The  CONTROL  CODE.  THE  CONSTRAINT  (ODE  IS 
C  1.0  FOR  CONSTRAINED  DISPLACEMENTS  AND  0.0  FOR  UNCONM RA I  NFC 

C  DISPLACEMENTS.  THE  CONTROL  CODES  are  the  same  A '  I  I  HAP T 

c 

130  READ  904, I ,FI <F2,KK 
PUNCH  90A  «  I . F 1 «  F2 
C 

J  =  2*  < I- )  ) 

I F ( F 1  ) 890  <  1 60  < 1  AO 
1  A  O  DO  1! O  K=I  <  N02 

STIFF  ( J+I  ,  K  )  =  0  •  0 
150  ST  IFF  (K  < J+  1  ) =0.0 

STIFF ( J+I < J+I 1=1.9 
160  IF (F? ) 890 <  220  *  170 
170  DO  I F 0  K  = I  < N02 

STIFF  ( J+2.K 1=0.0 
ISO  ST  IFF  (K  < J  +  ?  1=0.0 

STIFF ( J+2< J+2'=I .0 
220  IF (KK ) 890. 1 3O<230 
230  CONTINUE 
C 

600  DO  6C  A  I  =  I  <  NO 2 
stoR£=st:ff( i ,  i  > 

STIFF!  1,11  =  1.0 
DO  601  J= 1  » N02 

601  STIFF  <  I <J) =ST IFF(  I  .Jl/STORE 
DO  604  K= I , N02 

I  F ( K— !  1 602  <  604 , 602 

602  STORf =STIFF(K,  I  ) 

67  in- (*.«!/ =0.0 
DO  603  J= 1 , N02 

603  STIFF(K,J)=ST IFF(K.J) -STORE* ST  IFF (  I  , J) 

604  CONTINUE 
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TA8LE  A— I  (CONTINUED) 


235  DO  236  I  -  1  .N02 

236  P< I ) =0»0 
C 

PUNCH  LOAD  HEADINGS 

PUNCH  999 
PUNCH  999 
PUNCH  903 
PUNCH  999 
PUNCH  902 
PUNCH  999 
C 

C  READ  THE  NUMBER  OF  A  NODE  WITH  AN  EXTERNAL  LOAD*  THE  X  AND  Y 

C  LOADS.  AND  THE  CONTROL  CODE*  TWO  LOADS  MUST  BE  SPECIFIED* 

C  ALTHOUGH  ONE  OF  THEM  MAY  BE  ZERO.  THE  CONTROL  CODES  ARE  The  oAME 

C  AS  IN  PART  2. 

C 

238  READ  904. I ,P1 .P2.KK 
PUNCH  904* I , P 1 . P2 
C 

J*2*<  I  —  I > 

P( J+ 1 ) =P1 
P<  J+2  5  =P2 
I F ( KK  )  890  »  238  *  240 
2 40  CONTINUE 
C 

DO  250  1=1 . NO 2 
D  < I  ) *0.0 
DO  250  .J  =  l  ,  NO 2 

250  D< I »«0( I >+STIFF< I *J)«P< J) 

C. 

C  punch  DISPLACEMENT  headings 

c 

PUNCH  999 
PUNCH  999 
PUNCH  905 
PUNCH  999 
PUNCH  906 
PUNCH  999 
C 

C  DISPLACEMENTS  -  THE  NODE  NUMBER  IS  FOLLOWED  8Y  THE  X  AND  Y 

C  DISPLACEMENTS.  DISPLACEMENTS  ARE  PUNCHED  FOR  ALL  NODES.  EVEN  THE 

C  CONSTRAINED  ONES? 

C 

00  2 60  1=5. NO 
J*S* ( I -  1  ) 

260  PUNCH  904. i *0< J+l) «D( J+2) 

C 

IF (SENSE  SWITCH  2 >280. 270 

270  CONTROL  102 
PRINT  980 
PRINT  985 
CONTROL  i 02 
Pf?  I N  T  990 
PRINT  995 
280  PAUSE 
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C  TABLE  A-l  (CONTINUED) 

GO  TO  235 
890  STOP 

90 1  FORMAT (  1  1 HCONSTRA I  NTS ) 

902  FORMAT ( 5H  NODE • 6X ♦ 2H  X.14X.2H  V) 

903  FORMAT ( 5HL0A0S > 

904  FORMAT (  I  3. 2£ 16.8. 1 3) 

905  FORMAT ( 1 3HD I SPLACEMENTS ) 

906  FORMAT ( 5H  NODE . 6X « 2H  U.14X.2H  V) 

980  FORMAT ( 49H  FOR  A  NEW  PROBLEM  OR  FOR  A  CHANGE  OF  CONSTR A I  NTS  *  1 H .  ) 
985  FORMAT (45H  THE  ENTIRE  SET  OF  PROGRAMS  MUST  BE  RELOADED® > 

990  FORMAT ( 47H  HOWEVER  <  ANY  NUMBER  OF  ADDITIONAL  LOAD  VECTORS) 

995  FORMAT ( 45H  MAY  BE  LOAOED  at  THIS  TIME.  DEPRESS  *STAPT*> 

999  FORMAT ( 1 H  ) 

END 


c  c 


TABLE  A-2 


DIRECT  STIFFNESS  ANALYSIS  -  PANEL 
DORA  IV  -  25  NODES  -  2  DEGREES  OF  FREEDOM 


NO  OF  YOUNG*S 
NODES  MODULUS 


P01SS0N*S 

RATIO 


12  • 1 OOOOOOOE  06  .31600000E  00 


NODE  X 


Y 


1  .OOOOOOOOE-50 

2  .40000000E  01 

3  .6000000GE  01 

4  .OOOOOOOOE-SO 

5  « 4 OOOOOOOE  0  1 

6  .60000000E  01 

7  #0 OOOOOOOE "50 

8  .40000000E  01 

9  .60000000E  01 

10  .OOOOOOOOE-50 

11  « 40000000E  0  1 

12  • 60000000E  01 


.OOOOOOOOE-50 

•OOOOOOOOE-50 

.OOOOOOOOE-50 

•  80000000  E  01 
.80000000E  01 
•80000000E  01 

•  1 6000000E  02 
« 1 6000000 E  02 
. 1 6000000E  02 
•24000000E  02 
• 24000000E  02 
• 24000000E  02 


BARS 

NODES  AREA 

1  2  • 30 1 OOOOOE  00 

2  3  . 30 1 OOOOOE  00 

4  5  «  30 1 OOOOOE  00 

5  6  . 30 1 OOOOOE  OC 

7  8  . 30 l OOOOOE  00 

8  9  • 30 1 OOOOOE  Oo 

10  11  . 30 1 OOOOOE  Oo 

11  12  .301 OOOOOE  00 

1  4  .60200000E  00 

2  5  .60200000E  00 

4  7  .60200000E  Oo 

5  8  .60200000E  00 

10  .60200000E  00 

8  1 1  .60200000E  00 
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c  c 


TABLE  A— 2  (CONTINUED) 


RECTANGULAR  PLATES 

NODES  THICKNESS 

1245  • 1 OOOOOOOE  00 

2356  • 1 OOOOOOOE  00 

4578  « 1 OOOOOOOE  00 

5  to  8  9  .lOOOCoOOE  00 

7  81011  . 1 OOOOOOOE  00 

8  91112  . 1 OOOOOOOE  00 


CONSTRAINTS 

NODE  X  Y 

3  • 1 OOOOOOOE  01  .00000000E-50 

6  •  1  OOOOOOOE  01  «00000000E-50 

9  • 1 OOOOOOOE  01  .00000000E-50 

12  « 1 OCOOOOOE  01  . 00000000E-50 

2  • 0 OOOOOOOE -80  .10000000E  01 


LOADS 

NODE  X  Y 

10  .00000000E-50  • 50000000E  00 

11  • 000000 OOE -50  • 50000000E  00 


DI SPt  ACEMENTS 

MODE  U  V 


1  - ♦ 505526 47E-0 7 

2  -.28441 223E-07 

3  .OOOOOOOOE-50 

4  . 1 4233045E-06 

5  .61 20939 3E -07 

6  .00000000E-50 

7  .8031 6400E-07 

8  .2640631 9E-07 

9  .00000000E-50 

10  . 798 3 38 70 E-0 7 

11  .29270008E-07 

12  .00000000E-50 


• 752524 76E- 06 
.00000000E-50 
. 15515326E-06 
•89208273E-06 
.6891 2770E-06 
.6352647  IE-06 
. 1 2602640E-05 
. 1 1859289E-05 
. 1 2038422E-05 
. 1 7228336E-05 
. 1619891 1 E- 05 
. 1 522209 1 E-05 
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TABLE  A-2  (CONTINUED) 


NODES 


STIFFENER  STRAINS 

NODES  STRAINS 

!  2  .55278560E-08 

2  3  . 1422061  IE-07 

4  5  - • 20280265E-C7 

5  6  - • 3 060 469 8E -07 

7  8  1 3477520E-07 

8  9  13203159E-07 

10  II  -» 12640965E-07 

11  12  -« 1 4635004E-07 

1  4  . 1 7444747E-07 

2  5  .86 1 40962E-07 

4  7  . 460226628-07 

5  8  .62100 550E-07 

7  10  .57821200E-07 

8  1  1  .54245275E-07 


WEB  STRAINS 

NODES 

12  4  5 

2  3  5  6 

4  5  7  8 

5  6  8  9 

7  81011 

8  91112 


EXX 

-.73762040E-06 
- .8 1920420E-08 
- . 1 6878892E-0  7 
- .21903928E-07 
1 305924  3E-0  7 
- . 1 39 1 9082E-0  7 


EVY 

•51 792854E-07 
. 73077446E-07 
.5406 1 404E-07 
.66586 169E-07 
. 560332 30E-07 
. 47020569E-07 


EXY 

-. 101 77662E-0b 
.3C925730E-07 
-.4071 2350E-07 
-.11 162630E-07 
-.22010880E-07 
-. 19763200E-07 


c  c 


TABLE  A— 3 


DIRECT  STIFFNESS  ANALYSIS  ~  PANEL  CONFIGURATION  II 
DORA  IV  -  25  NODES  -  2  DEGREES  OF  FREEDOM 


NO  OF  YOUNG*S  POISSON*S 

NODES  MODULUS  RATIO 

24  • 1 OOOOOOOE  08  .31600000E  00 


NODE  X  Y 


1 

•OOOOOOOOE 

-50 

•OOOOOOOOE- 

-50 

2 

• 20000000E 

0  1 

. OOOOOOOOE- 

-50 

3 

.40000000E 

0  1 

•OOOOOOOOE- 

-50 

4 

.60000000E 

0  1 

. OOOOOOOOE- 

-50 

5 

.OOOOOOOOE 

-50 

. 3000000CE 

0  1 

6 

• 20000000E 

0  1 

. 30000000E 

0  1 

7 

.40000000E 

0  1 

• 30000000E 

01 

8 

• 60000000E 

0  1 

• 30000000E 

01 

9 

•OOOOOOOOE 

-50 

. 50000000E 

0  1 

1  0 

•20000000E 

0  1 

• 50000000E 

01 

1  1 

•  40000000E 

0  1 

• 50000000E 

01 

12 

•60000000E 

0  1 

• 50000000E 

01 

13 

•OOOOOOOOE 

-50 

• 80000000E 

01 

1  4 

•20000000E 

01 

•80000000E 

01 

15 

•40000000E 

0  1 

• 80000000E 

01 

16 

•60000000E 

0  1 

• 80000000E 

01 

17 

•OOOOOOOOE- 

-50 

- • 80000000E 

0  1 

18 

•  20000000E 

0  1 

- • 80000000E 

01 

19 

•  40000000E 

0  1 

- •8000000CE 

01 

20 

•60000000E 

C  1 

— • 80000000E 

01 

2! 

•00000000E-50 

-. 1 6000000b 

02 

22 

• 20000000E 

0  1 

1 6000000E 

02 

23 

•40000000E 

0  I 

-• 1 6000000E 

02 

24 

•60000000E 

01 

-. 1 6000000E 

02 
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c  C  TABLE  A- 

8APS 


NODES  AREA 


1 

2 

.30 1 OOOOOE 

00 

2 

3 

.30 1000G0E 

00 

3 

4 

.30 1 OOOOOE 

00 

13 

1  4 

.301 OOOOOE 

00 

1  4 

15 

•  30 1 OOOOOE 

00 

1  5 

16 

•  30  1  OOOOOE 

00 

1 

5 

. 60200000E 

00 

3 

7 

.60200000E 

00 

5 

9 

• 60200000E 

00 

7 

1  1 

• 60200000E 

00 

9 

1  3 

• 60200000E 

00 

1  1 

1  S 

. 60200000E 

00 

17 

ie 

.301 OOOOOE 

00 

18 

19 

.301 OOOOOE 

00 

19 

20 

. 30 1 OOOOOE 

00 

21 

22 

. 30 1 OOOOOE 

00 

22 

23 

•  30 1 OOOOOE 

00 

23 

24 

.30  1  OOOOOE 

00 

1 

17 

• 60200000E 

O'-' 

3 

19 

•60200000E 

00 

17 

21 

•60200000E 

00 

19 

23 

• 60200000E 

00 

rectangular  plates 

nodes  THICKNESS 


1 

2 

5 

6 

.  1  OOOOOOOE 

00 

2 

3 

6 

7 

.  1  OOOOOOOE 

00 

3 

4 

7 

8 

.  10000000E 

00 

5 

6 

9 

1  0 

. 1 OOOOOOOE 

00 

6 

7 

1  0 

1  1 

.  1  OOOOOOOE 

00 

7 

8 

l  1 

12 

.  1  OOOOOOOE 

00 

9 

1  0 

1  3 

1  4 

.  1  OOOOOOOE 

00 

1  0 

1  1 

14 

1  5 

.  1 OOOOOOOE 

00 

1  1 

12 

1,5 

16 

. 1 OOOOOOOE 

00 

17 

18 

1 

2 

. 1 OOOOOOOE 

00 

18 

19 

2 

3 

.  1  OOOOOOOE 

00 

19 

20 

3 

4 

. 1 OOOOOOOE 

00 

21 

22 

17 

18 

.  1 OOOOOOOE 

00 

22 

23 

13 

1  9 

.  i OOOOOOOE 

00 

23 

24 

19 

20 

.  1  OOOOOOOE 

00 

constraints 


NOOE 

X 

Y 

24 

» 1 OOOOOOOE 

0  1 

•00000000E-50 

20 

.  1  OOOOOOOE 

C  1 

.00000000E-50 

4 

.  1  OOOOOOOE 

01 

•00000000E-50 

8 

.  1  OOOOOOOE 

0  1 

•OOOOOOOOE— 50 

12 

.  1 OOOOOOOE 

0  1 

.OOOOOOOOE -50 

16 

•  1 OOOOOOOE 

c  1 

.00000000C-50 

23 

.OOOOOOOOE- 

■50 

. 1 OOOOOOOE  0 1 

(CONTINUED) 
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TABLE  A-3  ( CONT  I  NUED  > 


LOADS 

NODE  X  Y 

13  • OOOOO OOO E~ 50  *50000000E  00 
15  .0C1000000E-50  «50000000E  00 


displacements 


node  u  V 


1 

.80047330E-07 

.  127141 39E-05 

2 

•  5658 1 539E-0  7 

. 1 2398273E-05 

3 

.30157771E-07 

. 1 1 966348E-05 

4 

•00000000E-50 

. 1 2037975E-05 

5 

•  8457  5  695E-07 

. 1 4385289E-05 

6 

•51000175E-07 

. 1 3982069E-05 

7 

•23830666E-07 

. 1 3632036E-05 

8 

•OOOOOOOOE-50 

•  1 3570573E-05 

9 

• 6 1 306633E— 07 

•  ! 556960  IE-05 

1  0 

•  4 1 S34539E— 07 

. 1 4966693E-05 

1  1 

•2 1 93697 2E— 07 

.  1  4705235E-05 

12 

•00000000E-50 

•  1 4501 31 6E-05 

1  3 

•75304050E-07 

. 1 7583445E-05 

1  4 

•55530950E-07 

. 1 544621 4E-05 

15 

•2932791 6E-07 

. 16494427E-05 

16 

•00000000E-50 

. 1 4 99 37 56 E -05 

1  7 

•  1 3537860E-06 

•90784096E-06 

18 

. 1 1062759E-06 

•78620528E-06 

19 

•  601 6607  OE  — 07 

•  70 1 26854E-06 

20 

•  000 00 00 0E~ 50 

•64637866E-06 

21 

-.48537227E-07 

•  7 540803 7E- 06 

22 

— • 4 6085 29 4 E— 07 

•46944948E-06 

23 

-•27554959E-07 

•00000000E-50 

24 

•00000000E-50 

. 1 58571 47E-06 
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TABLE  A-3  (CONTINUED) 


STIFFENER  STRAINS 


NODES  STRAINS 


1 

2 

-.11 732895E-07 

2 

3 

-.1321 1 884E-07 

3 

4 

-. 1 5078885E-07 

13 

1  4 

- . 9886 5 50 OE -08 

1  4 

1  5 

-. 1310151 7E-07 

15 

16 

-. 1 4663958E-07 

1 

5 

.55705000E-07 

3 

7 

. 55 522 9 3 3E -07 

5 

9 

•59215600E-07 

7 

1  1 

. 53659950E-07 

9 

1  3 

•6712813 3E -07 

1  ! 

15 

. 59639733E-07 

17 

18 

-.  12375505E-07 

18 

19 

- • 25230  760 E -07 

1  9 

20 

- • 30083 035E -07 

21 

22 

. 1 2259665E-0Q 

22 

23 

• 9265 1 675E-08 

23 

24 

•  1 37774 79E -07 

1 

17 

.454466 1 2E-07 

3 

19 

.6! 920787E-0? 

17 

2  1 

. 1 9220073E-07 

19 

23 

•87658567E-07 

WEB 

STRAINS 

NODES  EXX 

E  VY 

EXY 

1 

2 

5  6  -. 14259328E-07 

•54 24 9090 E-07 

-. 18153320E-07 

2 

3 

6  7  -.  13398320E-07 

.54 1580S0E-07 

-.21533710E-07 

3 

4 

7  8  -.  13497109E-07 

•  53304 760E-07 

— .80041000E— 09 

5 

6 

9  10  -. 1 3260903E-07 

. 54 2234 00 E-07 

-•33260870E— 07 

6 

7 

10  11  -.1  1  766770E  — 0  7 

•  5 1 44  5 580 E— 0  7 

-. 18052120E-07 

7 

8 

11  12  -.  1  1441910E-07 

•50 098 5 50 E-07 

-.71 079800E-08 

9 

1  0 

13  14  -.961 12980E-08 

.4 1 556080E-07 

-.63887840E-07 

1  0 

1  1 

14  15  -.  1  1525151E-07 

.3761 1 S90E-07 

.23 1 8344  OE  —0  7 

i  1 

12 

15  16  -.  12816222E-07 

.  38027200E-07 

-.4 1382930E— 07 

1  7 

18 

1  2  -.  12054200E-07 

.51 074686E-07 

-.45141 640E-07 

18 

19 

2  3  -.  19221 322E-07 

.5931 1767E-07 

-.3721  3700E-07 

19 

20 

3  4  -.2258096 1 E-07 

. 65799069E— 07 

-. 1 380731 OE-07 

21 

22 

17  18  - .55747690E-C8 

•29407274E-07 

-.80277350E-07 

22 

23 

18  19  - . 7982 7960 E -08 

.6362652 1 E-07 

-. 12331944E-06 

23 

24 

19  20  -.S1527780E-08 

.7431 7233E-07 

•  3 1 4 0296 0E-0  7 

c  c 


TABLE  A-4 


DIRECT  STIFFNESS  ANALYSIS  -  PANEL  CONFIGURATION  III 
DORA  IV  -  25  NODES  -  2  DEGREES  OF  FREEDOM 


NO  OF  YOUNG»S  POISSON*S 

NODES  MODULUS  RATIO 


24 

•  1 OOOOOOOE 

08 

. 316000OOE 

00 

NODE 

X 

Y 

1 

•  OOOOOOOOE 

-50 

•OOOOOOOOE- 

-50 

2 

•  20000000E 

0  1 

.OOOOOOOOE- 

-50 

3 

•40000000E 

01 

.OOOOOOOOE- 

-50 

4 

•60000000E 

0  1 

•  OOOOOOOOE- 

-50 

.OOOOOOOOE 

-50 

• 30000000E 

01 

6 

•20000000E 

0  1 

. 30000000E 

01 

7 

.40000000E 

0  1 

. 30000000E 

Ol 

8 

.60000000E 

Cl 

. 30000000E 

Ol 

9 

.OOOOOOOOE- 

-50 

. 50000000E 

01 

1  0 

•  20000000E 

01 

• 50000000E 

Ol 

1  1 

•  40000000E 

01 

• 50000000E 

Ol 

12 

.  60000000E 

01 

. 50000000E 

0  1 

1  3 

.OOOOOOOOE- 

-50 

•80000000E 

01 

1  4 

•20000000E 

01 

• 80000000E 

Ol 

15 

.  40000000E 

0  1 

. 80000000E 

01 

16 

•60000000E 

01 

• 80000000E 

Ol 

1  7 

.OOOOOOOOE- 

-50 

. 1 6000000E 

02 

1  8 

•20000000E 

01 

. 1 6000000E 

02 

19 

.  4  OOOOOOOE 

01 

. 1 6000000E 

02 

20 

.60000000E 

01 

. 16000000E 

02 

21 

.OOOOOOOOE- 

-50 

-.80000000E 

Ol 

22 

•  20000000E 

01 

-.80000000E 

01 

23 

•  40000000E 

01 

- • 80000000E 

01 

24 

.60000000E 

01 

- • 80000000E 

Ol 

A-2  7 


' '  IT 
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C  C  TABLE  A-4 

BARS 


NODES  AREA 


1 

2 

•  30 1 OOOOOE 

00 

2 

3 

r 30 1 OOOOOE 

00 

3 

4 

•  30 1 OOOOOE 

00 

1  3 

1  4 

»  30 1 OOOOOE 

00 

1  4 

15 

.301 OOOOOE 

00 

1  5 

16 

.301 OOOOOE 

00 

1 

5 

•60200000E 

00 

3 

7 

.60200000E 

00 

5 

9 

. 60200000E 

00 

7 

1  1 

.6C200000E 

00 

9 

1  3 

.60200000E 

00 

i ; 

15 

.  30200009E 

00 

17 

18 

•  30 1 OOOOOE 

00 

18 

19 

.30 100000P 

00 

19 

20 

•  3C  1 OOOOOE 

00 

21 

22 

.30 1 OOOOOE 

00 

22 

23 

.301 OOOOOE 

CO 

23 

24 

.30 1 OOOOOE 

00 

1 

21 

.60200000E 

00 

3 

23 

.60200000E 

00 

1  3 

1  7 

.60200000E 

00 

1  5 

19 

.60200000E 

00 

rectangular  plates 


NODES  THICKNESS 


1 

2 

5 

6 

1 OOOOOOOE 

00 

2 

3 

6 

7 

1 OOOOOOOE 

00 

3 

4 

7 

8 

1 OOOOOOOE 

00 

5 

6 

9 

1  0 

1 OOOOOOOE 

00 

6 

7 

1  0 

1  1 

1 OOOOOOOE 

CO 

"7 

e 

1  1 

1  2 

1 OOOOOOOE 

00 

9 

1  0 

i  3 

1  4 

1 OOOOOOOE 

00 

1  0 

i  i 

1  4 

15 

ioooooooe 

00 

1  1 

12 

15 

1  6 

1 OOOOOOOE 

00 

21 

22 

1 

2 

o 

1 OOOOOOOE 

00 

22 

23 

2 

3 

1 OOOOOOOE 

00 

23 

24 

3 

« 

leoooeocE 

00 

1  3 

1  4 

I  7 

1  8 

1 OOOOOOOE 

00 

I  A 

1 

«  O 

a  *  j 

1  / 

I  OOCwOCCc. 

00 

1  5 

t  £ 

1  9 

kC 

i  OOOOOOOE 

00 

constraints 


node  X 


r 


24 

» 1 OOOOOOOE 

0  1 

.OOOOOOOOE-50 

8 

. 1 0000 ^00 K 

0  i 

vOOooOOOOE-SO 

4 

.  1  OOOOOOOE 

0  1 

.OOOCOOOOE -50 

1? 

. 1 OOOOOOOE 

o ; 

.0000C000E-50 

16 

.  1 OOOOOOOE 

0  1 

.OOOOOOOw. >50 

20 

.  i OOOOOOOE 

0  1 

•OOOOOOCOE-50 

23 

.OOOCOOOOE- 

■50 

. 1 OOOOOOOE  0 1 

(CONTINUED) 


A-28 


c  c 

LOADS 

NODS  X 


TABLE  A- 


Y 


17  .00000000E-50  .50000000E  00 
19  .00000000E-50  .50000000E  00 


D 1 SPLACEMENTS 


NODE  U  V 


1 

.  1 3206989E-06 

.9157781 GE-06 

2 

.  1  1  225373E-06 

. 7832 1 554 E— 06 

3 

.61 489530E-07 

.7020761 1 £—06 

A 

.00000000E-30 

6364  5268E-06 

5 

.  1  43200S9E-06 

.  1 033622 1 E-05 

6 

. 1 0357874E-06 

.9656 1 388E-06 

7 

•  5702004 7E  — 0  7 

.89863685E-06 

9 

.  00000  0 00 E —50 

.8878805 1 E— 06 

9 

.  1  25S8823E-06 

.1 1 2471 2 1 E-05 

1  0 

•87080298E— 07 

.  1 0823834E-05 

1  1 

.445021 05E-07 

.  1  024 1 376E-05 

12 

.ooooooooe-so 

.  1 0237987E-05 

1  3 

. 6258 798 6E— 07 

. 1 2747602E— 05 

1  4 

.43476770E-07 

. 1 2 60 662 9 E— 05 

1  5 

.  2256655 3E~ 07 

. 1 1967597E-05 

16 

.00000000E-50 

.121 30666E-05 

1  7 

.8381 0380E-07 

.  1 7546426E-05 

1  6 

.  59394 033E— 0  7 

. 1 5908938E-05 

1  9 

•  30876461  t  — 0  7 

. 1 643756 1 E— 05 

20 

.00000000E-50 

. 1 5422935E-05 

21 

-.50281  472E-07 

•75928189E-06 

22 

-.485518 1 6E-07 

.475 1 4069E-06 

23 

-.29089571 E-07 

eOOOOOOOOE— 50 

24 

.00000000E-50 

. 16121 557E-06 

(CONTINUED) 


A-29 


c  c 


TABLE  A-a  (CONTINUED) 


STIFFENER  STRAINS 
nodes  STRAINS 


1 

2 

-.99080800E-08 

2 

3 

-.253821005-07 

3 

4 

-.30744765E-07 

13 

1  4 

-  .  95556 OBOE -08 

1  4 

15 

-.1 04551 08E-07 

15 

16 

1  1 283276E-07 

1 

5 

.39281 333E-07 

3 

7 

• 65520 24 7E -07 

5 

9 

•45545000E-07 

7 

1  1 

•62750400E-07 

9 

1  3 

.5001 6033E-07 

1  1 

1  £ 

.57540700E-07 

1  7 

18 

-.  122081 73E-07 

1  8 

19 

-.  1  4258786E— 07 

19 

20 

-.  1  54  38230E-0  7 

2  1 

22 

.86482800E-09 

22 

23 

.9731 1 225E-06 

23 

24 

. 1 4  544  78 5E  —  0  7 

1 

21 

. 19562027E-07 

3 

23 

.877595 1 4E-07 

13 

17 

. 59985300E-07 

15 

19 

•55874550E-07 

WEB 

strains 

nodes  exx 

EYV 

EXY 

1 

2 

5  6  -. 1 4859578E-07 

.5O040380E-07 

— . 4973337 OE— 07 

2 

3 

6  7  - .24330723E-07 

.63159840E-07 

-.3921 9870E-07 

3 

4 

7  8  - .29627395E-07 

.74664760E-07 

-. 1 9S39860E-07 

5 

6 

9  10  -.  19532521 E-07 

.519648B0E-07 

-.361 1 2020E-07 

6 

7 

10  11  -.22284222E-07 

•60567570E-07 

-.3855981 0E-07 

7 

8 

11  12  - .25380S78E-07 

-65354740E-07 

- »  590 3280 OE- 08 

Q 

10 

13  14  — . 1 440478 7 E  —  07 

.54721 250E-07 

-.31 873790E-07 

1  0 

1  1 

14  15  -. 158721 04E-07 

.58483590E-07 

-.41 460430E— 07 

1  1 

1  2 

15  16  -. 16767164E-07 

.60315000E-07 

.  33607000E-r'9 

21 

22 

1  2  -.45216260E-08 

•29035691 E-07 

-.32728630E-07 

22 

23 

2  3  - .79254890E-08 

.63 1 34435E-07 

-. 1 2 33584 9E -06 

23 

24 

3  4  - .80999900E-C8 

.73582077E-07 

•  295592 30 E-0  7 

1  3 

1  4 

17  18  - . 1 088 1 flo 1 F-07 

,  sOMPomr-07 

-.421 407902-07 

1  4 

15 

18  19  -  .  1  2356948E-0"? 

-  455767  1  Ui;  -07 

- . , 2460200E -08 

15 

16 

19  20  - . 1 336 2754E-0V 

•  4  f)5  1  3964 C  -07 

-  .  2076956OE  -O'7 

A  -  30 
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c  C  TABLE  A-5 

DIRECT  STIFFNESS  ANALYSIS  -  PANEL  CONFIGURATION  IV 
DORA  IV  -  25  NODES  -  2  DEGREES  OF  FREEDOM 


no  of 

YOUNG*S 

PO I SSON*S 

NODES 

MODULUS 

RATIO 

24 

• i oooooooe  oa 

•  31 600000E  00 

node 

X 

Y 

1 

.00000000E-50 

• OOOOOOOOE 

-50 

2 

• 20000000E  0 1 

.OOOOGOOOE- 

-50 

3 

•  40000000E  01 

.OOOOOOOOE- 

-50 

4 

.60000000E  01 

.OOOOOOOOE 

-50 

5 

•00000000E-50 

. 30000000E 

01 

6 

• 20000000E  01 

. 30000000E 

0  1 

7 

.400C0000E  01 

. 30000000E 

0  1 

a 

•  60000000E  01 

. 3G000000E 

01 

9 

.00000000E-50 

.50000000E 

01 

1  0 

• 20000000E  01 

. 50000000E 

0  1 

1 1 

.400C0000E  01 

.50000000E 

0  1 

12 

.60000000E  01 

.50000000F 

0  1 

1  3 

•000C0000E-50 

.eoooooooF 

0  1 

14 

• 20000000E  01 

.aooouoooE 

0  1 

15 

.40000000E  01 

.aoooooooE 

01 

16 

. 60000000E  Cl 

»  80000000E 

0  1 

1  7 

.00000000E-50 

.  16000000E 

02 

18 

• 20000000E  01 

. 1 6000000E 

02 

1  9 

•  4 0000 OOOF  01 

. 1 6000000E 

02 

20 

• 60000000E  01 

.  1  6000000E 

02 

21 

.00000000E-50 

.24000000E 

02 

22 

• 20000000E  01 

.24000000E 

02 

23 

•  4 OOOOOOOE  01 

. 24000000E 

02 

24 

.60000C00E  01 

.24000000E 

02 

A-31 
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c  C  TA8LE  A-5  (CONTINUED) 

BARS 


NODES  AREA 


1 

2 

.30 1 OOOOOE 

00 

2 

3 

.30 100000E 

00 

3 

4 

.30  I OOOOOE 

00 

1  3 

14 

.30 1 OOOOOE 

00 

1  A 

1  5 

«  30 1 OOOOOE 

CC 

15 

16 

.30 ! OOOOoE 

00 

1 

5 

•6020OO00E 

00 

3 

7 

•602GO000E 

00 

5 

9 

.6020000CE 

00 

7 

1  1 

•60200000E 

oo 

9 

I  3 

•60200O00E 

00 

1  1 

1  5 

.602G0000E 

00 

I  7 

18 

.30  I OOOOOE 

00 

1  8 

I  9 

.30  I OOOOOE 

00 

1  9 

20 

.30 I OOOOoE 

oc 

21 

22 

.30  I OOOOOE 

00 

22 

23 

.30! OOOOOE 

00 

23 

24 

.30 I OOOOOE 

00 

1  7 

21 

•60200000E 

oo 

1  9 

23 

•60200000E 

00 

1  3 

1  7 

. 60200000E 

00 

1  5 

19 

.60200000E 

Oo 

RECTANGULAR  plates 

NODES  THICKNESS 


1 

2 

5 

6 

.  1  OOOOOOOE 

00 

2 

3 

6 

7 

„  1  OOOOOOOE 

00 

3 

4 

7 

3 

.  1 OOOOOOOE 

00 

? 

6 

9 

1  0 

.  i  ooooooce 

00 

6 

7 

10 

1  1 

. 1 OOOOOOOE 

00 

7 

8 

1  1 

1  2 

,  1 OOOOCOUE 

00 

9 

!  0 

1  3 

1  4 

.  1 OOOOOOOE 

CO 

1  0 

1  1 

1  4 

15 

.  10C00000E 

00 

1  1 

12 

1  5 

16 

.  1  OOOOOOOE 

00 

1  3 

1  4 

1  7 

1  8 

.  1  OOOOOOOE 

00 

1  4 

55 

!  3 

1  9 

.  1  OOOOOOOE 

00 

1  5 

1  6 

1  9 

20 

.  1  OOOOOOOE 

00 

1  7 

18 

21 

22 

«  !  OOOOOOOE 

00 

l  8 

!  9 

2  ? 

28 

•  1  OUOUOOOt 

oc 

1  9 

20 

23 

2« 

.  1  OOOOOOOE 

00 

CONSTRA 1  NTS 

NODE  X  Y 

20  » 1 OOOOOOOE  0!  . OOOOOOOOE-50 

8  . ICOOOOGCE  0!  .OOOOOOOOC-5C 

A  .1 OOOOOOOE  01  .00D00OC0E-5O 

12  .1 OOOOOOOE  01  .00000000E-50 

16  .-OOOOOOOE  0!  .OOOOOCOOE -SO 

24  .1 OOOOOOOE  01  .OOOOOOOOE-SO 

2  .OOOOOOOOE-SO  . 1 OOOOOCGE  01 


A-32 


c  c 

LOADS 

NODE  X 


TABLE  A- 


Y 


21  ,0000C000E-50  .50G00000E  00 
23  .OOOOOOOOE-SO  . 50000000E  00 


D I splacements 


NODE  U  V 


1 

10867604E-06 

•84300G32E-06 

2 

- *986322558 -07 

.55634454E-0S 

3 

-.55491 384E-07 

.ooooooooe-50 

a 

» 00000 000E-50 

. 2532477 1 E-06 

5 

. 1 1 0979506 -06 

.85770273E-06 

6 

. 1 081 923 1 E-06 

. 5966 4 25 3E— 06 

7 

c  52 l 4  7767E-0  7 

.33207062E-06 

s 

.00000000E-50 

. 3 1 62 094 5E— 06 

9 

. 1 5500238E-06 

. 88665 4 4 6E— 06 

1  0 

.12517964E-06 

.7000053CE— 06 

l  1 

.  70666772E-0"? 

.50966276E-06 

12 

•00000000E-50 

,490448096-06 

13 

.  10986507E-06 

.970017066-06 

14 

•  8564  7 1 65E-0  7 

.8581 67S6E-06 

15 

.4671 209 1 E~0  7 

.74349165E-C6 

16 

.OOOOOOOOE-50 

. 730 25990 E-C 6 

1  7 

.7640  31 65E-0  7 

. 1 32700 1 9E-05 

18 

.53168969E-07 

. 1 3C48632E-05 

19 

.27581 046E-07 

.  12456  1526-05 

20 

.ooooooooa-50 

. 1 2568298E-05 

21 

.831 39 1 30E— 07 

. 1803661 3E-05 

22 

. 59234  0to5E-0  7 

. 1 o438256E~05 

23 

.30847897E-07 

. 1 693902 1 E-05 

24 

.00000000E-50 

. 1 5950  362 E-05 

(CONTINUED) 
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C  C  TABLE  A-5  (CONTINUED) 

STIFFENER  STRAINS 


nodes  STRAINS 


I 

2 

.50218950E-08 

2 

3 

.21 570435E-07 

3 

4 

.27745693E-07 

13 

14 

121 08955E-07 

1  4 

IS 

1 9467537E-07 

15 

16 

-.2335604 5E- 07 

1 

5 

. 4900803 3E -08 

3 

7 

•  1  1 069020E -06 

5 

9 

.  1 4  475865E-0  7 

7 

1  1 

•88796070E-07 

9 

13 

• 2778753 3E-07 

1  1 

15 

.77942963E-07 

1  7 

18 

-. 1 1 61 7098E-07 

18 

19 

-. 1 279396 1 E-07 

19 

20 

-. 1 3790523E-07 

21 

22 

-. 1 1 952532E-07 

22 

23 

-. 14 193084E— 07 

23 

24 

-.1 542 394 8e -07 

17 

21 

. 59607425E— 07 

19 

23 

•56035862E-07 

13 

17 

.44623 100E-07 

1  3 

19 

.62765450E-07 

WEB 

strains 

nodes  exx 

EYY 

EXY 

1 

2 

5  6  . 18141490E-08 

.91 667280E-08 

- . 65848980E -0  7 

2 

3 

6  7  - • 32259 1 30E-08 

.62061 433E-07 

-. 15281851E-06 

3 

4 

7  8  • 83590 400E -09 

.65B38724E-07 

.77286494F-07 

5 

6 

9  10  -.81 5248 20 E -08 

. 33078640E-07 

- . 96674800E-07 

6 

7 

10  11  — .2763935 3E -07 

.70238730E-07 

-.  1  0485204E-06 

7 

8 

11  12  -.30703635E-07 

. 8 795 769 0E~ 07 

-.41392100E-08 

9 

10 

13  14  13510162E-07 

.40254 1 80E-07 

-.88736220E-07 

1  0 

1  1 

14  15  -.23361 985E-07 

.65331 900E-07 

- • 86835S90E-07 

1  1 

12 

15  16  -.2934471 6E-07 

.789401 10E-07 

-.12 1 04050E-07 

13 

1  4 

1718  - . 1 1 S63026E-07 

•502300 10E-07 

-.3761 8240E-C7 

1  4 

15 

18  19  -.  16130749E-07 

.59301 180E-07 

- . 4 670665 OE -07 

15 

16 

19  20  -.  18573285E-07 

.6429334  IE-07 

-.  1 6999 70 OE  —08 

1  7 

18 

21  22  -. 1 1784816E-07 

.50988860E-07 

— .44743540E— 07 

18 

19 

22  23  - . 1 3493523E-07 

.49203080E-07 

- • 1 (Uy6JUUt-UB 

19 

20 

23  24  -.  14607236E-07 

.491  55333E--07 

-.21 70 8660 E-07 
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c 


c 


TABLE  A-6 


STIFFENER  AND  WEB  STRAIN  PROGRAM  LISTING 

C  STIFFENER  and  WEB  STRAINS  FROM  NODAL  DISPLACEMENTS 

c  *********************************************************** 

c 

C  E*L.  cook  -  JUNE  1966 

C 

DIMENSION  X( 25 > • Yt 25) ,U<25)  ,V<25> *  I  A (25) , !B(25) . A < 25 J 
DIMENSION  JA(  16) » JB(  16) »  JC<  16) • JD<  16) «  f<  16) 

DIMENSION  DELTA(O) »SIGMA< S) 

C 

C  READ  NUMBER  OF  NODES,  YOUNGS  f  JOULUS,  AND  POISSONS  RATIO 

C 

801  READ  901, NO « E  «  VV 
EP-E/ < 1 ,0-VV**2 ) 

GS=E/<2.0*< 1.0+VV) ) 

READ  NODAL  COORDINATES 

DO  10  I = I . NO 

802  READ  901 ,K,X<  I > . Y<  I  ) 

I F ( K  —  I  1700. 10,700 

10  CONTINUE 

read  STIFFENFR  data 
IK=  1 

803  REAO  902, IA< IK) , IB( IK) ,A< IK) , ICODE 
IK=IK+I 

I F (  I  CODE- 1  >803,30,700 
30  IK=lk-l 

READ  WEB  DATA 

JK=  1 

804  REAO  90  3, JA( JK ) , JB< JK )  , JC ( JK ) , JD ( JK ) , T ( JK)  ,  I  CODE 
JK  =  JK-*  1 

I  r  <  I  CODE- 1  >804,50,700 
50  JK  =  JK -  1 

read  nodal  displacements 

DO  60  1=1, NO 

fine  REaD  90! : JtU! ! ! sV! ! ! 

IF (  I -J ) 700,60 , 700 
60  CONTINUE 


A-  3  5 


n  n 


* 


c  C  TABLE  A-6  (CONTINUED) 

C  CALCULATE  STIFFENER  STRAINS 

C 

DO  607  I  L  =■  1  •  1  K 
I  A  I  =  I  A  (  I  L  ) 

X  A  =  X (  I  A  I  ) 

V  A  =  Y (  I  A  I  ) 

UA  =  U<  I  A  I  ) 

V  A  =  V  (  I  A  I  ) 

I  B I  =  I B  (  ID 
XB“X (  IB  I  ) 

YB-Y (  I B I  ) 

U8=U ( IBI) 

VB  =  V(  IBI  ) 

ALX  =  XB  —  X A 
AL Y  =  YB—  Y  A 
AL2=ALX»*2+ALY**2 

ESTR= ( ALX# ( UB-UA ) +  ALY*  <  VB- VA )  >/AL2 
BO  7  PUNCH  902  .  I  A  J  ,  IBI  » ESTR 

CALCULATE  WEB  STRAINS 

DO  806  I L  = 1  . JK 
J  A  I  =  J  A (  IL) 

X  A  =  X ( JA I  ) 

YA  =  Y( JA I  ) 

DELTA (  1  )=U< JAI  ) 

DELTA(B)=V(  JAI  ) 

J8 I = J8  <  IL) 

XB  =  X ( JB I  ) 

DELTA ( 2 ) =U ( J5  I  ) 

DELT  A ( 6 ) =V ( JB1  ) 

JC I =  JC  <  I L ) 

XC  =  X<  JC I  ) 

YC  =  Y ( JC I  ) 

DELT  A ( 3) =U<  JC I  > 

DELTA ( 7 ) = V ( JC I  ) 

JO  I = JO (  I L  ) 

XD=X< JD1 ) 

DELTA ( 4 ) =U<  JD!  ) 

DELTA( 8  >  =V  <  JO  I  ) 

XBA  =  XB  —  X  A 
P= < XC-XA ) /XBA 
D  = ( XD-XA ) /XBA 
B  =  D-  1 .0 
TR-:  1  .0/  (  n-P  ) 

YCA=YC-YA 
E  T  A  =  0  •  5 
X  I  -0.5 

S 1 GMA (1  )  =  <ETA-1.0)/XBA 
S  I  GMA  (  ?  !  ~  -  c. !  CM  A  (  !  ) 

S IGMa (JJ=-TR*ETA/XBA 
c i GMA ( A  ;  r  —  c ■ G“A ( t ) 
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C  C  TABLE  A-6  (CONTINUED) 

EXX=0.0 

00  70  !  =  1  ,  A 

70  EXX=EXX+S 1 GMA (  I  )*DELTA(  I  ) 

SIGMA  <  5  >  r ( XI- 1 .0 )/YCA 
S 1 GMA ( 6 ) =-Xl/YCA 
SIGMAt  7  )  =  TR* ( D-X 1  l/YCA 
S I GMA (8)=TR*(X1-P)/YCA 
EYY=0.0 
EX Y  =  0  •  0 
DO  80  J=  1 . 4 

EYY=EYY  +  SlGMA(J  +  4>  *DELTA( J  +  4  ) 

80  EXY  =  EXY  +  SI  GMA(  J  +  4  )  *  DELTA  t  J  )  +  S1GMA(  J  )  *DEL.T  A  (  J  +  4  ) 
806  PUNCH  904 , JA I  , J8 I  4 JC 1  4 JD I  • EXX . EYY ,EXY 
GO  TO  801 
C 

700  STOP 
C 

901  FORMAT (  1 3, 4E  16.8 ) 

902  FORMAT (213.E16.8. 13) 

903  F ORM AT (413.E16.8t  13) 

904  FORMAT ( 4  I  3  »  3C 16.0) 

C 

END 
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APPENDIX  B 


THE  MATRIX  FORCE/DISPLACEMENT  PROGRAM 

Due  to  the  analogy  that  exists  between  the  matrix  force 
method  and  the  matrix  displacement  method  of  structural  analysis 
it  has  been  possible  to  write  a  program  that  can  be  used  to 
analyze  structures  by  either  method.  The  analogy  between  the 
two  methods  was  first  shown  by  Argyris  and  Kelsey  (Reference  1) . 
Pestel  and  Leckie  (Reference  5)  have  expanded  on  the  analogy 
and  have  developed  an  excellent  notation  which  will  be  used  here 

The  original  program  was  written  by  the  author  for  a  force 
analysis.  It  was  intended  that  the  experimental  panel  would  be 
analyzed  using  only  the  force  method  and  the  direct  stiffness 
method.  However,  one  of  the  author’s  graduate  students,  Mr. 
Gordon  E.  Lambert.,  elected  to  write  a  displacement  method 
program  for  his  masters  thesis  (Reference  3) .  A  study  of  the 
force  program  and  the  analogy  revealed  that  the  force  program 
could  be  modified  and  used  for  both  types  of  analysis.  The 
modifications  consist  of  one  additional  digit  in  a  control  card 
to  specify  the  type  of  analysis,  some  additional  headings,  and 
one  program  segment  which  is  used  only  in  the  displacement 
analysis . 

The  equations  for  both  methods  are  derived  in  References  1, 
3,  and  5;  therefore,  the  derivation  will  not  be  shown  here.  The 
equations  which  must  be  solved  to  perform  an  analysis  will  be 
stated,  the  analogy  will  be  discussed,  and  some  of  the  unique 
features  of  the  program  will  be  explained. 

The  Force  Method 


lilt.  a.  v_>  I.  v,  ctjua  txuiiD  ait:  . 

!D10!  =  [bi)T(fv'  tv 

(B-l) 

1Dll'  =  [Bl’ 

(B-2) 

[X]  =  -[V'liDi0J 

( B  -  3 ) 

[B]  =  [B01  +  [B1] [Xj 

( B — 4 ) 

B-l 


[Fd] 

=  [b]t[fv] [B] 

(B-5) 

{p} 

=  [B]  {  f  } 

(B-6) 

{ d } 

=  [FdHf} 

(B-7) 

Three  of  the  matrices  in  the  above  equations  must  be  deter¬ 
mined  from  the  physical  and  geometric  properties  of  the  structure: 

[Bq]  =  a  set  of  internal  forces  in  the  structure  due  to 
unit  values  of  the  external  forces.  Each  column 
represents  a  set  of  internal  forces  due  to  one 
of  the  external  loads. 

[B^]  =  a  set  of  internal  forces  in  the  structure  due  to 

unit  values  of  the  redundants.  Each  column  represents 
a  set  of  internal  forces  due  to  one  of  the  redundants. 
For  a  statically  determinate  structure  [B^]  =  [0] 
and  [B]  =  [Bq]  . 

[F  ]  =  the  flexibility  matrix  of  the  unassembled  structure. 
This  is  a  diagonally  partitioned  matrix  where  the 
submatrices  are  the  element  flexibility  matrices. 

Equations  (B-l)  through  (B-7)  can  be  easily  evaluated  once 
these  three  matrices  have  been  determined.  The  [Bg] ,  [B^] ,  and 

[F  ]  matrices  for  the  experimental,  panel  are  shown  later  in  this 
appendix.  The  physical  interpretation  of  the  other  matrices  in 
the  above  equations  are  as  follows: 

[D^g]  =  the  relative  displacements  at  the  redundants  due 

to  unit  values  of  the  external  loads.  Each  column 
represents  the  displacements  due  to  ore  of  the 
external  loads. 

[ D i 1 ]  =  the  relative  displacements  at  the  redundants  due 
to  unit  values  of  the  redundants.  Each  column 
represents  the  displacements  due  to  one  of  the 
redundants . 

[X]  =  the  values  of  the  redundants  due  to  unit  values  of 
the  external  loads.  Each  column  represents  the 
redundants  due  to  one  of  the  external  loads. 
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[B]  =  the  values  of  the  internal  forces  due  to  unit  values 
of  the  external  loads.  Each  column  represents  the 
internal  forces  due  to  one  of  the  external  forces. 
[F^]  =  the  flexibility  matrix,  i.e.,  the  values  of  the 
external  displacements  due  to  unit  values  of  the 
external  loads.  Each  column  represents  the  dis¬ 
placements  due  to  one  of  the  external  loads, 

p  -  the  column  matrix  of  internal  forces, 

f  =  the  column  matrix  of  external  forces, 

d  =  the  column  matrix  of  external  displacements. 

The  Displacement  Method 

The  displacement  equations  are: 


[C1q]  (AjJ  [Kp]  [Aq] 

( B  —  8 ) 

[C11]  "  [Al]T[Kp] [A1] 

(B-9 ) 

[Y]  =  -[C11]"1[C10] 

(B-10) 

[A]  =  [Aq]  +  [A1]  [Y] 

(B-ll) 

[Kf]  =  [A]T[Kd]  [A] 

(B-12) 

mi 

< 

ii 

> 

(B-13) 

{ f }  =  [KJ  d 

i. 

( B-14 ) 

The  three  input  matrices  required  for  a  displacement  analysis 

are : 

[ A q ]  =  the  deformations  due  to  unit  displacements  at  the 
external  loads. 

[A-jJ  =  the  deformations  due  to  unit  displacements  at  the 

kinematic  deficiencies.  A  kinematic  deficiency  exists 
for  each  degree  of  freedom  for  which  there  is  no 
external  load.  If  there  are  loads  corresponding  to 
all  unconstrained  external  displacements,  [A-^]  =  [0] 
and  [A]  =  [Aq] . 

[K  ]  =  the  stiffness  matrix  for  the  unassembled  structure. 

P 

The  physical  interpretation  of  the  other  matrices  are  as 
follows : 
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[C1Q]  =  the  net  forces  at  the  kinematic  deficiencies  due 

to  unit  values  of  the  displacements  at  the  external 
loads.  Each  column  represents  the  forces  due  to 
one  of  the  external  loads. 

[Cn]  =  the  net  forces  at  the  kinematic  deficiencies  due 
to  unit  values  of  the  displacements  at  the 
deficiencies.  Each  column  represents  the  forces 
due  to  a  displacement  at  one  of  the  deficiencies. 

[Y]  =  the  displacements  at  the  kinematic  deficiencies 
due  to  unit  displacements  at  the  external  loads. 

Each  column  represents  the  displacements  at  the 
deficiencies  due  to  a  displacement  at  one  of  the 
external  loads. 

[A]  =  the  deformations  due  to  unit  values  of  the  dis¬ 
placements  at  the  external  loads.  Euv-h  column 
represents  the  deformations  due  to  one  external 
displacement . 

[Kf]  =  the  stiffness  matrix,  i.e.,  the  external  forces 
required  to  produce  unit  displacements  at  these 
forces .  Each  column  represents  the  external  forces 
due  to  a  unit  displacement  at  one  of  the  forces. 

{v}=  the  column  matrix  of  deformations 
{ f }  and  { d >  were  previously  defined. 

Equations  (B-13)  and  (B-14)  give  the  deformations  and  external 
forces  for  prescribed  displacements.  Since  the  external  forces 
are  usually  known,  rather  than  the  displacements,  and  the  internal 
forces  are  needed,  rather  than  the  deformations;  two  additional 
equations  are  necessary  to  fully  implement  a  displacement  analysis. 
They  are 


[B]  =  [K  ]  [A]  (KJ 
[Fd!  =  [Kj]-1 

Once  these  two  equations  are  evaluated,  Equations 
and  (B-7)  can  be  used  to  calculate  the  internal  forces 
displacements . 


(B-15) 

(B-16) 

(B~6) 
and  the 
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The  Analogy  Between  the  Force  and  Displacement  Methods 


The  analogy  between  the  two  methods  is  evident  if  Equations 
(B-l)  through  (B-7)  are  compared  with  Equations  (B-8)  through 


(B-14) .  The  analogous  matrices 


Force 

Displacement 

tf } 

{d} 

<p> 

(v) 

[B0] 

[Aq] 

[Ax3 

[K  ] 

P 

References  1,  3,  and  5  all 
shown  above  with  varying  degree 
explanation. 


are  tabulated 

below 

Force 

Displacement 

O 
i — i 

Q 

O 
i — i 

V 

1 — 1 

1 — 1 
Q 

ICiil 

[X] 

[Y] 

[B] 

[A] 

[Fd] 

[Kf] 

present  tabulations  of  the  type 
of  detailed  operations  and 


The  Matrix  Force/Displacement  Program 


The  Matrix  Force/Displacement  Program  is  divided  into  six 
parts : 


Part  I 
Part  II 
Part  III 
Part  IV 

Part  V 
Part  VI 


Input 

Recursion  Analysis 

Calculation  of  Flexibility  or  Stiffness  Matrix 
Internal  Forces  and  Displacements-Displacement 
Method 

Internal  Forces  and  Displacements-Output 
Initial  Forces 


A  program  listing  is  shown  in  Table  B-l.  Parts  I  through 
V  will  be  discussed  here.  Part  VI,  which  is  used  only  for 
thermal  analyses,  is  discussed  in  Appendix  C.  Although  the  pro¬ 
gram  naturally  divides  into  the  six  parts  listed  above,  the 
primary  reason  for  the  division  is  the  effective  utilization  of 
the  computer.  The  program  is  written  in  PDQ  Fortran  for  the 
IBM  1620  described  in  Appendix  A. 

Part  I  -  Input .  The  inputs  to  the  program  consist  of  a 
control  card  followed  by  the  flexibilities,  or  stiffnesses,  of 
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the  unassembled  structure  and  the  nonzero  elements  of  [Bq]  and 

[B^] ,  or  [Aq]  and  [A^] .  The  control  card  contains  the  basic 
parameters  of  the  structure  such  as  the  number  of  elements; 
internal  forces;  external  forces;  redundants,  or  deficiencies; 
and  nonzero  elements  in  [B^j  and  [B^] ,  or  [Aq]  and  [A^] .  It  also 
contains  a  control  digit  which  specifies  the  type  of  analysis  to 
be  performed.  The  format  of  the  control  card  and  the  other  input 
cards  is  described  in  detail  in  Table  B-l. 

The  input  format  of  the  element  flexibilities,  or  stiffness, 
requires  explanation.  Since  [Fy] ,  or  [K  ] ,  has  an  order  equal 
to  the  number  of  internal  forces,  it  is  usually  very  large.  The 
storage  requirements  are  greatly  reduced  if  tae  element  flexibilities, 
or  stiffness,  are  stored  in  a  special  format  and  used  as  needed. 

The  minimum  storage  requirements  would  result  if  each  flexibility, 
or  stiffness,  matrix  was  given  a  variable  name  and  properly 
dimensioned  ,  or  if  blocks  of  storage  were  set  aside  for  each 
type  of  element-bars,  rectangular  plates,  etc.  This  would, 
however,  reduce  the  generality  of  the  program.  To  avoid  this, 
a  single  array  of  dimensions  £x5  is  used,  where  l  is  tne  number 
of  internal  forces.  Five  columns  are  specified  because,  to 
date,  the  largest  element  flexibility,  or  stiffness,  matrix 
used  has  been  the  5x5  rectangular  or  trapeqoidal  plate  matrix. 

The  element  flexibilities,  or  stiffnesses,  are  right- justified 
in  the  array  to  permit  identification  of  each  type  of  element 
flexibilities,  or  stiffnesses.  The  identification  procedure  is 
described  in  the  next  section. 

The  inputs  are  also  punched  into  cards,  with  appropriate 
headings,  for  problem  identification  and  debugging. 

Far  t  II  -  Recursion  Analysis .  The  solution  of  the  matrix 
force  or  displacement  equations  previously  presented  requires 
the  inversion  of  [D^]  or  [C^]  ♦  Since  the  order  of  [ D ^ ]  is 

equal  to  the  number  of  redundants  and  the  order  of  [C^]  is  equal 

to  the  number  of  kinematic  deficiencies,  the  required  inversion 
may  be  relatively  large.  Festel  and  Leckie  (Reference  5)  have 
presented  a  recurrence  method  where  the  redundants,  or  deficiencies, 
are  divided  into  a  number  of  groups  and  an  equal  number  of  lower- 
ordered  inversions  are  performed.  During  the  programming  of  the 
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recurrence  method  it  was  discovered  that,  not  only  is  the 
programming  simplified,  the  data  storage  requirements  are 
minimized  if  the  number  of  groups  is  equal  to  the  number  of 
redundants,  or  devidicnci.es .  Therefore,  if  there  are  n 
redundants,  or  deficiencies,  there  are  n  one-by-one  inversions 
or,  in  fact,  no  inversions  at  all.  As  is  shown  in  Part  IV,  a 
displacement  analysis  requires  the  inversion  of  the  stiffness 
matrix  [K  ] ;  however,  this  matrix  is  usually  not  large. 

The  recurrence  equations  for  the  force  method  are  as 
follows : 

(D01>  =  1Bi)T'Pv1(B011  (B-17) 

{X1}  =  -  {D^qJ/D^  (B-18) 

[B^1]  =  [bJ]  +  {B^}  { X 1 } T  (3-19) 

Equation  (B-15)  is  the  recurrence  method  equivalent  of  a  com¬ 
bination  of  Equations  (B-l)  and  (B-2) .  These  equations  '"ould 
be  written  as 


[D01]  =  [DlcJDll5  =  [B1]T[Fv][Bo'B1]  (B-20) 

For  each  recursion,  the  last  column  of  the  matrix  [Bq^]  is 
taken  as  { B^ } .  As  a  result,  both  {Dq^}  and  {X1}  are  vectors 

rather  than  rectangular  matrices.  The  unit  redundants  (X1}.are 

determined  by  dividing  all  of  the  elements  of  except  the 

i4*  1 

last,  by  the  negative  of  the  last  element.  The  matrix  [Bq^  ] 

has  one  less  column  than  [B^]  ,  and  after  n  recursions, 

is  the  unit  internal  force  matrix  (Sj  . 

The  recurrence  equations  for  the  displacement  method  can  be 
written  from  Equations  (B-17)  through  (P-19) ,  with  the  aid  of 
the  analogy,  as 

{C01'  =  {Al}T[Kp3  CA013  (B-21) 

{Y1}  =  -  (B-22) 
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(3-23) 


rAoi1]  =  [Ao]  +  tAi}{*1‘}T 

The  triple  matrix  product  in  Equation  (B-17)  and  (B-21) 

lvolves  the  matrices  (F  ]  and  (K  ]  which  are  stored  in  the 

v  p 

ompressed  format  described  in  Part  I.  The  operations  required 

i  T 

or  a  force  analysis  will  be  illustrated.  The  product  { B ^ }  [F  ] 
s  performed  first.  The  first  row  of  [F  ]  is  checked,  from  the 
eft,  for  the  first  nonzero  element.  If  the  first  element  in 
;he  row  is  nonzero,  the  element  f.exibility  matrix  is  a  5x5  matrix. 
Phis  matrix  is  then  premult ipliec.  by  the  first  five  elements  of 
{ B ^ }  to  obtain  the  first  five  elegants  of  the  product.  Then  the 

sixth  row  of  f F  ]  is  checked.  If  the  fourth  element  were  the 

v 

first  nonzero  element,  the  element  flexibility  would  be  a  2x2 
matrix.  The  sixth  and  seventh  element  of  {B^}  would  be  used  to 
obtain  the  next  two  elements  of  the  product.  This  procedure  is 
continued  until  all  elements  of  the  product  have  been  calculated. 
Finally,  [B^]  is  premultiplied  in  the  normal  manner  by  the 

product  (B^)T[Fv]  to  obtain  {D^}. 

Provision  is  also  made  in  this  part  of  the  program  for  punch¬ 
ing  intermediate  output  which  can  later  be  used  for  a  thermal 
analysis  (Appendix  C) . 

Part  III  -  Ca Iculation  of  Flexibility  or  Stiffness  Matrix . 

The  purpose  of  this  part  of  the  program  is  to  evaluate  the  flexi¬ 
bility  matrix,  Equation  ( B  —  5 ) ,  for  a  force  analysis  or  the 
stiffness  matrix,  Equation  (B-12) ,  for  a  displacement  analysis. 

It  should  be  noted  that  these  two  equations  are  not  given  by 
Pestel  and  Leckie.  It  can  be  shown,  however,  that  they  are 
equivalent  to  Pestel  and  Leckie :s  equations. 

A  procedure  similar  to  the  one  used  in  Part  II  for 

{ B^ } T [ Fv3 ,  or  (A^}T[Kp],  is  used  to  obtain  [F^] [B] ,  or  [K^] [A] . 

The  major  difference  is  that  the  products  are  rectangular  matrices 
rather  than  vectors. 

Part  IV  -  Internal  Forces  and  Pis placements -Pis placement 

Method .  This  part  is  used  in  a  displacement  analysis  to  evaluate 

Equations  (B-15)  and  (B-16) .  Since  the  product  [K  ] [A]  was 

P 

determined  in  Part  III  and  stored,  it  is  not  necessary  to  perform 
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this  operation  again  in  tnis  part. 

Part  V  -  Internal  Forces  and  Displacement-Output .  The 
only  purpose  of  this  part  of  the  program  is  to  punch  the  internal 
force  matrix  [B]  and  the  flexibility  matrix  [F^]  for  either  type 
of  analysis. 

Analysis  of  the  Test  Panel 

The  use  of  the  Matrix  Force/Displacement  Program  is  illus¬ 
trated  by  analysing  Panel  Configuration  I  (Figure  A-l)  by  both 
the  force  method  and  the  displacement  method.  The  other  panel 
configurations  could  not  be  analyzed  because  of  the  size 
limitations  of  the  program.  The  results  are  compared  with  those 
calculated  by  the  direct  stiffness  method. 

An  exploded  view  of  the  panel  is  shown  in  Figure  B-l.  The 
independent  internal  forces  selected  for  the  analyses  are  also 
shown.  There  are  twenty  elements,  forty-four  internal  forces, 
and  five  reactions. 

Force  Analysis ♦  The  first  requirement  of  a  force  analysis 
is  to  determine  the  number  of  redundants.  This  number  is  given 
by  the  equation 

R  =  P  +  E  -  D  (B-24 ) 

where  R  =  number  of  redundants. 

P  =  number  of  independent  internal  forces. 

E  =  number  of  external  reactions. 

D  =  number  of  nodal  degrees  of  freedom.  For 

Panel  Configuration  I,  there  are  twelve  nodes, 
each  with  two  degrees  of  freedom, 
therefore  R  =  44  +  5  -  2  x  12  =  25 

The  twenty-five  redundants  must  now  be  chosen.  Rather  than 
attempting  to  define  a  statically  determinate  structure  by 
removing  selected  internal  or  external  constraints,  Argyris ' 
method  (Reference  1)  of  self-equilibrating  redundant  force 
systems  will  be  used.  Figure  C-2  shows  the  four  types  of 
systems  selected  for  this  analysis.  All  possible  combinations 
of  these  systems  must  be  used  to  obtain  the  required  twenty- 
five  redundants.  The  number  of  each  type  of  system  is 


Rib-Web  Systems  =  12 
Stringer-Web  Systems  =  9 
Four  Panel  Systems  =  2 
Two  Panel  Systems  =  _2 
Total  Systems  =  25 


The  unit  redundant  matrix  [B^]  can  be  determined  using  Figures 
B-l  and  B-2.  Each  column  of  [B^]  consists  of  the  internal 
forces  due  to  one  of  the  redundant  systems.  The  seventy-six 
nonzero  elements  of  [B^]  are  given  in  Table  B-2  along  with  their 
row  and  column  designation. 

The  internal  forces  due  to  unit  values  of  the  external 
forces  are  shown  in  Figure  B-3.  Each  of  the  two  systems  of 
forces  is  statically  equivalent  to  its  corresponding  unit  external 
load.  The  nonzero  elements  of  [B^]  are  given  in  Table  B-2  with 
their  row  and  column  designation.  Column  1  corresponds  to  the 
outboard  load  and  column  2  to  the  inboard  load. 

The  element  flexibilities  were  determined  by  inverting  the 
stiffnesses.  The  calculation  of  the  stiffnesses  is  discussed  in 
the  following  section.  The  flexibilities  are  shewn  in  Table  B-2 
in  the  format  in  which  they  are  stored  in  the  computer.  Although 
all  of  the  ribs  and  stringers  are  followed  by  the  plate  elements 
in  this  example,  any  order  could  have  been  used. 

The  results  of  the  analysis  are  also  shown  in  Table  B-2. 

The  two  columns  labeled  "Internal  Forces"  are  the  matrix  [B]  and 
the  "Displacements"  are  the  matrix  [F^] .  These  results  are 
discussed  in  a  later  section. 

Displacement  Analysis .  The  number  of  kinematic  deficiencies 
for  a  displacement  analysis  is  given  by  the  equation 

K  -  D  -  F  -  E  (6-25) 

where  K  =  number  of  kinematic  deficiencies. 

D  =  number  of  nodal  degrees  of  freedom. 

F  =  number  of  external  forces. 

E  -  number  of  external  reactions, 
therefore  K  =  2  x  12  -  2  -  5  =  17 

Each  unconstrained  nodal  displacement  must  be  given  a  unit  value, 
and  the  corresponding  element  deformations  must  be  obtained  to 
determine  the  matrix  [A^j .  Figure  B-4  shows  the  required  dis- 
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placements  for  the  example.  The  unit  displacements  at  the  loads 
(**)  are  also  shown.  The  deformations  resulting  from  these 
displacements  are  used  to  determine  the  matrix  [Aq] .  The 
reference  system  used  to  determine  the  plate  element  deformations 
is  shown  below.  The  X-axis  is  fixed  to  the  lower  edge  of  the 
plate  and  thy  Y-axis  remains  perpendicular  to  the  X-axis. 

This  system  is  consistant  with  the  independent  internal  forces 


(Figure  B-l) .  The  nonzero  elements  of  [A^]  and  [AQ]  are  shown 
in  Table  B-3. 

The  rib  and  stringer  stiffnesses  were  determined  from  the 
equation 

Kg  =  AE/L  (B-26) 

The  plate  element  stiffnesses  were  calculated  from  the  same 
rectangular  plate  stiffness  matrix  used  in  DORA.  The  8x8 
singular  stiffness  matrix  was  reduced  to  a  5x5  nonsingular  matrix 
by  striking  out  the  rows  and  columns  corresponding  to  the  con¬ 
straints  shown  in  the  sketch  above.  The  plate  flexibility  matrices 
used  in  the  force  analysis  were  obtained'  by  inverting  the  stiff¬ 
ness  matrices.  The  element  stiffnesses  are  tabulated  in  Table 
B-3. 

The  results  are  tabulated  in  Table  B-3.  As  in  the  force 
analysis,  the  unit  internal  force  matrix  [B]  is  labeled  "Internal 
Forces"  and  the  flexibility  matrix  [F^l  is  labeled  "Displacements". 

Comparison  of  Results 

A  comparison  of  both  the  "Internal  Forces"  and  the  "Displace¬ 
ments"  calculated  by  the  force  and  displacement  methods  shows 
that  the  maximum  difference  is  one  digit  in  the  seventh  significant 


digit.  This  difference  can  be  considered  to  be  negligible. 

A  comparison  of  the  deflections  at  the  loads  can  also  be 
made  between  the  force/displacement  method  and  the  direct 
stiffness  method.  Using  the  force  method  flexibility  matrix, 
the  displacements  at  the  loads  are 


outboard 

.21905393 

.12551318 

0.5 

=  10  5 

. ^inboard  . 

.12551320 

.1986546 

0.5 

s 


.17228355 

.16198933 


The  above  results  agree  with  and  of  Table  A-2 

through  the  sixth  significant  digit.  Although  the  difference 
here  is  greater  than  the  difference  between  the  force  and  dis¬ 
placement  method  results,  it  is  still  considered  negligible. 

The  Matrix  Force/Displacement  Program  could  be  easily 
modified  to  include  three-dimensional  structures.  The  number 
of  columns  in  [ Fv ]  would  be  increased  and  the  procedures  for  the 
multiplications  involving  [F  ]  would  be  expanded. 

The  major  disadvantage  of  the  present  program  is  that  the 
three  input  matrices  must  be  determined  by  hand.  Preliminary 
investigation  indicates  that,  at  least  for  the  displacement 
method,  all  of  these  matrices  can  be  machine  generated.  This 
possibility  should  receive  additional  investigation. 


t 
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c  c 


TABLE  B-I 


MATRIX  FORCE/D! SPLACEMEMT  PROGRAM  LISTING 
E.L.  COOK  ANO  G.E.  LAMBERT 


~  ww  .  4  /  JO 

THIS  PROGRAM  may  s£  used  to  ANAuYZE  STRUCTUTES  by  EITHFR  the 
matrix  FORCE  method  OR  THE  MATRIX  DISPLACEMENT  method.  THE  structure 
."AY  HAVE  any  COMBINATION  of  DIFFERENT  TYPES  Of  ELEMENTS  AS  long  as  no 
ELEMENT  HAS  A  FLEXIBILITY  OR  STIFFNESS  MATRIX  LARGER  THAN  F I VE-BY-F I VE • 
OTHER  RESTRICTIONS  AR£  — 

I.  MAXIMUM  NUMBER  OF  ELEMENTS  =  55 

a.  maximum  number  of  internal  forces  =  55 
3.  maximum  number  of  external  forces  =  ]0 
A.  MAXIMUM  NUMBER  cf  external  forces 

PLUS  REDUNDANTS  OR  KINEMATIC  DEFICIENCIES  -  45 


THE  INPUTS  TO  THE  PROGRAM  ARE  AS  FOLLOWS- 


*#**  ****  **** 
*  * 

*  card  i  * 

*  * 
*#*■**  ******* 


parameter  caro  -  number  of  elements  <k>,  number  of 

INTERNAL  FORCES  (L).  NUMBER  OF  EXTERNAL  FORCES  (M). 
NUMBER  OF  REDUNDANTS  OR  DEFICIENCIES  (N>.  NUMBER  OF  NON¬ 
ZERO  ELEMENTS  IN  BO  OR  AO  (NZEBO).  NUMBER  OF  NONZERO 
ELEMENTS  in  BI  OR  A!  (NZEBI).  and  A  DESIGNATION  FOR  THE 
METHOD  BEING  USED  <KK>.  FOR  THE  FORCE  METHOD.  KK=I,  AND 
FOR  THE  DISPLACEMENT  METHOD.  KK  =  2.  THESE  PARAMETERS  MUST  BE  IN  FIXED 
POINT  NOTATION  (NO  DECIMAL  POINTS)  AND  MUST  BE  SEPARATED  BY  AT  LEAST 
ONE  BLANK ,  I  .E  .  . 


12  24  2  I  1  7  54  I 


************ 
*  * 

*  CARDS  2  * 

*  ThRu  l+i  * 

*  * 
************ 


ELEMENT  FLEXIBILITIES  OR  STIFFNESSES  -  these  CARDS 
CONTAIN  THE  ELEMENT  FLEXIBILITIES  OR  STIFFNESSES  IN  A 
SPECIAL  FORMAT.  each  CARD  MUST  CONTAIN  FIVE  NUMBERS  IN 
FLOATING  POINT  notation  (DECIMAL  POINTS  REQUIRED).  A 
TWO-FORCE  MEMBER  WILL  HAVE  A  SINGLE  CARD  WITH  FOUR  ZEROS 
FOLLOWED  BY  THE  FLEXIBILITY  OR  STIFFNESS  OF  THE  ELEMENT- 


0.0  0.0  0.0  0.0  .666666675-06 


A  BEAM  ELEMENT ( SAY  pinned  AT  The  ends  with  MOMENTS  at  both  ENDS)  will 
have  TWO  CARDS  WITH  THREE  ZEROS  FOLLOWED  BY  The  ELEMFNT  FLEXIBILITY 
OR  ST IFFNESS.  I .E. . 


0.0  0.0  0.0  .  I  I  I  I  I  I  ! IF-06  -.05555555C-06 
0.0  0.0  0.0  -.0555555NE-O6  .  I  I  I  I  1  I  I  1  E-Of, 

AND  SO  ON  UP  TO  A  PLATE  ELEMENT  WHICH  WILL  HAVE  FIVE  CARDS  CONTAINING 
AN  ELEMENT  FLEXIBILITY  OR  STIFFNFSS  MATPIX  OF  THF  FORM  - 

> 4 0000E -05  .29999E-05  .99999E-D6  -.33333E-06  -.33333E-06 
. 29999E -05  .86666E-05  .56666E-C5  .29999E-05  -.29999E-05 
• 99999E-06  .56666E-05  .66666E-05  .2666&E-95  -.33333E-05 
- • 33333E-0 6  .29999E-05  .26666E-05  .39999E-05  -.I9999E-05 
-.33333E-C6  -.29999E-05  -.33333E-05  -.I9999E-05  .39999E-05 
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TABLE  B-i  (CONTINUED) 


************  NONZERO  ELEMENTS  IN  BO  OP  AO  -  EACH  CF  THESE  C ARCS  CONTAINS 

*  next  *  one  nonzero  element  of  bo  or  ao  preceeded  by  numbers 

*  NZEBO  *  INDICATING  the  ROW  AND  COLUMN  IN  WHICH  IT  APPEARS. 

*  CARDS  *  THE  ROW  AND  COLUMN  ARE  IN  FIXED  POINT  NOT  A  T  1  ON  AND  THE 

************  ELEMENT  is  in  FLOATING  POINT  notation,  1 . e . < 

10  3  -.7071 0680E+01 

************  NONZERO  elements  IN  B1  OR  A1  -  THESE  CARDS  ARE  REQUIRED 

*  NFXT  *  ONLY  FOR  STATICALLY  INDETERMINATE  OR  KINEMATICALLY 

*  NZEBl  *  DEFICIENT  STRUCTURES.  THEY  CONTAIN  THE  NONZERO  ELEMENTS 

*  CARDS  *  OF  Bl  OR  A 1  IN  THE  SAME  FORMAT  AS  THE  CARDS  IN  THE 

»##****#****  PRECEDING  GROUP. 


ci  n  non  non  on 


TABLE  B-l  (CONTINUED) 


C  matrix  force/displacement  program 

C 

C  PART  I  -  INPUT 

C  **■#***##**■#■*•*#**•* 

c 

C  E.L.  COOK  AND  G.E.  LAMBEPT  -  REVISED  APRIL  1966 

C 

COMMON  KtLtMiNi  MN  »KK,FV(55«5>  «  BO  ;  (55*45) 

C 

read  and  punch  parameters 

1  READ  I  0 0 » K » L , M » N » NZEBO , NZEB I  *KK 
GO  TO( 50 1 . 502  > *KK 

501  PUNCH  301 
GO  TO  502 

502  PUNCH  302 

503  PUNCH  300.K.L»M.N»NZEB0,NZE8I 
read  AND  Punch  FV  OR  KP  by  ROWS,  five  VALUES  per  card 

GO  TO( 5 1  I  .512)  .KK 

51 1  PUNCH  31 1 
GO  TO  513 

512  PUNCH  312 

513  DO  2  I  =  1  .  L 

READ  101. FV( I .1 ).FV< I«2)«FV( 1.3) .FV(I ,4)«FV( 1.5) 

2  PUNCH  1 01 ,FV( I , 1 ) ,FV( I .2) .FV( I .3) ,FV( I .4 ) .FV( I .5 ) 

CLEAR  BO  1  OR  AO  I 

MN- M+N 
DO  3  I  -  1  .  L 
DO  3  J= 1 ,MN 

3  BO  1  (  1 . J )=0.0 

REAO  AND  PUNCH  NONZERO  ELEMENTS  IN  BO  OR  AO 

GO  TO ( 52 1 .522) .KK 

521  PUNCH  321 
GO  TO  523 

522  PUNCH  322 

523  DO  4  I  A  = 1  . NZEBO 
READ  1 02,  I  . J.B01  (  I  , J) 

4  PUNCH  1 02. 1 . J.BOl ( I . J) 

IF (N) 14 ,7,530 


B-19 


4 


c  c 


TABLE  B-I  (CONTINUED) 


C 

r 

530 

531 

532 

533 


1  00 
101 
I  02 

300 

301 

30  2 

31  1 
31  2 

321 

322 

331 

332 
C 


read  and  punch  nonzero  elements  in  bi  or  a; 

GO  T0( 531  « 532  > .KK 
PUNCH  331 
GO  TO  533 
PUNCH  332 
DO  6  I  A  = 1  .NZEBI 
READ  102, I , J.Bl I J 
PUNCH  I  02.  I  .  J.B  1  I  J 
JM= J+M 

BO  1  (  I  .JM)=b1  I J 
ALINK-L!NK< I .0) 

STOP 

F ORMAT  *  15) 

FORMAT (52 16. 8) 

FORMAT ( 2 1 5.E 16.8 ) 

FORMAT  <615/1 

FORMAT (26H  ^ORCE  METHOD  -  PARAMETERS/) 

FORMAT < 33H  DISPLACEMENT  METHOD  -  PARAMETERS/) 
FORMAT ( 22H  ELEMENT  FLEXIBILITIES/) 

FORMAT ( 20H  -LEMENT  STIFFNESSES/) 

F0RMAT(/23H  NONZERO  ELEMENTS  in  B0/> 

FORMAT </23H  NONZERO  ELEMENTS  IN  AO/) 

FORMAT ( /23H  NONZERO  ELEMENTS  IN  Bl/) 

FORMAT </23H  NONZERO  ELEMENTS  IN  A  1 /  ) 

FND 
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C  C  TABLE  B- 1  (CONTINUED) 

C  MATRIX  force/displacement  PROGRAM 

C 

C  PART  II  -  RECURSION  ANALYSIS 

Q  *•»#***■******•»******#***#**•**■#•*■* 

C 

C  E.L.  COOK  AND  G.E.  LAMBERT  -  REVISED  JULY  1966 

C 

COMMON  K ,L ,  M, N, MN.KK ,FV( 55  «  5  )  «  BO  1  ( 55*  A 5 ) 
DIMENSION  B1 Tfv( 55 > . DO  1 ( 45 ) 

C 

PRINT  102 
PAUSE 
C 

DO  19  I R  = 1 ,N 
C 

C  B 1  TRANSPOSE  TIMES  FV  OR  A 1  transpose  times  kp 
C 

DO  2  1  =  1  iL 

2  B 1 TF  V  <  I  )=0.0 
J=  I 

00  15  1 S= 1  iK 

1F<FV<  J  «  1  )  ) 20 .5*  3 

3  J4=J+4 

DO  4  I  A= 1 , 5 
1 AJ= 1A+J-1 
DO  4  JA~Ji J4 

4  B 1 TFV (  1 AJ ) =B 1 TFV<  1  A J )+80l  <  JA . MN )  *F V ( J A .  I  A  ) 

J  =  J  +  5 

GO  TO  IS 

5  IF  <FV< J»2)  120.8*6 

6  J3= J+3 

DO  7  I  A  =  2 i 5 
I  A  J  =  1A  +  J-2 
DO  7  JA=J,J3 

7  6  I TFV!  I  A J ) =B 1 TFV (  I  A J ) +BOl  <  J A  « MN ) »F V ( J A ,  1  A ) 

J  =  J  +  4 

GO  TO  15 

8  IF (FV( J i 3)  ) 20  ♦  1  1  <9 

9  J2=J+2 

DO  10  14=3,5 

1 Aj= 1 A+J-3 
DO  10  JA=J,J2 

1 0  B 1 TFV (  1  A J ) =B 1 TF V (  1 AJ  J+BOl  ( JA i MN )  *F V ( J A ,  1  A  ) 

J  =  J  +  3 

GO  TO  15 

11  lF(FV(Jt4)  )20« 14, 12 

1 2  J 1 =J+1 

DO  13  1  A  =  4 , 5 
I  A J  = I A  + J-4 
DO  13  J  A  =  j , j i 

13  B1TFV(IAJ)=B1TFV(1AJ;+B01(JAiMN)*FV(JA»1A) 

J  =  J  +  2 

GO  TO  15 

14  B1TFV(J)=B01 (J»MN!*FV(J<5) 

J  =  J+  1 

15  CONTINUE 
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TABLE  B- I  (CONTINUED) 


C 


C 

C 

C 


C 

C 

C 


C 

C 

C 


C 

C 

C 


C 


CALCULATION  OF  DOKBlfFV  TIMES  BO  1  OR  A1TKP  TiMeS  AO  1  ) 

DO  16  !=I.MN 
DO  I  (  I ) =0.0 
DO  16  J= I *L 

16  DO  1(1) =00 1  <  I  >+81TFV< J)*SOI  (  J  .  I  ) 

punch  intermediate  output  for  initial  strains 

I F ( SENSE  SWITCH  1>21«23 

21  PUNCH  101 .D01 <MN> 

DO  22  I  =  1 »  L 

22  PUNCH  I  01  ,001  (  I ,MN> 

23  CONTINUE 

CALCULATION  OF  REDUNDANTS  OR  DEF I C I ENC I ES ( S TORED  IN  D01> 

DO  17  I = ! . MN 

17  DO  I  (  ! ) = -DO  1(1) /DO  I  ( MN ) 

CALCULATION  OF  NEW  BO  I  OR  A01 

DO  18  I  =  I  .  L 
DO  18  J= I .MN 

18  BO  I  (  I . J ) =801 <  1 . J )+B0 1 (  I .MN) *D01 ( J ) 

19  MN=MN-1 
MN=M+N 

AL I NK  =L I NK (1.0) 

20  STOP 

101  FORMAT ( 5E 1 6 • 8 ) 

102  FORMAT (38H  TURN  SWITCH  I  ON  FOR  THERMAL  ANALYSIS/ I 2H  PRESS  START) 
END 
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C  C  TABLE  B-l  (CONTINUED) 

C  MATRIX  FORCE/DISPLACEMENT  PROGRAM 

C 

C  PART  III  -  CALCULATION  OF  FLEXIBILITY  OR  STIFFNESS  MATRIX 

C  ****•#*■***#***•#*#*•**•**■**■*#**•**■*■*■*■*■**■****«••***■*****#*■****#***** 

C 

C  E  «L  •  COOK  AND  G.E.  LAMBERT  -  REVISED  APRIL  1966 

C 

COMMON  K.L  .M.N.MN.KK ,FV( 55 . 5  )  «B<  55  « 45 ) «FD<  10. 10  ) 

C 

C  CALCULATION  OR  FV  TIMES  B  OR  KP  TIMES  A ( STORED  IN  0) 

C 

MN=M+M 
M 1 =M+ 1 
DO  25  I  =  1  .L 
DO  25  J  =  M 1 »  MN 
25  B ( 1 <J)=0.0 
DO  15  1R=1,M 
I RM= I R+M 
J=1 

DO  15  I S= 1  ,K 
I F ( FV ( J « 1 )  )  20  » 5 «  3 

3  J4=J+4 

DO  4  I A= 1 «  5 
I AJ= IA+J-1 
DO  4  JA=J,J4 

4  B (  I AJ. IRM ) =B(  1 AJ,  IRM  >+B<  JA,  1R)*FV( JA,  1  A ) 

J-J  +  S 

GO  TO  15 

5  I F  <  F V<  J i 2 )  120.8.6 

6  J4  =  J+3 

DO  7  1 A  =  2  «  5 
I AJ= 1A+J-2 
DO  7  JA=J,J4 

7  B<  1AJ,  IRM) =B(  1AJ,  IRM)+B( JA«  1R)*FV( JA.  I  A) 

J  =  J  +  4 

GO  TO  15 

8  I F ( F V ( J  »  3 )  120.1  1  ,9 

9  J4=J+2 

DO  10  1 A=3 . 5 
1  A J= I A  +  J-3 
DO  10  JA=J,J4 

10  B (  1  A J . 1 RM ) =B (  IAJ. 1RM)+3(JA. lR)*FV ( JAt 1A) 

J  =  J  +  3 

GO  TO  15 

11  IF(FV< J.4) )20. 1A. 12 

12  J4=J+1 

DO  13  1 A  =  4 . 5 
1  AJ= 1 A+ J-4 
DO  13  JA=J,J4 

13  B(  I  A- .  IRM) r0(  I AJ. 1RM)+B( JA,  IR)*FV( JA. I  A ) 


GO  TO  15 

1  A  B( J.  IRM)=B( J,  IR)#FV( J.5> 
J  =  J+  1 

15  CONTINUE 
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TABLE  B-l  (CONTINUED) 


C  C 

C  CALCULATION  OF  FD  OR  KP 

C 

DO  18  I  =  I tM 
IM= I+M 
DO  18  J=1,M 
FD<  1  * J ) =0.0 
DO  18  K-i.L 

FO (  1  *  J ) =FD  <  1  »J)+B(K.  1 M  >  *B ( K • J ) 

1 F  <  ABS ( FD  <  l.J)  >-.10**20)17.  17, 16 
FD (  I  *  J ) =0.0 
CONTINUE 
C 

AL 1 NK=L I NK ( 1 .0 ) 

C 

20  STOP 
C 

END 
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T 


non  no 


C  c 


TABLE  B-I  (CONTINUED) 


C  MATRIX  FORCE/DISPLACEMENT  PROGRAM 

C 

C  PART  IV  -  INTERNAL  FORCES  AND  DISPLACEMENTS  -  DISPLACEMENT  METHOD 

C  ***-*•#**#•****■*  •******#■******#*******##****','■•***■&*•*■*#■*•*•*■» •******■•**■*##**•** 

c 

C  E«L.  COOK  AND  G.E.  LAMBERT  -  REVISED  APRIL  1966 

C 

COMMON  K ,L .M.N.MN.KK ,FV( 55 «  5 )  i B ( 55  *  45  >  *  FD<  I  0.  10 ) 

C 

GO  Y0< 1 8.600 ) «KK 
C 

INVERSION  OF  KF  TO  OBTAIN  FD 

600  DO  604  I  =  1  . M 
STORE=FD ( I  ,  I  ) 

FD( I • I  )  =  1 .0 
DO  601  J=  1  .M 

601  FD( I . J) =FD ( I . J ) /STORE 
DO  604  K  -  I  »M 
I F ( K— I  1602.604.602 

602  STORE=FD ( k  ,  I  ) 

FD(K. I ) =0.0 
DO  603  J= 1  ,M 

603  FO(K.J) =FO(K. J)-ST0Re*F0(  I  .J) 

6C4  CONTINUE 

CALCULATION  OF  B 

00  4 1  I = 1 .L 
DO  41  J= 1 .M 
B( I . J) =0.0 
DO  41  K= 1 .M 
I m=K+M 

41  6(  I  «J)=B(  I  .J) +B (  I .  I M ) *FD<  K.J) 

C 

1 8  ALI NKrL I NK ( 1 . 0 ) 

C 

END 
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non  n  n  <1 


c  c  TABLE  B-l  (CONTINUED) 

C  MATRIX  force/displacement  PROGRAM 

c 

C  PART  V  -  INTERNAL  FORCES  AND  DISPLACEMENTS  -  OUTPUT 

C  ****#**#***###*#*#**#♦*#****«*♦♦#**»♦****»**»#*»*♦#), 

c 

C  E.L.  COOK  AND  G.E.  LAMBERT  -  REVISEO  APRIL  !966 

C 

COMMON  K .L  *M,N*  MN.KK ,FV( 55*5  >  *B< 55*45 ) * FD(  I 0< 1 0  > 
INTERNAL  FORCES 
PUNCH  401 

GO  T0(6I 1 .621 ,631  .641 ,651  )  ,M 

611  DO  612  1=1, L 

612  PUNCH  400, B< I , I ) 

GO  TO  700 

62  I  DO  622  I  =  I  . L 
622  PUNCH  400, B< I , I ) ,B( I ,2) 

GO  TO  700 

631  00  632  I = I . L 

632  PUNCH  400, 3< I , 1 ) ,B< I .2) ,B< I ,3) 

GO  TO  700 

641  DO  642  1=1 ,L 

642  PUNCH  400, B< I . 1 ) ,B( I .2) ,B( I ,3) *B( I ,4) 

GO  TO  700 

651  DO  652  I = I ,L 

652  PUNCH  40O,B(  I , I  ).B( I  ,2) ,B<  I ,3),B< I ,4) ,B(  I ,5) 
DISPLACEMENTS 

700  PUNCH  402 

GO  TOC  71  I ,721  ,731 ,74  1  ,751  )  ,M 

711  DO  712  1=1, M 

712  PUNCH  400, FD( I 
GO  TO  800 

72  I  DO  722  I  =  1  , M 
722  PUNCH  400, FD( I 
GO  TO  800 

731  DO  732  1  =  1 , M 

732  PUNCH  400, FD( I 
GO  TO  800 

74 1  DO  742  I  =  1  , M 

742  PUNCH  400, FD< I 
GO  TO  BOO 

751  DO  752  1=1 , M 
7A?  PUNCH  400, "D(l 
C 

800  ALINK=LINK( 1 .0) 

C 

400  FORMAT ( 5E I  6 • B ) 

401  FORMA T  (  /  1  6H  INTERNAL  FORCES/) 

402  FORMAT ( / I 4H  DISPLACEMENTS/) 

C 

END 


,  I  > 

, I ) ,FD< 1,2) 

,  1  ) ,FD<  I  ,2) ,FD<  I *3  > 

,  1  )  ,FD(  I *2)  ,FD<  I *3) ,FD<  I  ,4  ) 

*  I  !  ;FD!  I -2!  >  F  0 !  I  ,3) ,FDC  I ,4 } , F D ‘  i  ,5) 


B-1'6 


(CONT INUEO) 


C  C  TABLE  8- 1 

C  MATRIX  FORCE/Ol SPLACEMENT  PROGRAM 

C 

C  PART  VI  -  INITIAL  FORCES 

c  •**■*****■**#*###*•&####•*•*■***#* 

c 

C  E.L.  COOK  -  JULY  1966 

C 

COMMON  K.L.M.N 
0 1  MENS  ION  BH ( 55 . 55 )  . B I  ( 55  1 
C 

00  10  I  =  I  .  L 
00  10  J=  I  «  L 
10  9H< 1 « J) =0*0 
00  14  1  A= 1  .N 

READ  1 0 1  <0  1  1 
DO  12  J=  1  .L 
12  READ  101 .81 (J) 

00  14  1  =  1 ,L 

00  14  J= 1 .L 

14  BH ( I . J ) =BH ( 1 . J ) -B 1  <  1) *  B 1 ( J ) /D 1  1 
DO  16  1-l.L 
DO  16  J“ 1 «  L 
16  PUNCH  1 0 1 . BH ( I . J ) 

C 

101  FORMAT ( 5E 16*8) 

C 

END 
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c  c 


TABLE  8-2 


FORCE  ANALYS i 

FORCE  METHOD  -  PARAMETERS 

20  4  4  2  25  !  0  76 

ELEMENT  FLEX  1 B ! L 1 T 1 ES 

•00000000E-50  .00000000E-50 

•00000000E-50  • OOOOOOOOE- 50 

.00000000E-50  .OOOOOOOOE-50 
• OOOOOOOOE -50  « OOOOOOOOE -50 

.OOOOOOOOE-50  « 00000000E-50 
.  OOOOOOOOE— 50  • 00000000E-50 
•00000000E-50  .00000000E-50 

• OOOOOOOOE-SO  .P^^00000E-50 
« 00000000E-50  <•  w000000E-50 

•00000000E-50  . wo000000E-50 

•00000000E-50  « OOOOOOOOE- 50 

«  00 00 00 OOF -50  • OOOOOOOOE-SO 
•OOOOOOOOE-SO  .OOOOOOOOE-50 
•OOOOOOOOE -50  . 00000000E-50 
. 1 7438655E-05  . 50000000E-06 

.50000000E-C6  . 1 4887080E-04 

.  1 2438654E-05  . 1 438708 1 E-04 

. 3 1600000E-06  . 48775398E-05 

. 3 1600000E-06  -. 42455399E-05 

•91 097958E-06  . 24999967E-06 

.24999967E-06  . 3748 1 823E-04 

•66097930E-06  • 3723 1 824E-04 
»  3 1 599989E -06  . 699 1 9585E-05 

. 31 599999E-06  - . 63599639E-05 

. 1 7436655E-05  . 50000000E -06 

. 50000000E-06  . 1 4887080E-04 

. 1 2438654E-05  . 1 438708 1 E-04 

•  3 1 600000E -06  . 48775398E-05 

. 3 1600000E-06  - . 42455399E -05 

.9 1 097958E-06  . 24999967E -06 

• 24999967E-06  . 3748 1 823E-04 

.66097930E-06  . 3723 < 824E-04 

. 3 1 599989E-06  . 699 1 Y585E -05 
. 3 1 599999E-06  - .63599639E-05 

. 1 7438655E— 05  . 50000000E-06 
. 50000000E-06  . 1 4887080E-04 

« 1 2438654E-05  . 1 438708 1 E~0 4 
s  3  ! 6O0OOOE  — 06  -/10775390c-  —  C5 
, 3 1 6000C0E-06  - . 42455399E-05 

*  9 1 097958F-06  . 24999967E-06 

.  24999967E-06  , 3748 1 823E-04 
.66097930E-06  . 3723 1 824E-04 
.  3 1 599989E-06  .699 1 9585E-05 

. 3 1 599999E-06  - .63599639E-05 


-  PANEL  CONFIGURATION  1 


•OOOOOOOOE-50  .00000000E-50 
.00000000E-50  . 0000000CE-50 

. OOOOOOOOE- 50  .00000000E-50 

.00000000E-50  .00000000E-50 

.00000000E-50  .0000D000E-50 

.OOOOOOOOE-50  • OOOOOOOOE -50 

.OOOOOOOOE-50  •OOOOOOOOE ~SC 
•OOOOOOOOE-50  .00000000E-50 
.OOOOOOOOE- 50  .00000000E-50 

. 00000000E-50  .00000000E-50 

•00000000E-50  .00000000E-50 

.00000000E-5C  • OOOOOOOOE “50 

.OOOOOOOOE-5C  .00000000E-50 
•00000000E-5C  .OOOOOOOOE-SO 
. 12438654E-05  . 3 1 6OOOOOE-C0 

. 14387081E-04  .48775398E-05 

. 1 5630947E-04  .45615403E-05 

.45615403E-05  . 42807699E-05 

-.4561 5405E-05  -.28076989E-06 

.66097930E-06  . 3 1 599989E-06 

. 37231824E-04  . 699 1 9585E~05 

♦37S92806E-04  . 66759589E “05 

.66759589E-05  . 5668990 1 E~05 

-.66 75964 3E-05  . 2331 0084E~05 

. 12438654E-05  . 3 1 600000E-06 

. 1 433708 1 E-04  • 48775390E-O5 

. 1 5630947E-04  . 456 1 5403E~05 

.456 15403E-05  . 42807699E-05 

-.45615405E— 05  - . 28076989E-06 

. 66097930E-06  • 3 1 599989E-06 

.37231824E-04  . 699 1 9585E-05 

. 37892806E-04  . 66759589E-05 

•66759589E-05  . 5668990 1 E -05 

-.66759643E-05  . 233 1 0084E-05 

. 12438654E-05  . 3 1 6C0000E-06 

. 1 438708 1 E-04  .48775398E-05 

. 1 5630947E-C4  . 456 1 5403E“05 

ACC  1  c  a  A^r  nr  a  nnn-»^  a- 

-.45615405E-05  - . 28076989E-06 

.66097930E-06  . 3 1 599989E -06 

.37231324E-04  , 699 1 9585E-05 

. 37892306E-04  . 66759589E -05 

.66759589E-05  . 5668990 1 E~05 

-.66759643E-05  . 233 1 0084E -05 


. 1 3289036E-05 

•  664  45 1 80E-06 

•  1 32890 36E-05 
. 66445 1 80E-06 
. I  32890 36E -05 
•66445180E-06 
. 1 3289036E-05 
.664451 80E-06 
. 1 32890 36E-05 
. 1 32890 36E— 05 
. 1 32890 36E-C5 
. 1 3289036E -05 
. 1 3289036E-05 
. 13289036E-05 
. 31600000C-06 

-.42455399E-05 
- . 456 1 5405E-05 
-. 28076989C-06 
.42807699E-05 
. 31599999E-06 
-.63599639E-05 
- «  6675964 3E- 05 
. 23310084E-05 
. 56689920 E -05 
. 31600000E-06 
-.42455399E-05 
-. 45615405E-05 
- • 28076989E-06 
• 42807699E-05 
. 31 599999E-06 
- . 6 35996 39 E -05 

—  . 6675964  3E-05 
.23310084E— 05 
•  56689920 E-05 
. 31600000E-06 

-. 42455399E-05 
- . 456 1 5405E-05 

—  *  £ 8 07 6 9 8 9 1 - 0 o 
. 42807699E-05 
. 3 1 599999E-06 

- • 63599639E-05 
- . 6675964 3E -05 
•23310084E-05 
. 56689920 E -05 
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c  c 


TABLE  B-")  (CONTINUED) 


nonzero  elements  IN  BO 


1  0 

1 

•  1 0000000E 

0  1 

12 

1 

•  1 OOOOOOOE 

01 

37 

1 

- • 50000000E 

00 

38 

1 

. 1 OOOOOOOE 

01 

40 

1 

•  50000000E 

00 

41 

1 

• 50000000E 

00 

42 

1 

- • 50000000E 

00 

1  0 

2 

•  1 OOOOOOOE 

0  1 

12 

2 

•  1 OOOOOOOE 

0  1 

1  4 

2 

. 1 OOOOOOOE 

01 

ONZERO  elements  IN  B1 

1 

1 

•  '.OOOOOOOE 

0  1 

IS 

1 

•  1 OOOOOOOE 

01 

2 

2 

•  1 OOOOOOOE 

01 

20 

2 

•  1 OOOOOOOE 

0  1 

3 

3 

•  1 OOOOOOOE 

01 

25 

3 

. 1 OOOOOOOE 

01 

4 

4 

. 1 OOOOOOOE 

01 

30 

4 

. 1 OOOOOOOE 

0  1 

5 

5 

. 1 OOOOOOOE 

01 

35 

5 

• 1 OOOOOOOE 

0  1 

6 

6 

•  1 OOOOOOOE 

01 

40 

6 

« 1 OOOOOOOE 

01 

3 

7 

. 1 OOOOOOOE 

01 

16 

7 

. 1 OOOOOOOE 

01 

1  7 

7 

1 OOOOOOOE 

01 

4 

8 

« 1 OOOOOOOE 

01 

2  1 

8 

. 1 OOOOOOOE 

0  1 

22 

8 

1 OOOOOOOE 

01 

5 

9 

, 1 OOOOOOOE 

0  1 

26 

9 

•  1 OOOOOOOE 

0  1 

27 

9 

1 OOOOOOOE 

01 

6 

1  0 

. 1 OOOOOOOE 

01 

31 

1  0 

•  1 OOOOOOOE 

0  1 

32 

1  0 

1 OOOOOOOE 

0  1 

7 

1  1 

. 1 OOOOOOOE 

0  1 

36 

1  1 

•  1 OOOOOOOE 

0  1 

37 

1  1 

1 OOOOOOOE 

01 

6 

1  2 

. 1 OOOOOOOE 

01 

4  1 

12 

. 1 OOOOOOOE 

01 

42 

1  2 

1 OOOOOOOE 

0  1 

9 

1  3 

. 1 OOOOOOOE 

0  1 

•  O 
l  V/ 

4  ^ 

1  — ' 

—  w  u  ^  o  o  o  e. 

0  i 

1  0 

1  4 

•  1 OOOOOOOE 

0  1 

23 

1  4 

- . 1 OOOOOOOE 

0  1 

1  1 

15 

. 1 OOOOOOOE 

0  1 

28 

15 

1  OOOOOOOE 

0  1 

1  2 

16 

.  1  OOOOOOOE 

01 

33 

16 

1  OOOOOOOE 

Cl 

1  3 

1  7 

.  1  OOOOOOOE 

01 

38 

17 

1  OOOOOOOE 

01 

1  4 

IS 

.  1  OOOOOOOE 

0  1 

43 

18 

1 OOOOOOOE 

0  1 
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TABLE  B-2  (CONTINUED) 


10 

19 

•  1 0000000E 

0  1 

19 

19 

1 OOOOOOOE 

0  1 

12 

20 

•  1 OOOOOOOE 

0  1 

29 

20 

1 OOOOOOOE 

0  1 

14 

21 

*  1 OOOOOOOE 

0  1 

39 

21 

-• ! OOOOOOOE 

0  1 

17 

22 

1 OOOOOOOE 

0  1 

18 

22 

• 20000000E 

0  1 

19 

22 

- • 20000000E 

0  1 

20 

22 

1 OOOOOOOE 

0  1 

21 

22 

. 1 OOOOOOOE 

01 

23 

22 

- • 40000000E 

01 

24 

22 

•40000000E 

0  1 

27 

22 

. 1 OOOOOOOE 

01 

30 

22 

•  1 OOOOOOOE 

0  1 

31 

22 

1 OOOOOOOE 

0  1 

27 

23 

1 OOOOOOOE 

0  1 

28 

23 

•20000000E 

01 

29 

23 

-.2000C000E 

0  1 

30 

23 

1  OOOOOOOE 

01 

31 

23 

« 1 OOOOOOOE 

0  1 

33 

23 

-.40000000E 

01 

34 

23 

.40000000E 

01 

37 

23 

. 1 OOOOOOOE 

01 

40 

23 

«  1  OOOOOOOE 

01 

41 

23 

-.1  OOOOOOOE 

01 

22 

24 

. 1 OOOOOOOE 

0  1 

23 

24 

-.40000000E 

0  1 

24 

24 

c  40000000E 

01 

32 

24 

1 OOOOOOOE 

0  1 

32 

25 

. 1 OOOOOOOE 

01 

33 

25 

—  •  4 OOOOOOOE 

01 

34 

25 

•40000000E 

0  1 

42 

25 

1 OOOOOOOE 

01 

3-30 


internal  forces 

.197227035-01 
.48524270E-01 
-.S7857788E-01 
-.79863338E-01 
-• 169051 67E-01 
,248940308-01 
-.89653643E-01 
-.12518682E  00 
,11 63281 7E  00 
•50709291 E  00 
• 34600593E  00 
*31 432608E  00 
•60246070E  00 
•  1 2306753E  00 
* 19722703E-0 1 
— • 1 0884 1 70E  00 
-.61 164867E-01 
.22368496E  00 
-.24874727E-02 

—•121 48229E  00 

-.33272969E-01 
•56951259E-01 
•60668359E-01 
.9A713037E-01 
. 509839 18E -01 

946381 21E— 01 
-•291 07679E-0 1 
, 24 1 49878E  00 
-.15221922E-01 

- , 329597 70E-0 3 

-•69947840E— 0 1 
•52983390E-0 1 
•86535900E-01 
•26855220E-01 
,  77732955E-0 1 
-•89653643E-01 
-•11659404E  00 
.39753934E  00 
-•15461664E  00 

. l 23400  27E -0 1 
.81 060831E-0 1 
-•87774645E-0 1 
. 31 549070E-01 

• 000 00 00 OE -50 
05Spv_acements 

•2190539 3E -05 
. 1 2551320E-05 


TABLE  B-2 


.  1  3555068E-0 1 
• 37084025E-0 1 
-.6422960 1 E-0  1 
— , 1 0 4 377 24 E  00 
- o64229588E-0 1 

_ • 1 04377 1 9E  00 

. 13555014E-01 
. 37084 1 37E-0 1 
.937068 1 5E-0 1 
.53004589E  00 
.  2  08107  4 1 E  00 
.43336064E  00 

•  9370682 2E -01 

.53004593E  00 

. 1 3555068E-0  1 
-.842391 14E-01 
-.551 06279E-0  1 
. 1 8498399E  00 
.3341 0857E -0 1 
-.1 02261 37E  00 
-. 783774  2  IE -02 
.31 159175E-  01 
•6456685 1 E-0 1 

.93285627E-0 1 

•  200095 1 6E-0 1 

.  20009494E-0 1 
-.20009506E-01 
.  70583396E-0 1 

.1362C602E  00 

.42805897E-0 1 
• 42805S57E-0 1 
- . 42805858E-0 1 
. 58456955E-0 1 
•93285599E— 0 1 
-.842  39085E-0 1 
.  1 355501 4E-0 1 
•  1 2579042E  00 
-.93706022E-01 
.31 2 1 0 1 7 1 E  00 

-.78376081E-02 

- • 1 02261 28E  00 
. 789  39869E-0  1 
.15705242E  00 
.O0OOOOOOE-5O 


. 1 255 1 31 8E-05 

, 19846546E-05 


(CONTINUED* 
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c  C  TABLE 

DISPLACEMENT  ANALYSIS  - 
DISPLACEMENT  METHOD  -  PARAMETERS 
20  44  2  17  5  84 


ELEMENT  STIFFNESSES 


• OOOOOOOOE 

-50 

.00000G00E-50 

c OOOOOOOOE-50 

. OOOOOOOOE- 50 

. 75250000E 

06 

.OOOOOOOOE 

-50 

.00000000E-50 

.OOOOOOOOE-50 

.00000000E-50 

• 15050000E 

07 

.OOOOOOOOE 

-50 

•OOOOOOOOE-50 

.OOOOOOOOE-50 

. 00000000E-50 

. 75250000E 

06 

.OOOOOOOOE 

-50 

•OOOOOOOOE- 50 

.00000000E-50 

• OOOOOOOOE— 50 

•  1 5050000E 

07 

.OOOOOOOOE 

-50 

.00000000E-50 

•00000000E-50 

•  000000 OOF— 50 

. 75250000E 

06 

.OOOOOOOOE 

-50 

.OOOOOOOOE-50 

.00000000E-50 

.00000000E-50 

. 1 5050000E 

07 

.OOOOOOOOE 

-50 

•OOOOOOOOE-50 

. 0000000OE-50 

.OOOOOOOOE— 50 

. 75250000E 

06 

.OOOOOOOOE 

-50 

• OOOOOOOOE-50 

.OOOOOOOOE-50 

.00000000E-50 

. 1 5050000E 

07 

.OOOOOOOOE 

-50 

•OOOOOOOOE-50 

.OOOOOOOOE-50 

.00000000E-50 

. 75250000E 

06 

.OOOOOOOOE 

-50 

.OOOOOOOOE-50 

•OOOOOOOOE-50 

. 00000000E-50 

. 75250000E 

06 

.OOOOOOOOE 

-50 

•OOCCOOOOE— 50 

•OOOOOOOOE-50 

• OOOOOOOOE- 50 

. 75250000E 

06 

.OOOOOOOOE 

-50 

•OOOOOOOOE-50 

.OOOOOOOOE— 50 

. 00  300000E— 50 

. 75250000E 

06 

.OOOOOOOOE 

-50 

.OOOOOOOOE-50 

.000GG000E-50 

. 00000000E-50 

. 75250000E 

06 

. OOOOOOOOE 

-50 

.OOOOOOOOE-5G 

.OOOOOOOOE— 50 

.00000000E-50 

. 75250000E 

06 

.S0394546E 

06 

.30698789E  06 

—  .40  19727  OE  06 

.7221 067 ! E  04 

-. 18274854E 

06 

.30698789E 

C6 

♦ 80394546E  06 

-.70896066E  06 

- • 1 8274854E  06 

. 7  22 1 067 1 E 

04 

-.40 197270E 

06 

-  •  ' . 896066E  06 

. 8o  394546E  06 

-. 722 1 067 1 E  04 

. 18274854E 

06 

. 7  22  1  C67 ' E 

04 

- . i Q274S54E  06 

—  .722  1 067 l E  C  4 

•  4  3844839E  06 

-. 16071504E 

06 

18274854E 

06 

•  722 1 067 1 E  O" 

. 1 8274854E  06 

—  .16  07 1 5C4E  06 

• 43844839E 

06 

. 1 51 29061E 

07 

. 70896070E  06 

- • 756453 1 OE  06 

. 722 1 067 1 E  04 

-. 18274854E 

06 

. 70896070E 

06 

.151 2906 1 E  07 

-. 14654137E  07 

- • 1 8274854E  06 

•  722 1 0671 E 

04 

- . 7564531 OE 

06 

—  •  i 4654 1 37e  07 

. 15 12906  IE  07 

- . 7 22 1 06  7 1 E  04 

. 1 8274354E 

06 

.7221087  IE 

04 

— . 1 0274854E  06 

- . 722 1 067 1 E  04 

. 599 1 634  OE  06 

-.46029672E 

06 

-. 1 8274854E 

06 

. 722 1 067 1 E  04 

. 1 8274854E  06 

-.46029672E  06 

•  599 1 6340E 

06 

•  80  394546E 

06 

. 30698789E  06 

— . 40 1 97270E  06 

. 722 1 0671 E  04 

-. 1 8274854E 

06 

. 30698789E 

06 

. 80394546E  06 

- • 70896066E  06 

-. 18274854E  06 

. 722 1 0671 E 

04 

-.40 197270E 

06 

-.7089606C-E  06 

. 80  394  546E  06 

- • 722 1 067 1 E  04 

. 1 9274854E 

06 

.7221 067 1 E 

04 

- . 1 8274854  E  06 

- . 722 1 067 ! E  04 

. 4 3Q44839E  06 

-. 16071504E 

06 

-. 1 8274854E 

06 

•722I0671E  04 

. 1 827  4854E  06 

-.16071504E  06 

. 4  3844  839E 

06 

. 1 5 1 290b 1 E 

07 

• 70896070E  06 

-t 756453 10E  06 

. 722 1 067 1 E  04 

-. i 82 74 854 E 

06 

.70896070E 

06 

.151 2906 1 E  07 

—  . 146541 37  E  07 

- . 1 8274854E  06 

.  1 22 1 0671 E 

04 

— . 756453 1 OE 

06 

— . 1 4654 1 37E  07 

.151 2906 1 E  07 

-.7221 067 1 E  04 

. 1 8274854E 

06 

•  722 1 0671 E 

04 

- . 1 8274854  E  06 

— .72210671E  04 

. 599 1 6340E  06 

- »  46029672E 

06 

-. 1 8274854E 

06 

.722 1067  IE  04 

. 1 8274854E  06 

-.46029672E  06 

. 599 1 6340c 

06 

.80394546E 

C6 

.30c  3709E  06 

-.401 9727  OE  06 

•  722 1 0671 E  04 

-. 18274854E 

06 

. 306987Q9E 

06 

.80  4546E  06 

- . 70896066E  06 

-• 1 8274654E  06 

•  722 1 0671 E 

04 

-.401 97270E 

06 

-.708  -  j066E  06 

. 80394546E  06 

-. 722 1 067 1 E  04 

. 1 9274854E 

06 

. 722 1 067 1 E 

04 

- . 1 827  4854E  06 

-.7221067  IE  04 

•  4 3844839E  06 

-. 16071504E 

06 

-. 18274354E 

06 

.7221 067 1 E  04 

.18274654E  06 

- . 16071504E  06 

. 43844839E 

06 

. 1 5 1 2906 1 E 

07 

•  7  08960  TOE  06 

- . 756453 1 OE  06 

. 7 22 1 067 1 E  04 

-.182 /4854E 

06 

. 70896070E 

06 

.151 2906 ! E  07 

-.14554137E  07 

- . 1 8274054E  06 

«  7 22 1 0671 E 

04 

- . 7564531 OE 

06 

-.1465413'tE  07 

.151 2906 1 E  07 

- • 722 1 0671 E  04 

. 18274854E 

06 

.7221 0671E 

04 

- . 1 8274854E  06 

-. 722 1 067 1 E  04 

.5991,63405  06 

- . 46029672E 

06 

-. 18274854E 

06 

•  722 1 067 1 E  04 

. 18274854E  06 

- • 46029672E  06 

•  599 1 6340e 

06 
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TABLE  B- 3  (CONTINUED) 


NONZERO  ELEMENTS  IN  AO 


1  3 

38 
14 

39 
43 


1  . 1 OOOOOOOE 

1  .1 OOOOOOOE 

2  • 1 OOOOOOOE 

2  • 1 OOOOOOOE 

2  .1 OOOOOOOE 


NONZERO  elements  IN  A1 


1  1  - 

15  1 

9  2  - 

16  2  - 

17  2  - 

18  2  - 

1  3 

2  3  - 

15  3  - 

16  3  - 

17  3  - 

20  3 

21  4 

22  4 

24  4 

3  5  - 

16  5 

25  5 

9  6 

11  6  - 

18  6 

2  j  6  - 

27  6  - 

28  6  - 

3  7 

4  7  - 

17  7 

21  7 

25  7  - 

26  7  - 

27  7  - 

30  7 

1 0  8 

12  8  - 

19  8 

?7  n 

26  8 

27  8 

29  8  - 

31  8  - 

32  8  - 

33  8  - 

24  9 

31  9 

32  9 

34  9  - 


1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
20000000E 
20000000E 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
40000000E 
40000000E 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
20000000E 
20000000E 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
1 OOOOOOOE 
l 0Qn0009E 
20000000E 
200C0000E 
1 OOOOOOOE 
40000C00E 
40000000E 
1 OOOOOOOE 
1 OOOOOOOE 
40000000E 
4 OOOOOOOE 
1 OOOOOOOE 


0  1 
01 
01 
01 
01 


0  1 
0  1 
01 
01 
0  1 
01 
01 
01 
01 
01 
01 
0  1 
0  1 
01 
0  1 
01 
01 
01 
0  I 
01 
0  1 
0  1 
01 
0  1 
01 
0  1 
0  1 
0  1 
01 
0  1 
0  1 
0  1 
0  1 
01 
0  1 
0  1 
0  1 
0  1 
01 
01 
0  1 
01 
01 
01 
0  1 
O  1 
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TABLE  B-3  (CONTINUED) 


5 

10 

1 OOOOOOOE 

0  1 

26 

10 

. 1 OOOOOOOE 

0  1 

35 

10 

. 1 OOOOOOOE 

0  1 

1  1 

1  1 

•  1 OOOOOOOE 

0  1 

13 

1  1 

1 OOOOOOOE 

0  1 

28 

1  1 

•  1 OOOOOOOE 

01 

36 

1  1 

- • 20000000E 

01 

37 

1  1 

- »  20000G00E 

0  1 

38 

1  1 

1 OOOOOOOE 

01 

5 

12 

•  1 OOOOOOOE 

0  1 

6 

12 

-• 1 OOOOOOOE 

0  1 

27 

1  2 

*  1 OOOOOOOE 

01 

31 

12 

•  1 OOOOOOOE 

0  1 

35 

12 

1 OOOOOOOE 

01 

36 

1? 

-• 1 OOOOOOOE 

0  1 

3"r 

12 

1 OOOOOOOE 

0  1 

AO 

1  2 

*  1 OOOOOOOE 

0  1 

1  2 

13 

. 1 OOOOOOOE 

0  1 

14 

1  3 

1 OOOOOOOE 

0  1 

29 

1  3 

. 1 OOOOOOOE 

01 

33 

13 

» 1 OOOOOOOE 

01 

36 

1  3 

.20000000E 

01 

37 

13 

1 20000000E 

01 

39 

13 

1 OOOOOOOE 

0  1 

41 

1  3 

-.40000000E 

0  1 

42 

1  3 

-•40000000E 

01 

43 

1  3 

1 OOOOOOOE 

0  1 

34 

14 

•  1  OOOOOOOE 

0  1 

4  1 

1  4 

•  40000000E 

0  1 

42 

1  4 

•  4  OOOOOOOE 

0  1 

44 

1  4 

1  OOOOOOOE 

0  1 

7 

15 

1 OOOOOOOE 

0  1 

36 

15 

. 1 OOOOOOOE 

c  1 

7 

16 

. 1 OOOOOOOE 

0  ! 

8 

16 

-• 1 OOOOOOOE 

0  1 

37 

16 

. 1 OOOOOOOE 

0  1 

4  1 

16 

. 1 OOOOOOOE 

01 

44 

1  7 

. 1 OOOOOOOE 

0  J 
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TABLE  B-3  (CONTINUED) 


INTERNAL  FORCES 

. 1 9722669E-0 1 
.48524245E-01 
-.57857780E-01 
- • 79863380E-0 1 
-♦ 1 6905 1 55E-0 1 
.24894 1 1 0E-01 
- . 896536 50E-0 1 
- . 1 25 1 8663E  GO 
. 1 1 632805E  00 
•5Q709297E  00 
.34600586E  00 

•  3 1 4326 1 OE  00 
,6o246065E  00 
. 1 2306755E  00 
. 1972271  IE-01 

10084 170E  00 
-.61 164824E-01 
. 22368487E  00 
— . 24874490E— 02 
-.1214B231E  00 

—  . 33 2 7 30 4 8E  —0 1 
. 5695 1 298E-0 I 

•  60668 326E  —  0 1 
.94713022E— 01 
.50983897E-0I 

-.946381 30E-01 
-.291 07669E-0I 
«  24 1 49880E  00 
-.  1  5222C/30E-0  1 
-.32953000E-03 
— • 69947790E-0 1 
.52983370E-01 
.86535826E-01 
• 26855288E-0 I 
.77732990E-0I 
-.89653670E-01 
-. 1  1  659407E  00 
• 39753937E  00 
- • 15461 66 1 E  00 
. 123401 56E-0 1 
.8 1060900E— 01 
-.87774710E-01 
. 3 1 549050E-01 

—  .  I91.22933E-0"7 

D 1 splacements 

. 2 1 905389E -05 
. 1 255 1 3 1 6E-05 


. 1 3555044E-0 I 
.3708401 5E-0 1 
-.64229596E-01 
-♦ 1 0437723E  OC 
-.64229580E-0 1 

—  • 1 04 3772 1 E  00 

•  1 35551 20E-0 1 
• 37084060E-0 1 
•93706734E-01 
• 53004577E  00 
•208 1 0740E  00 
. 4 3336056E  00 
. 93705820E-0 1 
.53004587E  00 
. 1 3555080E-0 1 

- .84239 1 26E-0 1 
-.551 06275E-0 1 
. 1 8498394E  00 

•  334 1 0889E— 0 1 
-.1 02261 39E  00 
- • 78378230E-0  2 

.31 159214E-01 
•64566807E-0 1 
•93285584E-0 1 
»  20009525E-0 1 
. 200094 70E-01 

—  »  20009491 E-0 1 

•  70583420E-C  1 
. 1 3620587E  00 
»  42806090E-0 1 
•428061 90E-0 1 

-•428061 50E-0 1 
•5845683SE-0 1 
•93285660E-01 
- • 8  4  239030E -0 1 
. I 3555060E-0 1 

•  1 2579024E  00 
-.93706730E-0 1 

•312101 63E  00 
-.78378750E-02 
- • 1 0226 1 48E  00 
•78940080E-0 1 

•  1 5785242E  00 
.  otc  *  1  934E  0 7 


. 1 255 1 3 1 7E -05 
•  1 9846539E-05 
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APPENDIX  C 


RECURRENCE  METHOD  FOR  INITIAL  STRAINS 

Pestel  and  Leckie  have  developed  the  matrix  force  equations 
for  structures  with  initial  strains  in  Reference  5,  Section  9-3. 
They  have  not,  however,  presented  a  recurrence  method  for  solv¬ 
ing  these  equations.  The  following  pages  present  the  derivation 
of  such  a  method  and  an  example  of  its  use. 

The  basic  equations  for  the  internal  forces  and  the  dis¬ 
placements  in  the  absence  of  external  loads  are 

{p}  =  -[B13  [D11]"1lB1jT{h}  (C-l) 

{d}  =  [ B ] T { h }  (C-2 ) 

where  {pi  =  the  internal  forces  due  to  the  initial  strains. 

{d}  =  the  displacements  at  the  external  load  points 
due  to  initial  strains. 

(h)  =  the  deformations  in  the  unassembled  structure 
due  to  the  initial  strains. 

[B^]  =  the  internal  forces  in  the  structure  due  to 
unit  values  of  the  redundants.  This  is  the 
same  as  the  unit  redundant  matrix  used  in  an 
isothermal  analysis. 

[D-^J  =  the  matrix  of  displacements  at  the  redundants 

for  unit  values  of  the  redundants.  This  matrix 
is  also  used  in  the  isothermal  analysis. 

[B]  =  the  internal  forces  due  to  unit  values  of  the 
external  forces.  This  matrix  is  one  of  the 
primary  results  of  an  isothermal  analysis. 

The  only  quantity  in  Equations  (C-l)  and  (C-2)  that  is  not 
available  from  an  isothermal  analysis,  when  the  recurrence  method 
is  used,  is  the  matrix  [ 3  .  This  means  that  [D.^]  must  either 
be  calculated  for  thermal  analyses  or  that  a  recurrence  method 
must  be  derived  which  obviates  the  need  for  [D^]  •  The  equa¬ 
tions  which  can  be  used  to  derive  the  recurrence  equations  are 


C-l 


(C-j) 


{p}  =  [B1Hx} 

{v}  =  { h }  +  [Fv]  { p }  ( C  —  4 ) 

(B^]T{v)  =  10}  (Compatibility)  (C-5) 

where  {x}  =  the  redundants  due  to  the  initial  strains, 

{v}  =  the  deformations  at  the  internal  forces. 

[F  ]  =  the  flexibility  matrix  for  the  unassembled 
structure , 

As  in  the  isothermal  recurrence  method,  the  n  redundants  are 
eliminated  one  at  a  time  with  n  recursions.  Assume  that 
Equation  (C-3)  can  be  rewritten  as 

{p}  -  [Bq] { f 1 }  +  { }  x1  +  [Bp { h } 

For  the  first  recursion 


(p}  [RjHf1}  +  (B*)  J  +  [B*  ]  { h }  (C-6 ) 

where  [bj]  =  first  n-1  columns  of  actual  [B.  ]  . 

u  l 

{ B* }  =  last  column  of  actual  [ B ^ ]  . 

[Bp  =  first  approximation  of  initial  stress  matrix 

[B.  ]  . 

n 

{ f 1 }  =  all  redundants  except  the  fir  t  one  to  be  elimi¬ 
nated  . 

x*  =  the  first  redundant  to  be  eliminated, 
n  =  the  number  of  redundants. 

[ Bp  =  (0]  since  the  structure  is  originally  assumed  to  be 
determinate . 

Substituting  Equation  (C-6)  in  Equation  ( C - 4 ) 

{V1}  =  th}  +  [ F v ] ( [Bq ] { f 1 }  +  IB*}  x1  +  [Bp { h } )  (C-7 ) 

The  relative  displacement  at  the  first  redundant  x*  can  be 
eliminated  by  using  Equation  (C-7)  in  Equation  (C-5)  giving 


where 


TT  ITT  TT  TT 

{Bp  {h}  +  { D^0  >  f  f  1  +  Dpx1  +  fDp>  {h 

(d}o|T  = 

D11  -  (b})TIFv"B1» 


=  0  ( C  —  8  ) 


C-2 


(C-9 ) 


( D  ^  h } T  =  { >T [ Fv)  [B*]  =  {0}  since  [B*]  =  [0] 

Solving  Equation  (C-8)  for  x1 

x1  =  {X1  }T{  f 1 }  -  {  B?  } T {  h } 

D11 

where  { X1  }T  - - -  {D*0)T 

D11 

The  first  redundant  x1  can  now  be  eliminated  from  Equation  (C-6) 
by  using  Equation  (C-9)  . 


{ p>  =  ([B*|  +  {B*HXI}T)  {f1}  -  (—— { B^ }  {  b|  } 1 )  { h) 

D11 


,1.  f  T>  1 1 T , 


(C-10) 


The  force  equation  can  also  be  written  in  terms  of  the  second 
redundant  as  the  unknown,  giving 


{ P }  =  [Bq1]  {  f 11 }  +  { B ^ 1 } x 1 1  +  (B^lth) 
Comparing  Equations  (C-10)  and  (C-ll) 


(C-ll) 


[Bq1]  =  the  first  n-2  columns  of  [Bq]  +  {B^HX1}* 


{B^}  =  the  last  column  of  [Bq]  +  {B^}{XI}T 
IX. 


(B 


D 


~  {B^}{B*}T 


{f11}  =  the“iirst  n-2  elements  of  {  f 1 } 
x11  =  the  last  element  of  {f1} 

The  [ Bq 1 ]  and  { B ^ 1 }  matrices  are  defined  exactly  as  in  a  iso¬ 
thermal  analysis  with  zero  loads. 

The  second  recursion  can  be  continued  by  substituting 
Equation  (C  — 11)  in  liquation  — 4)  . 


{v11}  =  { h }  -f  [Fv]  (  [bJ1]  (f11)  +  { B^ 1  lx1 1  +  [B^Jfh})  (C-12 ) 


Now  substituting  Equation  (C-12)  in  Equation  (C-5) 

{ B^1  }T{ h }  +  (D^l^f11)  +  D“  x11  +  {D^T{h>  =  0 

where  {D**}T=  < B*1 )T [Fv] (B*1 ] 


(C.-13) 


D 

(D 


II 

11 

II 


=  {B*T}T[Fv] {B*1' 


lh>  =  (BI1I)T[Fv][B;T] 


C-3 


The  last  expression  above  can  be  shown  to  be  zero  as 
follows : 


tB“|B“)T(Fv][B“] 


UbJ]  +  {B^}(XI)T)T(FV]  <■ 


T  T  T  1  T  T 

-(B0]  (FV]{B1}  IT"  {BI} 


D 


11 


Ko1  . (b{}t  -(X1)  DJ± 


D 


11 


'11 


=  [0  j  0] 


— - 

Dn 

-  {X1}  {  B^-  >  T  [F  ]  {B*}  ~~{B^T 

U11 

IT  I  IT  I  IT 

{B^}1  =  {  X1  }  {  Bp  1  -{X1HB1}1 


Solving  Equation  (C— 13 )  for  x11 

x11  =  (XII)T(f11. - V-  (B,II)T(h}  (014) 

“ll*1  1 

Substituting  Equation  (C-14)  in  Equation  (C-ll) 

{p}  -  ([bJ1]  +  {BjI}{XII!T){f11}  +  ([B*1] - j — {B^IHBBI}T{h} 

D11 

(C-15 ) 

Writing  the  force  equation  now  in  terms  of  the  next  redundant 
{ p )  =  [Bq11  3  {  f 111  >  +  {B^II}x111  +  [  B  ^ 1 1  ]  {  h  }  (016) 

Comparing  Equations  (C-15)  and  (C-16) 

[Bq11]  -  all  but  the  last  column  of  [ B ^ 1  ]  +  {B^IHXII}T 
{ B J 1 1 }  =  last  column  of  [bJ1]  +  f B^1 > { X1 1 ) T 


C-4 


1 


(B-1,  -  |B«,  - 


-III 


X 


III 


D 


II 

11 


,nII  w_II,T 
■{31  HB1  j 


{fL‘l"L}  =  all  but  the  last  element  of  {t11} 


=  the  last  element  of  {f11} 


The  third  recursion  can  now  be  continued  by  substituting 
Equation  (C-16)  in  Equation  (C-4) 


(v111} 


(h)  +  £Fv]  (  [B^11  ]  {  f 111 }  +  { B^11  ]  x111  f  [BjHi(h}) 


(C-17) 


Substituting  Equation  (C-17)  in  Equation  (C-5) 


{ B  ? 1 1 }T{h}  +  {D 


iii,t/4:iii 


10  ; 


{f- 


III  III 


}  +  D1]L  x 


T TT  t 

+  (Djp'th} 


=  0 


(C-18 ) 


,  /r.III,T  r  DIII  -.T ,  f-III, 

where  {D10  }  =  {B1  }  iFy] [BQ  ] 


D 


III 

11 


<IlT  = 


The  last  equation  above  can  be  shown  to  be  zero  as  follows 

[Bq11 I B^11]  lFv] IBJ11] 

=  (  [Bq11]  +  (B^}(Xn}T)  [Fv]  (  [B^1  ] - ^-{B^HbJ1}1) 

D11 

=  ([Bp]T  +  {XII}{Bp}T)  [FvJ  (Bp]  -  — 

=■  [ B q  1 }  T  [ F v ]  [B^1  ]  +  {X11 }  {  B^1  }T  [Fvl  [bJ1] 


-  [BqI]t[fv3  (b^1}— ^-i-(bi1i}t 


11 

,11,  1 


-  {XII}{B*I}T[Fv] (B1  i — T1 

D11 


{  B ? 1  }T 


C-5 


It  has  been  previously  shown  that 

[bJI|E^I]T[Fv] [B*1]  =  l 0 | 0] 

therefore,  the  first  two  terms  in  the  above  equation 
giving 


(bJi|b”]t[fv1(bJi] 


■  -[b5I)T(fv1(b^>4i  <4 >T 

U11 


-{xiiHb|i}t[fv1  {eJI}-"TI 


{b?i}t 


D 


11 


{n11}  d11 

- — e^1  )T  -  IX11}  (B^)T 


'11 


D 


11 


=  {XII}{B*I}T  -  {XII}{B*I}T  «  r  0  I  0 ] 


Solving  Equation  (C-18)  for  x 


III 


xm  =  {xni}T{fin}  _ 


1  ,-III.T,.  , 
-\B,  }  {h} 


D 


III 1  1 
11 


Substituting  Equation  (C-19)  in  Equation  (C-16) 

(P>  -  UbJ11]  +  IB}IIKX11I)T)l£m)  +  UB”1] 


1  ,0III,T.  ,,  , 

ill  ,B1  )!B1  >  > (h> 


D 


11 


The  recursions  are  continued  until  all  redundan 
eliminated.  For  instance,  if  there  were  only  three 
the  next  force  equation  would  be 

(p)  -  (BFV]t£IV>  +  [Bfnx1V)  +  [bJV!  (h) 

Since  all  redundants  have  been  eliminated 


(£IV)  =  { XIV }  =  0 


are  zero, 


(C-19) 

(C-20) 

to  ci  ITG 

redundants , 

(C-2 1) 

( C  —  2  2 ) 
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Therefore 


{ P  >  =  [B*V){h}  =  ([B*11]  -  { B*11 }  { B*11 }T)  {hr 

D11 


=  ([B*1]  -  -  { B  J 1 }  { B  ^ 1 } T  -  {B^II}{B^II}THh} 


D 


11 


D 


1  1 


=  ( - ~{B*}{B*}T - Yf  {  B^1 }  {  B^1  }T 

°11  °11 


_4_  {BIH}{BIII}T){h} 


~11 

It  can  be  concluded  that 


(C-23) 


{p}  = 


n 

<-  l 


i=l  D 


(bJ}(bJ}T) { h } 


[Bh] {h} 


11 


(C-24 ) 


where 


[Bh]  = 


n 

-I 


i-i 


{Bp{BpT 


(C-25) 


The  means  for  evaluating  Equation  (C-25)  has  been  incorpor¬ 
ated  into  the  Matrix  Force/Displacement  Program  (Appendix  B)  as 
follows : 

1.  In  Part  II  (Recursion  Analysis),  and  { B^ } 

may  be  punched  into  cards  during  each  recursion. 

2.  In  Part  VI  (Initial  Forces),  this  intermediate 
output  is  reread  and  [B^]  is  determined. 

The  procedure  is  illustrated  for  the  simple  bar  structure 
shown  in  Figure  C-l.  The  results  are  shown  in  Table  C-l.  The 
"Internal  Forces"  m  the  table  are  the  unit  load  matrix  [B]  and 
can  be  used  in  Equation  (C-2)  to  obtain  the  thermal  displace¬ 
ments  at  the  load  points.  The  "Displacements"  shown  in  the 
table  are  the  isothermal  flexibility  matrix  [F^]  and  cannot  be 

used  for  thermal  displacements.  The  "Initial  Forces"  are  the 
thermal  unit  loads  [B^]  and  can  be  used  in  Equation  (C-24)  to 

calculate  the  internal  forces  due  to  the  initial  strains. 
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As  a  check  on  the  validity  of  the  recurrance  method,  the 
same  structure  was  also  analyzed  by  the  direct  method.  The 
program  written  for  the  check  is  shown  in  Table  C-3.  This 
program  is  written  in  1620  Fortran  II  and  utilizes  the  IBM 
System/360  Scientific  Subroutine  Package.  The  two  special 
subroutines  NOZERO  and  MATPCH  are  listed.  The  remaining  ones 
were  taken  directly  from  the  package.  The  results  are  shown 
in  Table  C-2 .  All  of  the  forces  and  displacements  agree  to 
six  significant  figures. 

The  test  panel  was  not  used  for  this  comparison  because 
the  dimensions  in  the  Fortran  II  program  could  not  be  made 
large  enough.  The  Force/Displacement  Program  can,  however,  be 
used  for  the  thermal  analysis  of  the  test  panel. 

An  initial  stress  analysis  for  the  displacement  method  has 
not  been  incorporated  into  the  Force/Displacement  Program. 

Pestel  and  Leckie  have  derived  an  equation  for  this  case 
(Reference  5,  p.  314)  and  the  displacement  analogy  of  Equation 
(C-25)  appears  in  their  equation  along  with  ocher  known  quantities. 
Therefore,  the  Force/Displacement  Program  can  be  extended  to 
include  displacement  thermal  analyses. 
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A  =  0.10  in. 

E  =  10  x  ] 06  psi 
L/AE  =  L  x  10~6  in. /lb. 

X1  =  pl'  x2  =  P3 '  and  x3  =  P5 
Figure  C-l.  Structure  for  Initial  Strain  Example 


c  c 


table:  c-  i 


EXAMPLE  -  INITIAL  STRAINS 


MATRIX  FORCE/DISPLACEMENT  PROGRAM 
FORCE  METHOD  -  PARAMETERS 

5  5  2  3  2  9 

ELEMENT  FLEXIBILITIES 


.OOOOOOOOE-50 

.00000000E-50 

•OOOOOOOOE-50 

•OOOOOOOOE-50 

•OOOOOOOOE-50 


.00000000E-50 

•00000000E-50 

•OOOOOOOOE-50 

.ooooooooe-so 

.00000000E-50 


# OOOOOOOOE-50 
.00000000E-50 
.OOOOOOOOE-SO 
.OOOOOOOOE-SO 
.OOOOOOOOE-50 


.00000000E-50 
.OGGOOOOOE-50 
.00000000E-50 
• 00000000E-50 
.OOOOOOOOE-SO 


•  50000 OOO E-0 A 
•aOOOOOOOE-04 

•  5000 00 00 E“0 4 
• 30000000E-04 
.50000000E-0A 


NONZERO 

ELEMENTS  IN  BO 

2 

1 

•  1 OOOOOOOE 

0! 

4 

2 

-.  1 OOOOOOOE 

01 

nonzero 

ELEMENTS  IN  81 

1 

1 

.  1  OOOOOOOE 

01 

2 

1 

-•80000000E 

00 

4 

1 

•60000000E 

00 

2 

2 

— .eopoooooE 

00 

3 

2 

.  1 OOOOOOOE 

0  1 

4 

2 

-•60000000E 

00 

2 

3 

.80000000E 

00 

4 

3 

-«60000000E 

CO 

5 

3 

•  10000000E 

01 

INTERNAL  FORCES 

•  22524931 E  00 

•  4050  3980E  00 
•29320 1 62E  00 
•9437B202E-0I 

- • 2252493 1 E  00 

displacements 

•  1 620 1 592E-04 
-.2831 3459E-05 


•  I 7908264E  00 
- *7078364 1E-0 1 
- »  26968L72E  00 
- • 62328940E  00 
- • J  7908265E  00 


-.2331 3458E-05 
•  1 B698682E-04 


c 


TABLE  C- 


f CONTINUED) 


INITIAL  FORCES 


ROW  1 

-.  14247019E 

05 

*  5631 2329E 

04 

a  1  4549974E 

04 

596942 15E 

04 

-.57529809E 

04 

ROW  2 

.56312329E 

04 

-. 1 4674005E 

05 

•  7  3300406E 

04 

•  2359455 1 E 

04 

- «  563 1 2329E 

04 

ROW  3 

. 1 4549974E 

04 

•  7  3300406E 

04 

1 2072546E 

05 

a  89B95242E 

04 

-.1 4549973E 

04 

ROW  4 

-.5969421 5E 

04 

.2359455  IE 

04 

• 8989S242E 

04 

1 2557020E 

05 

. 596942 13E 

04 

ROW  5 

- . 57529809E 

04 

- •> 5631 2329E 

04 

-. 1 4549973E 

04 

. 596942 13E 

04 

-. 1 4247020E 

05 

c  c 


TABLE  C-2 


EXAMPLE  ~  INITIAL  STRAINS 
GENERAL  MATRIX  FORCE  ANALYSIS  PROGRAM 

INTERNAL.  FORCES 

1  .22524929E+00  » 1 7908262E +00 

2  » 4050 3990 £+00  70783650E-0 1 

3  .29320I58E+00  - . 26968568E+G0 

4  .94378200E-0  1  -.62328950E  +  00 

5  -.22524929E+00  1 7908263E+00 

01 SPLACEMENTS 

1  . I620I 592E-04  - . 283 1 3465E -05 

2  -.2831 3463E-05  . 1 8698683E-04 

INITIAL  FORCES 


1  14247018E+05  . 563 1 2 330E+04  . 1 4549972E+04  - . 59694208E+04 

2  .5631 2330E+04  - . 1 4874 003E+05  » 73300397E+04  . 23594 552E+04 

3  . 1 4549972E+04  , ^3300396E+04  - . 1 2072544E+05  . 8969523 1 E+04 

4  -.59694208E  +  04  . 23594 549E +04  . 89895230E  +  04  - . 1 25570  I 8E+05 

5  -.57529805E+04  - . 563 1 2320E+04  - . 1 4549973E+04  „ 59694208E+04 


-.57529805E+04 
-.56312320E+04 
-. 14549974E+04 
•  59694  209E  +  04 
-. 1 42470 1 8E+05 
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a 


no  oo  oo  oo  oo  oo  oo  oo  on  on  on  on  o  ooooooo 


C  C 


TABLE  C-3 


GENERAL  matrix  force  analysis  program  listing 

C  GENERAL  MATRIX  force  analysis  without  storage  COMPRESSION 

c 

C  E.L.  COOK  -  JULY  1966 

C 

01  MENS  I ONS  FV  (  (LL+L  )  /2 )  .  BO ( LM ) .B 1  (LN) , 61 TFV ( NL ) . 0 1 0 ( NM ) , D 1  I  ( NN )  , 
1LWV(N) , MWV(N)  . X(NM) ,B1X(LM) ,8(LM)  .BTFV(ML) i FD ( MM )«B1D11I(  LN ) , 

26 1 T( NL ) »BH(LL )  ,RwV(L  ) 

NOTE  -  FV  AND  BO  ARE  STORED  IN  BH  ,  THEREFORE*  L  MUST  BE  EOUAL  TO 
OR  GREATER  THAN  1+2M. 

DIMENSION  F V ( 325) *8^(50) »  B  1 ( 125) *B1TFV( 125) »D10( 10) *D1 1 (25) « 
1LWV(5 )  *MWV (5 )  *  X  (  1  0 )  «  Bl X  <  50 ) .B( 50  )  , BTFV( 60 )  .  FD ( 4 ) «  B 1 T (  1 25 )  * 
2810111(125) *  BH ( 625 ) *RWV(25) 

EQUIVALENCE  ( BH (  1  ) .FV(  1  )  )  *  <BH<  352 )  *  BO ( »  )  ) « ( B1 TFV. Bl T ) *  ( 8. BO ) * 

1 (B1D1 1 1 *B1 ) 

READ  NO.  OF  INTERNAL  FORCES.  EXTERNAL  FORCES.  AND  REDUNDANTS 
1  READ  900.L.M.N 

READ  NONZERO  ELEMENTS  IN  FV  (SYMMETRIC  FORM) 

CALL  NO ZERO (FV.L.L. 1 ) 

READ  NONZERO  ELEMENTS  IN  80  (GENERAL  FORM) 

CALL  NO ZERO ( BO  «  L  «  M «  0 ) 

READ  NONZERO  ELEMENTS  IN  Bl  (GENERAL  FORM) 

CALL  NOZERO ( B 1 iL.N.O) 

CALCULATE  Bl  TRANSPOSE  TIMES  FV 
CALL  TPRD (Bl .FV.B1 TFV.L.N.O, 1 «L) 

CALCULATE  DIO  =  31  TRANSPOSE  TIMES  FV  TIMES  BO 
CALL  GMPRD (Bl TFV, BO. D1 O.N.L .M) 

CALCULATE  Dll  -  Bl  TRANSPOSE  TIMES  FV  TIMES  Bl 
CALL  GMPRD (01 TFV.Bl , D1 1 .N.L.N) 

INVERT  D1 1 

CALL  MINV(D1 1 .N.DET.LWV.MWV) 

CALCULATE  -X  =  D1 1  inverse  times  010 
CALL  GMPRD(D1 1 .DlO.X.N.N.M) 

CALCULATE  Bl  TIMES  “X 
CALL  GMPRD (31 ,X ,B1 X .L.N.M) 

CALCULATE  AND  PUNCH  B  =B0  -  BI<-X) 

CALL  GMSUB(B0,S1X->B.L,M) 

CALL  MATPCH( B*L , M, 0 »RWV  ) 

CALCULATE  B  TRANSPOSE  TIMES  FV 
CALL  TPRD(B,FV,QTFV»L.M,0, 1  ,L) 
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TAElE  C-3  (CONTINUED) 


C  C 

C  CALCULATE  AND  PUNCH  FD  =  B  TRANSPOSE  TIMES  FV  TIMES  B 

CALL  GMPRD(BTFV<QiFD<M.L<M) 

CALL  MATPCH(FD«M,M.0,RWV> 

CALCULATE  Qt  TRANSPOSE 
CALL  GMTRAC  . iGIT iLtN) 

CALCULATE  B1  TIMES  Dll  INVERSE 
CALL  TPRD ( B  1  T  ,  D  1  1«B1D111,N.L«0,0,N) 

CALCULATE  “8H  =  B I  TIMES  Dll  INVERSE  TIMES  Bl  TRANSPOSE 
CALL.  GMPRD  (3101  1  1  *31  T,3H,L  .N«L) 

C 

C  CALCULATE  AND  PUNCH  BH  =  -BH (-1.0) 

call  SMPY ( BH» - I •0*8H)LiLi 0 ) 

CALL  MATPCH(BH,L.L.OiRWV) 

C 

GO  TO  I 
C 

900  FORMAT (313) 

C 

ENC 


subroutine  nozeroca,n,m,msi 

C 

DIMENSION  A ( 1 ) 

c 

CALL  SCLA( AiOiOiNiMiMS) 

READ  900, NO 

DO  2  I  A  =  I  «  NO 

READ  900, 1 ,J,Al J 

CALL  LOC ( I ,J, IR , N.M, MS ) 

2  A  (  I  R  )  =  A  I  J 
C 

900  FORMAT ( 2  I  3  *E I  A .8 ) 

C 

RETURN 

C 

END 
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on  ("too 


C 

C 

C 

C 

1  0 

20 
30 

40 

50 

60 

70 

00 

90 

I  00 

1  10 

120 
130 


150 


160 

170 


C 

4Q0 

c 

900 
91  0 
920 


TABLE  C-3  < CONTINUED! 

SUBPOUTINt  MATPCH ( A .N. M.MS. R > 

0 1  MENS  1 CN  A<  1 >  »R< 1  > 

I F ( MS-2 ) 20  *  10.400 

DIAGONAL  FORM 
PUNCH  920 

PUNCH  900 . ( A ( 1 ) . I = l , N ) 

RETURN 

GENERAL  OR  SYMMETRIC  FORM 

I F ( M-5 ) 30  *  30 «  150 

OUTPUT  in  COLUMN  ^ORMAT 

PUNCH  920 

00  130  1 = ) . N 

CALL  LOC! I . 1 . 1R1 .N.M.MS) 

1 F  <  M— 1  >400.6  0.40 

CALL  LOC(  1  .2.  1R2.N,M<.MS> 

I  F  !  M-2  )  400 . 90  .SC 

CALL  LOC ( I . 3. IR3.N. M.MS > 

I  F  (  M'~3  )  4  00 . 100,60 

CALL  LOCI  1  ,4,  1RA.N.M.MS) 

IF (M-4) 400,1 10.70 

CALL  LOC< 1 .5. !R5,N,M,MS> 

GO  TO  120 

PUNCH  910, 1 »A( 1R1 ) 

GO  TO  130 

PUNCH  910. I ,A( IR1 ) ,A( 1R21 
GO  TO  i 30 

PUNCH  910.  I  «A( lRl  > .A(  IR2J  ,A(  1R3) 

GO  TO  I3C 

PUNCH  910,  1  »  A (  I R 1  )  , A (  I R2 ) , A (  IR3),A(  lR4> 

GO  TO  130 

PUNCH  910, 1 <  A ! IR1 ) ,A( IR2) ,A< 1R3) ,  A  ! 1R4) ,A( IR5) 

CONTINUE 

RETURN 

OUTPUT  IN  ROW  FORMAT 

DO  170  1 = 1 , N 

DO  160  J= 1 ,M 

CALL  LOC ( I  J.1R.N.M.MS/ 

R  (  J  )  =  A  (  1  R  ) 

PUNCH  920 

PUNCH  900, (R< J) , J=1 ,M) 

RETURN 


C  T 

I  W* 


FORMAT ( 5X . 5E 1  5  »  8  > 
FORMAT! I5,5ElS®a> 
FORMAT! 1H  ) 


ENO 
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APPENDIX  D 


ANALYTICAL  PREDICTION  OF  THE  TEMPERATURE 
DISTRIBUTION  IN  AN  INTEGRALLY  STIFFENED  PANEL 

D.L.  Hull 

This  report  has  to  do  with  one  phase  of  a  research  project 
which  is  to  conduct  an  analytical  and  experimental  analysis  of 
an  integrally  stiffened  panel  subjected  to  combined  mechanical 
and  thermal  loading.  The  phase  to  be  studied  here  will  be 
the  analytical  prediction  of  the  temperature  distribution  which 
would  result  from  heating  the  edges  of  a  panel  which  is  cooled 
on  the  upper  and  lower  surfaces  by  forced  convection.  The 
panel  to  be  studied  is  shown  in  Figure  D-l. 

NUMERICAL  ANALYSIS 

To  obtain  a  general  idea  of  what  the  effect  would  be  of 
changing  the  geometry  of  the  panel  and  the  values  for  the  co¬ 
efficients  of  thermal  conductivity  and  convective  heat  transfer, 
a  section  was  taken  out  of  the  panel  which  was  one  inch  in 
depth  with  the  configuration  shown  in  Figure  D-2.  This  section 
was  assumed  to  have  no  heat  flowing  from  the  sides  or  the  right 
end  (center  line  of  the  plate)  . 

Making  use  of  the  symmetry  of  the  panel,  one-half  of  the 
section  was  divided  into  26  elements  and  finite  difference 
equations  were  written  for  the  heat  transfer.  These  equations 
were  programmed  for  the  IBM  1620  digital  computer  and  solved 
for  the  temperature  distribution  along  the  length  of  the  section. 
The  program  used  was  one  which  makes  use  of  matrices  for  the 
solution  of  simultaneous  equations.  The  equations  are  of  the 
form  shown  below  for  element  number  one. 


kA 


kA. 

1 


hAi  <T.  -  Ti>  -  nr  (Ti  +  i  -  Ti>  +  if  <Ti  +  3  -  V  + 


Ai  <V  =  0 


where:  A^  is  the  area  between  adjacent  elements  perpendicular 

to  a  line  drawn  between  the  elements 
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L  is  the  distance  between  the  elements 


(i  =  1  to  26)  is  the  temperature  of  element  i 

is  ambient  temperature 

k  is  thermal  conductivity  BTU/hr-°F-in . 
h  is  convective  heat  transfer  coefficient 

BTU/hr-°F-in^ 

substituting  the  geometric  variables  the  equation  becomes 
h  |  (T„-Tx)  -k  <§)  (|)  (T2-Tl)-k(|)  (|)  (T4-Tl)-  |qR  =  0 


Rearranging,  dividing  by  k  and  setting  T  =0  gives 


hK  +  H  3B  H  3B  =_H 

4k  3B  4H  1  3B  2  4H  4  6k  4R 


These  equations  can  be  written  in  the  following  matrix  form: 

[A]  [T]  =  [R] 

where:  A  is  the  temperature  coefficients 

T  is  the  temperature 

R  is  the  applied  heat 

solving  for  T 

[T]  =  [ A]  _1  [ R] 

The  above  equations  were  solved  using  thirteen  different 
sets  of  variables.  One  set  of  variables  was  selected  as  a  base 
and  then  various  variables  were  changed  and  the  results  compared 
with  the  results  from  the  original  set.  This  comparison  gives 
an  indication  of  the  effects  of  varying  plate  geometry  and 
coefficients  of  thermal  conductivity  and  convective  heat  transfer. 
These  effects  are  illustrated  in  Figures  D-3  through  D-15. 

The  variables  used  in  the  thirteen  cases  were  as  follows: 


Set 

I 

Set 

.  II 

Set  III 

H 

= 

.875  in . 

H  = 

.  875 

H 

= 

.875 

t 

= 

.125  in . 

t  = 

.125 

t 

= 

.125 

B 

= 

C  =  .6  in 

1  . 

B  = 

C  =  .6 

B 

= 

C  =  .6 

W 

= 

3.4  in . 

W  = 

3.4 

W 

= 

3.4 

k 

= 

8.0  BTU/hr-F-in . 

k  = 

8.0 

k 

= 

8.0 

h3 

= 

.167  BTU/hr-F-in. 2 

h3 

.167 

h3 

= 

.167 

hl 

= 

W0 

*h2  - 

h4=. 167 

*hi 

= 

h2=h4  = 

hl  " 

0 

Set  IV 

Set  V 

Set  VI 

Set 

VII 

H 

= 

.875 

*H  =  .90 

*H 

=  .750 

H 

= 

.875 

t 

- 

.125 

*t  =  .10 

t 

=  .125 

t 

= 

.125 

B 

= 

C=  .  6 

B  =  C=.  6 

B 

=  C=.6 

B 

= 

C= .  6 

W 

= 

3.4 

W  -  '.4 

W 

=  3.4 

*W 

= 

3.0 

*k 

- 

6.0 

k  =8.0 

k 

=  8.0 

k 

= 

8.0 

h3 

= 

.167 

h3  =  .167 

h3 

=  .167 

h3 

= 

.167 

h. 

i 

= 

h2=h4=0 

hl  =  h2=h4= 

0  h± 

=  h2=h4=o 

hl 

= 

h2=h4  = 

Set  VIII 

Set  IX 

Set  X 

Set  XI 

H 

= 

.875 

H  =  .875 

H 

=  .875 

H 

= 

.875 

t 

= 

.125 

t  =  .125 

t 

=  .125 

t 

= 

.125 

*B 

= 

C= .750 

*B  =  C= . 50 

B 

=  C= . 60 

B 

= 

C= .  60 

W 

= 

3.4 

W  =3.4 

W 

=  3.4 

W 

= 

3.4 

k 

= 

8.0 

k  =8.0 

k 

=  8.0 

k 

= 

8.0 

h3 

— 

.  167 

h2  =  .167 

*h 

3 

=  .050 

*h3 

= 

.083 

hl 

= 

h2=h4=o 

hl  =  h2=h4= 

0  hx 

*"*2-^4-^ 

hi 

= 

h2=h4= 

Set  XII 

Set 

XIII 

H  = 

.875 

H  = 

.875 

t  = 

_  125 

x~  _ 

L.  — 

i  n  r 
•  X  J 

B  = 

C=.  60 

B  = 

C= .  60 

W  = 

3.4 

W  = 

3.4 

k  = 

8.0 

k  = 

8.0 

*h3  " 

.100 

*h3  " 

.200 

hl  = 

h2=h4=0 

hl  = 

h2~‘n4~ 

=  0 

indicates  the  variable  which  has  been  changed. 
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The  program  used  for  the  solution  of  the  abo\/e  sets  of 
equations  made  use  of  matrix  inversion  and  multiplication. 

This  requires  a  large  number  of  storage  locations  in  the  com¬ 
puter  and  a  considerable  amount  of  computer  time.  The  original 
program  written  for  this  purpose  had  to  be  divided  into  four 
separate  programs  so  as  not  to  exceed  the  memory  storage 
capacity  of  the  1620.  The  first  of  these  programs  calculates 
the  constants  to  be  used  in  the  second  and  third  programs  which 
calculate  the  non-zero  elements  of  the  coefficient  matrix  [A] . 
The  fourth  program  inverts  the  matrix  [A]  and  multiplies  it  by 
the  column  matrix  [R]  (the  heat  added  to  the  panel)  to  obtain 
the  temperatures  of  the  twenty-six  elements.  The  approximate 
computer  time  for  the  four  programs  is  shown  below. 


Program 

Input 

Calculation 

Print 

Punch 

No. 

Time 

Time 

Time 

Time 

1 

1 . 5  min . 

.4  min. 

.05  min . 

2 

1 . 5  min . 

.3  min. 

.4  min. 

3 

1 . 5  min . 

.3  min. 

.4  min. 

4 

2 . 5  min . 

15.3  min. 

1 . 1  min . 

7 . 0  min  . 

16.3  min. 

1 . 1  min . 

. 85  min . 

This  gives  a  total  time  fcr  each  run  of  approximately  25.25  min¬ 
utes  for  each  set,  or  a  total  time  of  approximately  5.5  hours. 

The  above  two  dimensional  analysis  gives  a  good  indication 
of  what  to  expect  in  a  three  dimensional  analysis.  Making  use 
of  the  fact  that  the  panel  is  symmetric  about  the  center  lines 
it  was  decided  to  use  only  one-fourth  of  the  panel  in  the  three 
dimensional  analysis.  As  a  first  attempt  at  obtaining  the 
temperature  distribution  for  the  quarter  panel  it  was  decided 
that  the  same  system  would  he  used  that  was  used  for  the  two 
dimensional  analysis.  The  quarter  panel  was  divided  into  364 
elements,  these  elements  were  arranged  in  14  sections  of  two 
rows  each  with  13  elements  in  each  row.  This  created  14  sets 
of  simultaneous  equations  which  can  be  solved  for  temperature. 
The  temperatures  obtained  for  a  given  section  are  dependent  upon 
the  temperatures  of  the  adjacent  sections  therefore,  it  is 
necessary  to  make  repeated  iterations  through  all  14  sections 
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in  order  to  obtain  the  correct  temperature  distribution.  It 
was  found  for  this  system  that  the  temperatures  for  each 
section  had  to  be  punched  out  by  the  computer  then  used  to 
find  14  new  sets  of  column  matrices  {Rj  (the  heat  supplied 
to  each  section)  then  these  column  matrices  were  used  to  solve 
for  a  new  temperature  distribution  in  the  next  iteration.  Due 
to  the  relative  slow  input/output  procedure  on  the  1620  it 
required  approximately  35  minutes  for  a  complete  iteration  of 
the  364  equations  and  it  was  found  after  30  iterations  that  it 
would  take  an  estimated  100  iterations  before  the  temperature 
would  converge  cn  an  equilibrium  set  of  values,  the  time  re¬ 
quired  for  this  would  be  excessive  (approximately  50  hours) 
therefore,  it  was  decided  to  use  another  method  to  solve  for 
the  temperature  distribution. 

For  the  second  attempt  at  solving  for  the  temperature 
distribution  the  quarter  panel  was  divided  into  112  elements 
as  shown  in  Figure  D-16,  and  finite  difference  equations  were 
written  for  these  elements.  Since  the  solution  of  a  large  numbe 
of  simultaneous  equations  by  the  matrix  method  was  found  to  be 
impractical  on  the  available  1620,  the  Gauss-Seidel  Iteration 
method  was  chosen  as  a  means  of  solving  the  equations.  Using 
this  method  and  an  IBM  1620  with  40,000  units  of  storage 
instead  of  20,000  units  it  was  found  that  it  required  apnroxi- 
mately  21  seconds  for  each  iteration.  This  allows  a  large 
number  of  iterations  in  a  relatively  short  time  as  compared  to 
the  method  used  in  the  previous  attempt.  Using  this  system  six 
sets  of  variables  were  chosen  to  investigate  the  effects  of 
changing  panel  geometry  and  the  convective  heat  transfer  co¬ 
efficient.  These  six  sets  are  shown  below.  These  six  sets 
of  variables  required  a  total  of  596  iterations  and  a  total 
time  of  3.6  hours. 


Set  I 

Set 

II 

Set 

III 

Set  IV 

A 

- 

7.70 

*A 

= 

9.70 

A 

= 

7.70 

A 

=  7.70 

B 

= 

3.40 

B 

= 

3.40 

B 

= 

3.40 

*B 

=  3.0 

C 

= 

.30 

C 

= 

.300 

C 

= 

.30 

C 

=  .30 

D 

= 

.60 

D 

= 

.600 

D 

= 

.60 

D 

=  .60 

E 

= 

1.0 

E 

- 

1.0 

E 

= 

1.0 

E 

=  1.0 

t 

= 

.125 

t 

= 

.125 

*t 

= 

.10 

t 

=  .125 

Q 

= 

100. 

Q 

- 

100. 

Q 

= 

100. 

Q 

=  100. 

k 

= 

8.0 

k 

= 

8.0 

k 

= 

8.0 

k 

=  8.0 

h 

= 

.  167 

h 

= 

.167 

h 

= 

.167 

h 

=  .167 

Set  V 

Set 

VI 

A 

= 

7.70 

A 

= 

7.70 

B 

= 

3.40 

B 

= 

3.40 

C 

= 

.30 

C 

= 

.30 

*D 

= 

.750 

D 

= 

.750 

E 

= 

1.0 

E 

= 

1.0 

t 

= 

.125 

t 

= 

.125 

Q 

= 

100. 

Q 

= 

10  0. 

k 

= 

8.0 

k 

= 

8.0 

h 

.167 

h 

— 

.083 

The  results  of  the  three  dimensional  analysis  for  Set  I 
are  shown  pictorially  in  Figure  D-17  and  the  results  for  all 
cases  are  tabulated  in  Table  D-l. 

Using  the  temperatures  obtained  from  the  numerical  analysis 
the  problem  of  solving  for  the  thermal  stresses  in  the  plate 


due  to  these  temperature  distributions  will  be  undertaken.  After 
this  problem  has  been  solved  an  experiment  will  be  set  up  to 
obtain  the  actual  stresses  due  to  a  temperature  gradient  in  the 
‘chanicoi  loads  applied  at  various*  points  on  the 


panel . 
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Figure  3-1.  Panel  Geometry 
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Figure  D-5.  Temperature  Distribution  for  Variable  Set  III. 
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Figure  D-8.  Temperature  Distribution  for  Variable  Set  VI. 
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Figure  D-10.  Temperature  Distribution  for  Variable  Set  VIII. 
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Figure  D-Xl.  jrcr.pejretur^  Distribution  for  variable 
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Figure  D-X2.  T^niparaturA  Distribution  for  Variable  Set  X. 
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TABLE  D-l 


Three-Dimensional  Temperature  Distribution 


SET  I 


Columns 


12  3 


1 

166.73 

137.00 

106.33 

0 

17b. 18 

i^b.bo 

92.72 

3 

lbb.hii. 

126.71 

65.10 

)i 

191.31 

127.26 

62. 70 

5 

193.93 

32b.  06 

62,02 

6 

19b. bo 

12b. Ii9 

62.06 

Rows  7 

I9li  .02 

12b. 99 

62.69 

6 

191.3b 

120.79 

bli .  72 

9 

lb.9.11 

131.32 

90.66 

10 

179.96 

139.61 

101.23 

11 

lb 9. 3b 

131. (|b 

90.72 

12 

192.09 

129.21 

bl.i,  67 

13 

199.29 

129.27 

02.92 

Hi 

19^  .69 

129.96 

62.36 

1 

2 

3 

1 

a. 

170.02 

13q.LiO 

107.96 

p 

161.27 

120.17 

(J1.70 

3 

190.69 

126.1!' 

6!i  .09 

b 

199.91 

lr’r  .6  7 

62.93 

9 

196.96 

130.96 

62.7 

6 

199.3).! 

1 iO.09 

62.76 

Rows  7 

19o.37 

1 301.7)'. 

63.10 

0 

199.33 

1 50.32 

01-.91 

9 

loo.  1), 

132.00 

90.15 

10 

10l.2li 

1)  1.9- 

10"- .  56 

11 

109.32 

132.10 

90.19 

12 

199.61 

130.60 

6b.  62 

13 

199.23 

131. 96 

63.29 

111 

200.73 

131.60 

63.06 

h 

s 

6 

7 

0 

61.16 

61.9b 

60.20 

)i3.79 

60.1-3 

66.70 

53.01 

66.01 

37.55 

33.3b 

60.29 

bn. 20 

6  J  .L-3 

32. 30 

27.52 

56.66 

)  2.65 

36.20 

29.19 

2,6  .62 

59.26 

IlO  .90 

36.29 

27.36 

22.72 

9)'.  96 

bO.  33 

35.66 

26.67 

22.06 

99.60 

60.63 

35.u7 

2n  ,90 

2  ;.  33 

90.06 

h'\  69 

36.65 

26.13 

23.65 

63.9!' 

)h  .02 

30.00 

30.'  9 

26.95 

72.70 

52.96 

60.10 

36.09 

30.37 

63.91 

66.70 

3o.  59 

30.51 

2  b  .60 

97.97 

1-2.35 

3'  .26 

27.65 

23.23 

99. b9 

bO. 23 

3)i 6 

2n.05 

21.57 

5b.  61 

39.39 

33*  °9 

25. 

20.05 

S'bT  II 

Columns 

6 

5 

6 

7 

6 

02.00 

■  2.)-5 

1-9.32 

63.95 

bo .  61 

67.3b 

5 1.6b 

1-5.06 

36.53 

32.21 

50.77 

66.62 

3CV9 

30.51 

25.61 

5l/.70 

60.96 

35.90 

27.13 

20.51 

56 . 66 

3*  .35 

33.93 

25.  36 

20.66 

56  .'ll 

30 .60 

33 .20 

26.76 

20.31 

56  .90 

3W.!  0 

3  *  • '  - 

25.16 

20.69 

97.01 

61.29 

35.2.' 

2  r.  .  66 

22.19 

-  3.06 

'  ‘  n .  0  7 

37.-7 

29.73 

2b. n.  3 

75.-0 

-2.75 

3C  .7'- 

}j>.60 

30.12 

"3.03 

>1 "  .01 

37  . -  '  ' 

20.  r  3 

25.65 

5to.97 

61.10 

36.93 

2o  •  36 

21.9b 

56.69 

39.05 

33.01 

2)i .  61 

20.20 

56  .2o 

36.23 

32.10 

23.02 

19.  bo 
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TABLE  D-l  (Continued) 


ShiT  III 


1  2 

1  168.37  ll.c.oQ 

2  19?  Jil  136.69 

3  205.UU-  133.26 

4  211.55  133.76 

?  211;.  93  131J.71; 

6  216,10  135.32 

Row  7  215.18  135 J;9 

8  211.93  136.05 

9  205,51;  ll.o.  Oli 

10  197.85  152.53 

11  205,95  1),0.20 

12  212.99  I36.-6 

13  216.99  13,  .33 

Hi.  218.83  136.85 


Columns 

3  4 

llli.71;  8-:. 67 
95.97  69.80 
85.79  58.83 
82.52  5U. 29 

81.68  52.55 

81.76  52.25 

82.68  53.22 

85. 1 1.6  56.37 

93.77  61-  .3); 
108.8);  76.6I 

93*86  6)1.31 

85.69  56. 2 u 

83.02  53.05 

82.27  51.92 


5  o 

65.01;  50.  Id; 
53*21  86, 89 

ldi.59  8-1.08 
)i  0 . 26  37.I6 

38.16  3)-. 76 

37.51  33.74 
38.15  33.98 

I' 0.52  35.39 
45.95  37.76 
53.76  35.8.1 

8-5.60  37.53 

1.0.12  38.67 

37.88  32.71 

36.36  31.72 


7  8 


44' 

.91 

81 

.47 

3o, 

.96 

32 

.85 

30, 

.60 

25 

•87 

26, 

.91 

21 

.75 

24, 

.87 

19 

.82 

24, 

.05 

19 

.10 

28. 

•  33 

19 

1  - 
.4J) 

25. 

■  79 

21 

.00 

28. 

,92 

28 

» 56 

33. 

27 

29 

.46 

28. 

71 

28, 

.39 

25. 

23 

20, 

.53 

23. 

34 

lo, 

.59 

22. 

88 

17, 

.77 

SiT  IV 


2 

3 

8. 

5 

6 

Row  7 
8 

9 

10 

11 

12 

13 

18- 


170.72 

160.01 

188.07 

192.79 

195.32 
196.13 

195.33 
192.71 
187.50 
181.02 
167.78 

193.33 
196.8.2 
197.82 


2 

182.8  5 
133.62 
132.28 
13  3*  30 

133.86 
138-.26 
138.23 
13);. 25 
136.33 
11l)i.80 

136.87 

13)i.61 

138.85 

135.21 


Columns 


3  8 

ID-. 28  91. n9 

100.63  76.01 

93.63  69.67 

91.51  66.23 
90.99  nit  .08 
91.02  68,59 

91*52  65,28 

93.28  67.51 

98,88  7^#l.  I. 

109.65  82  .*88 
98.92  73.80 
93.39  67.39 

91.70  65.02 

91.26  68.20 


5  6 


72 

.80 

59.  o2 

'  3' 

.01 

65.93 

55. 

.01 

50.86 

52, 

.13 

1-7.28 

50, 

.32 

45.06 

89. 

>78 

)•»<  .18 

5o. 

,25 

88  ..47 

52. 

20 

45.87 

?L 

w  • 

Ln 

V_*  f 

1.  O  1  r\ 
LLV  #  i.  ^ 

63. 

15 

89.77 

56. 

63 

47.97 

51. 

84 

45.28 

89. 

61 

83. Ill 

80. 

70 

82.89 

7  8 

6)1.32  51.17 
85.96  82.70 

)'0.88  j$6.V7 

37.29  32.28 

35.29  30.34 

38.52  29.61 

за. 88  29.97 

зб. 33  31.58 

39.44  35.14 

83.62  39.97 

39.25  38.97 

35.62  31.11 
33.92  29.13 
33.06  20.27 
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TABLE  D-l  (Continued) 


Row 


Row 


SET  V 


Columns 


* 

2 

3 

4 

5 

6 

7 

6 

1 

176.77 

Hi  1.19 

108.20 

82.46 

63.02 

Ij6.6l 

42 . 0).| 

39.20 

2 

18)1.63 

130.  h2 

94.06 

69.77 

53.86 

45.66 

36.61 

32.74 

3 

192.1b 

12:. .96 

85.73 

60.93 

1-6.93 

41.93 

31.64 

27.21 

4 

196.69 

126.85 

62.79 

57.09 

43.36 

36.96 

20.69 

21.16 

5 

199.16 

127.37 

61.92 

59 . 95 

61.'  1 

37.06 

27.21 

22 .  -4 

6 

200.07 

127.79 

81.93 

99.25 

41.04 

36.22 

26.50 

21.69 

n 

i 

129.48 

128.06 

82 ,6C 

96.03 

)il. 23 

9-..30 

26.65 

22.10 

8 

197.23 

120.61 

65.04 

58.56 

43.39 

37.23 

27.73 

23.31 

9 

:  92.74 

132.33 

91.52 

6lj  .71 

47.61 

30.61 

30.11 

2c.  07 

10 

167.30 

141.85 

102.24 

73.72 

53.20 

39.02 

33.39 

29.77 

11 

1°3«05 

132.21 

91.5? 

64.67 

47.1.i6 

36.55 

29.90 

25.69 

12 

196.06 

129.27 

65.19 

58.42 

43.01 

30.56 

27.19 

22.63 

13 

200.89 

126.63 

82.93 

22.81 

40.85 

39.12 

25.70 

21 ,26 

14 

202.19 

128.91 

82.23 

54*66 

39.97 

3^.39 

25.02 

20.59 

1 

2 

3 

SET  VI 

Columns 

U 

r* 

6 

7 

6 

1 

242.67 

210.20 

17^.93 

l)i7.09 

124.0)1 

108.22 

1^1.20 

96.75 

2 

222.75 

203.90 

163.96 

134.94 

114.55 

104.6  5 

93.24 

87.42 

3 

261.81 

203.01 

156.92 

12.  .12 

106.76 

96.95 

8c  .05 

79.17 

4 

267.24 

204  .Oli 

124.41.1 

122.01 

102.41 

94. 8k 

81.53 

74.29 

5 

270.20 

202.07 

153.70 

120.21 

100.09 

92.19 

76.66 

71.56 

6 

271.20 

202.28 

153.70 

119.77 

99.25 

91.02 

77.76 

70.45 

7 

270.38 

202.60 

156.35 

120.54 

99.77 

91.19 

76.03 

70.62 

Q 

U 

267.51 

crU^.OU 

15<  .37 

123.13 

101.91 

92.56 

79.76 

72.84 

9 

261.62 

207.29 

161.95 

129.24 

10.  .55 

94  .91 

83.37 

77.18 

10 

254.96 

213.69 

171.32 

137.68 

112.74 

96.43 

07.97 

62.62 

11 

262.14 

207.48 

161.99 

129.16 

106.30 

94.54 

63.01 

76.64 

12 

268.33 

206.08 

15C  .1*5 

1^2.85 

101.25 

91.59 

78.82 

71.93 

13 

271.80 

206.62 

154.20 

120.04 

98.63 

69.47 

76.41 

69.24 

14 

273.38 

206,82 

153. 90 

116.98 

97.50 

88. ho 

75.28 

66.05 
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